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Abstract 

The concept of entangled quantum states is considered in the context of 
systems of identical particles, based on the requirement that in order to repre- 
sent physical states both for the overall system and the sub-systems which may 
be entangled, the density operators must satisfy the symmetrisation principle 
and conform to super-selection rules that prohibit states in which there are co- 
herences between differing total particle numbers. These requirements, their 
justification and recent challenges to them are fully discussed in the paper. In 
the second quantisation approach used in this paper, both the system and the 
sub-systems are modes (or sets of modes) rather than particles, particles being 
associated with different occupancies of the modes. The definition of entan- 
gled states is based on first defining the non-entangled states - after specifying 
which modes constitute the sub-systems that may or may not be entangled, 
entanglement being a relative concept. Although multimode systems are also 
considered, this paper is mainly focused on two mode entanglement for massive 
bosons, such as bosonic atoms, but the results also apply for photons. A number 
of inequalities involving variances and mean values of operators based on the 
mode annihilation and creation operators have been proposed as tests for two 
mode entangled states, including the inequalities that define spin squeezing - 
the spin operators are bilinear combinations of the mode operators. The crite- 
ria for spin squeezing in two mode systems is examined, and it is found that 
spin squeezing is best considered for principle spin operator components, where 
the covariance matrix defining the spin fluctuations is diagonal. The modes 
related to the principle spin components are then orthonormal linear combi- 
nations of the original ones. It is shown that the presence of spin squeezing 
in at least one of the spin components requires entanglement of the relevant 
pair of modes. The paper also considers several of the other proposed tests for 
entanglement, including ones based on the sum of the variances for two spin 
components. However, whilst all of the tests are still valid when the present 
concept of entanglement based on the symmetrisation and super-selection rule 
criteria is applied, additional tests have been obtained here. Sometimes the new 
tests are less stringent than those obtained in other papers. In one case where 
the link between spin squeezing and entanglement was based on a concept of 
entanglement which is inconsistent with the symmetrisation principle, a revised 
treatment with modes rather than identical particles as the sub-systems does 
lead to spin squeezing being a test of entanglement if the sub-systems are pairs 
of modes - each pair however only being associated with a single boson. 
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1 Introduction 



Since the work of Einstein et al [T] on local realism, the famous cat paradox of 
Schrodinger [2] and the derivation by Bell et al [3J, [1] of inequalities based on 
treating measured quantities via a classical hidden variable theory which certain 
quantum quantum systems violated, entanglement has been recognised as being 
one of the key features that distinguishes quantum physics from classical physics. 
It arises in the context of quantum systems composed of distinct components 
or sub-systems. Such sub-systems are usually associated with sub-sets of the 
physical quantities applying to the overall system, and in general more than one 
choice of sub-systems can be made. The formalism of quantum theory treats 
pure states for systems made up of two or more distinct sub-systems via tensor 
products of sub-system states, and since these product states exist in a Hilbert 
space, it follows that linear combinations of such products could also represent 
possible pure quantum states for the system. Such quantum superpositions 
which cannot be expressed as a single product of sub-system states are known 
as entangled (or non-separable) states. The concept of entanglement can also 
be extended to mixed states, where quantum states for the system and the sub- 
systems are specified by density operators. The detailed definition of entangled 
states is set out in Section [2] This definition is based on first carefully defining 
the non-entangled (or separable) states such that all non-entangled states must 
be possible physical states. The entangled states are then just the physical states 
which are not non-entangled states. The notion of physical states, the nature 
of the systems and sub-systems involved and the specific features required in 
the definition of non-entangled states when identical particles are involved is 
discussed in detail in Section [2] Entangled states underlie a number of effects 
that cannot be interpreted in terms of classical physics, including spin squeezing, 
non-local measurement correlations - such as for the Einstein-Podolski-Rosen 
(EPR) paradox and violations of Bell Inequalities. Recent expositions on the 
role of entanglement in these effects and more generally in quantum information 
science include [5], [BJ, [7], 0, [9J, [ID] . 

It is now generally recognised that entanglement is a relative concept 
[12]) [33] )j [2], 03] > [H] and not only depends on the quantum state under dis- 
cussion but also on which sub-systems are being considered as entangled (or 
non-entangled) . A quantum state may be entangled for one choice of the sub- 
systems but may be non-entangled if another choice of sub-systems is made, 
an example being for the hydrogen atom [13] where energy eigenstates are non- 
entangled if the sub-systems refer to the centre of mass and the relative position 
of the electron from the proton, but which would be entangled if the sub-systems 
were the positions of the electron and proton. An example involving two differ- 
ent choices of single particle states in a two mode Bose condensate is given in 
Section [4] 

For a general quantum state various measures of entanglement have been 
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defined - see [7J , [5] , [TJ] , [15] , [TB] , [T7J , [IB] , for details of these, and are aimed 
at quantifying entanglement to determine which states are more entangled than 
others. This is important since entanglement is considered as a resource needed 
in various quantum technologies, such as teleportation and quantum information 
processing. Calculations based on such measures of entanglement confirm that 
for some choices of sub-systems the quantum state is entangled, for others it 
is non-entangled. For two mode pure states the entanglement entropy - being 
the difference between the entropy for the pure state (zero) and that associated 
with the reduced density operator for either of the two sub-systems - is a useful 
entanglement measure. As entropy and information changes are directly linked 
[7], [5], this measure is of importance to quantum information science. Another 
entanglement measure is particle entanglement entropy, defined by Wiseman 
et al [19], [20], [16] for identical particle systems and based on projecting the 
quantum state onto states with definite particle numbers. 

Although not directly relatable to the various quantitative measures of en- 
tanglement, the results for certain measurements can play the role of being 
signatures or witnesses or tests of entanglement [14) . |15j . |16) . These are in the 
form of inequalities for variances and mean values for certain physical quanti- 
ties, which are consequent on the inequalities that would apply for non-entangled 
quantum states. If such inequalities are violated then it can be concluded that 
the state is entangled for the relevant sub-systems. It cannot be emphasised 
enough that these tests provide sufficiency conditions for establishing that a 
state is entangled. The failure of a test does not guarantee non-entanglement 
- sufficiency does not imply necessity. The violation of a Bell inequality is an 
example of such a signature of entanglement, and the demonstration of spin 
squeezing is regarded as another. However, the absence of spin squeezing (for 
example) does not guarantee non-entanglement, as the case of the NOON state 
in SubSection 14.31 shows. A significant number of such inequalities have now 
been proposed and such signatures of entanglement are the primary focus of the 
present paper, which is aimed at identifying which of these inequalities really 
do identify entangled states, especially in the context of two mode systems of 
identical bosons. 

At present there is no clear linkage between quantitative measures of entan- 
glement (such as entanglement entropy) and the quantities used in conjunction 
with the various entanglement tests (such as the relative spin fluctuation in 
spin squeezing experiments). Results for experiments demonstrating such non- 
classical effects cannot yet be used to say much more than the state is entangled, 
whereas ideally these experiments should determine how entangled the state is. 

This paper deals with identical particles - bosons or fermions. In the second 
quantisation approach used here the system is regarded as a set of quantum 
fields, each of which may be considered as a collection of single particle states or 
modes. Hence both the system and sub-systems will be specified via the modes 
that are involved, so here the sub-systems in terms of which non-entangled (and 
hence entangled) states are defined are modes or sets of modes, not particles 
[IT] . [i~2] . [13]), [7J, [21], [22]. In this approach, particles will be described via 
the occupancies of the various modes, so that situations with differing numbers 
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of particles will be treated as differing quantum states of the same system, not 
as different systems - as in the first quantisation approach. Note that the choice 
of modes is not unique - original sets of orthogonal one particle states (modes) 
may be replaced by other orthogonal sets. An example is given in Section (3) 
Modes can often be categorised as localised modes, where the corresponding 
single particle wavefunction is confined to a restricted spatial region, or may 
be categorised as delocalised modes, where the opposite applies. Single parti- 
cle harmonic oscillator states are an example of localised modes, momentum 
states are an example of delocalised modes. This distinction is significant when 
phenomena such as EPR violations and teleportation are considered. 

Although multi-mode systems are also considered, in this paper we mainly 
focus on two mode systems of identical bosonic atoms, where the atoms at most 
occupy only two single particle states or modes. For bosonic atoms this situa- 
tion applies in two mode interferometry, where if a single hyperfine component 
is involved the modes concerned may be two distinct spatial modes, such as in a 
double well magnetic or optical trap, or if two hyperfine components are involved 
in a single well trap each component has its own spatial mode. Large numbers 
of bosons may be involved since there is no restriction on the number of bosons 
that can occupy a bosonic mode. For fermionic atoms each hyperfine compo- 
nent again has its own spatial mode. However, if large numbers of fermionic 
atoms are involved then as the Pauli exclusion principle only allows each mode 
to accommodate one fermion, it follows that a large number of modes must 
considered and two mode systems would be restricted to at most two fermions. 
Consideration of multi-mode entanglement for large numbers of fermions is out- 
side the scope of the present paper (see [23] for a treatment of this) , and unless 
otherwise indicated the focus will be on bosonic modes. Although oriented to- 
wards identical bosonic atoms the paper also applies to photons, though details 
of the analysis will differ. 

The work presented here begins with the fundamental issue of how an en- 
tangled state should be defined in the context of systems involving identical 
particles. In the commonly used mathematical approach for defining entangled 
states, this requires first defining a general non-entangled state, all other states 
therefore being entangled. It is contended that the density operators both for 
the overall system states and for the sub-system states of non-entangled states 
must represent physical states and in some other work (discussed below) this is 
not the case. A key feature required of all physical states for systems involving 
identical particles, entangled or not is that they satisfy the symmetrization prin- 
ciple. This places restrictions both on the form of the overall density operator 
and also on what can be validly considered to be a sub-system. In particular 
this rules out individual identical particles being treated as sub-systems, as is 
done in some papers (see below). In addition, super- selection rules (SSR) [24j 
only allow density operators which have zero coherences between states with 
differing total numbers of particles to represent valid physical states, and this 
will be taken into account for all physical states of the overall system, entangled 
or not. This is referred to as the global particle number super- selection rule In 
non-entangled or separable states the density operator is a sum over products 
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of sub-system density operators, each product being weighted by its probability 
of ocurring (see below for details). For the non-entangled or separable states, 
a so-called local particle number super-selection rule will also be applied to the 
density operators describing each of the sub-systems. These sub-system density 
operators must then have have zero coherences between states with differing 
numbers of sub-system particles. This additional restriction excludes density 
operators as defining non-entangled states when the sub-system density opera- 
tors do not conform to the local particle number super-selection rule. Conse- 
quently, density operators where the local particle number SSR does not apply 
would be regarded as entangled states. This viewpoint is discussed in papers 
by Bartlett et al [25], [2(5] as one of several operational approaches for defining 
entangled states. However, other authors such as [27], [28] state on the contrary 
that states when the sub-system density operators do not conform to the local 
particle number super-selection rule are still separable, others such as [25], [50] 
do so by implication, so in this paper we are advocating a revision to the widely 
held notion of entanglement in identical particle systems, .the consequence being 
that the set of entangled states is now much larger. This is a key idea in this 
paper - not only should super-selection rules on particle numbers be applied to 
the the overall physical state, entangled or not, but it also should be applied 
to the density operators that describe states of the modal sub-systems involved 
in the general definition of non-entangled states. The reasons for adopting this 
viewpoint are set out below. Apart from the papers by Bartlett et al [25J, [25] 
we are not aware that this definition of non-entangled states has been invoked 
previously, indeed the opposite approach has been proposed [27], [28]. However, 
the idea of considering sub-system states which satisfy the local particle num- 
ber SSR has been presented in several papers - [27], [28], [25], [26], [31], [32] . 
[3"3"] . mainly in the context of pure states for bosonic systems, though in these 
papers the focus is on issues other than the definition of entanglement - such 
as quantum communication protocols [27], multicopy distillation [25], mechan- 
ical work and accessible entanglement [3T], [32] and Bell inequality violation 
[33j . The consequences for entanglement of applying this super-selection rule 
requirement to the sub-system density operators are quite significant, and in 
the present paper important new entanglement tests are determined. Not only 
can it immediately be established that spin squeezing requires entangled states, 
but though several of the other inequalities (see below) that have been used as 
signatures of entanglement are still valid, additional tests can be obtained which 
only apply to entangled states that are defined to conform to the symmetrisation 
principle and the super-selection rules. 

It is worth emphasising that requiring the sub-system density operators sat- 
isfy the local particle number SSR means that there are less states than other- 
wise would be the case which are classed as non-entangled, and more states will 
be regarded as entangled. It is therefore not surprising that additional tests for 
entanglement will result. If further restrictions are placed on the sub-system 
density operator - such as requiring them to correspond to a fixed number of 
bosons again there will be more states regarded as entangled, and even more 
entanglement tests will apply. A particular example is given in SubSection 15.31 
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where the sub-systems are restricted to one boson states. 

The symmetrisation requirement for systems involving identical particles is 
well established since the work of Dirac. There are two types of justification for 
applying the super- selection rules for systems of identical particles. The first 
approach is based on simple considerations and will be outlined here. The sec- 
ond approach is more sophisticated and involves linking the absence or presence 
of SSR to whether or not there is a suitable reference frame in terms of which 
the quantum state is described [31], [35], [3J], [27], [25], [37], [3S], [3J], [23], [3T], 
[52] , [TB] . This approach will be described in SubSection 12.71 and Appendix [11] 
the key idea being that SSR are a consequence of considering the description 
of a quantum state by an external observer (Charlie) whose phase reference 
frame has an unknown phase difference from that of an observer ((Alice) more 
closely linked to the system being studied. Thus, whilst Alice's description of 
the quantum state may violate the SSR, the description of the same quantum 
state by Charlie will not. In the main part of this paper the density operator 
/5 used to describe the various quantum states will be that of the external ob- 
server (Charlie). Note that if the relationship between the phase references is 
known, then the SSR can be challenged (see SubSection 12.81 and Appendix [II]) . 
Returning to the more simple reasons referred to for invoking the superselcction 
rule to exclude quantum superposition states with differing numbers of identical 
particles (both massive and otherwise), these may be summarised as: 

1. No way is known for creating such states. 

2. No way is known for measuring all the properties of such states, even if 
they existed. 

3. Coherence and interference effects can be understood without invoking 
the existence of such states. 

4. The stability of such states against decoherence processes may not be 
great, so even if they could be created, they could rapidly change to other states. 
However, decoherence time scales that are not too short would be acceptable, 
so this last reason is of lesser importance. 

Invoking the physical existence of states that as far as we know cannot 
be made or measured, and for which there are no known physical effects that 
require their presence seems a rather unnecessary feature to add to the non- 
relativistic quantum physics of many body systems or to quantum optics, and 
considerations based on the general principle of simplicity (Occam's razor) would 
suggest not doing so until a clear physical justification for including them is 
found. Furthermore, experiments can be carried out on each of the mode sub- 
systems considered as a separate system, and essentially the same reasons that 
justify applying the super-selection rule to the overall system also apply to 
the separate mode sub-systems in the context of defining non- entangled states. 
Hence, unless it can be justified to ignore the super-selection rule for the overall 
system it would be inconsistent not to apply it to the sub-system as well. The 
onus is on those who wish to ignore the super-selection rule for the separate 
modes to justify why it is being applied to the overall system. In addition, joint 
measurements on all the sub-systems can be carried out, and the interpretation 
of the measurement probabilities requires the density operators for the sub- 
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system states to be physically based. The general application of super-selection 
rules has however been challenged (see SubSection l2.7[) on the basis that super- 
selection rules are not a fundamental requirement of quantum theory, but are 
restrictions that could be lifted if there is a suitable system that acts as a 
reference for the coherences involved. In Section [5] and related Appendix [T^] 
an analysis of these objections to the super-selection rule is presented, and in 
Appendix [TT] we see that the approach based on phase reference frames does 
indeed justify the application of the SSR both to the general quantum states 
for multi-mode systems of identical particles and to the sub-system states for 
non-entangled states of these systems. 

The other focus of this paper is on spin squeezing. Spin squeezing was first 
defined in the paper by Kitagawa et al [40] for general spin systems. It was sug- 
gested in this paper that correlations between the individual spins was needed 
to produce spin squeezing, though no quantitative proof was presented and the 
more precise concept of entanglement was not mentioned. For the case of two 
mode systems the earliest paper linking spin squeezing to entanglement is that 
of Sorensen et al j?T], which considers a system of identical bosonic atoms, each 
of which can occupy one of two internal states. This paper states that spin 
squeezing requires the quantum state to be entangled, with a proof given in the 
Appendix. A consideration of how such spin squeezing may be generated via 
collisional interactions is also presented. The paper by Sorensen et al is often 
referred to as establishing the link between spin squeezing and entanglement - 
see for example Micheli et al [42] , Toth et al [43] , Hyllus et al [44] . However, the 
paper by Sorensen et al [41j is based on a definition of non-entangled states in 
which the sub-systems are the identical particles, and this is inconsistent with 
the symmetrization principle. The present paper establishes the link between 
spin squeezing and entanglement based on a definition of entanglement consis- 
tent with the system and sub-system density operators representing physical 
states. 

It is also important to consider which components of the spin operator vector 
are squeezed, and this issue is also considered in the present paper. In the 
context of the present second quantisation approach to identical particle systems 
the three spin operator components for two mode systems are expressed in 
terms of the annihilation, creation operators for the two chosen modes. Spin 
squeezing can be defined (see Section [3]) in terms of the variances of these 
spin operators, however the covariance matrix for the three spin operators will 
in general have off-diagonal elements, and spin squeezing is better defined in 
terms of rotated spin operators referred to as principal spin operators for which 
the covariance matrix is diagonal. The principal spin operators are related to 
new mode annihilation, creation operators in the same form as for the original 
spin operators, where the new modes are two orthogonal linear combinations 
of the originally chosen modes. In discussing the relationship between spin 
squeezing and entanglement, the modes which may be entangled are generally 
those associated with the definition of the spin operators. 

The plan of this paper is as follows. In Section [2] the key definitions of 
entangled states are covered, along with detailed discussion on why the sym- 
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metrisation principle and the super-selection rule is invoked. Challenges to the 
necessity of the super-selection rule are outlined, with arguements against such 
challenges dealt with in Appendices [Til H2I an d H31 The next Section [3] sets 
out the definitions of spin squeezing and in the following Section 0] it is shown 
that spin squeezing is a signature of entanglement, both for the principle spin 
operators with entanglement of the two new modes and for the original spin 
operators with entanglement of the original modes. A number of other tests 
for entanglement proposed by other authors are considered in Section O with 
details of these treatments set out in Appendices HH [15l [lj)] Two key mathe- 
matical inequalities are derived in Appendix [5] The final Section [7] summarises 
and discusses the key results. 
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2 Entanglement 



2.1 Physical States 

The standard quantum theory notions of physical systems that can exist in var- 
ious states and have associated quantities on which measurements can be made 
are presumed in this paper. The measuring system made be also treated via 
quantum theory, but there is always some component that behaves classically, 
so that quantum fluctuations in the quantity recorded by the observer are small. 
The term physical state refers to a state that can either be prepared via a pro- 
cess consistent with the laws of quantum physics and on which measurements 
can be then performed and the probabilistic results predicted from this state 
(prediction), or a state whose existence can be inferred from later quantum mea- 
surements (retrodiction) . In quantum theory, physical states are represented by 
density operators for mixed states or state vectors for pure states, which must 
satisfy symmetrisation and other basis requirements in accordance with the laws 
of quantum theory. The quantum state, the system it is associated with and 
the quantities that can be measured are viewed here as entities that are viewed 
as being both ontological and epistimological. The observer is important, but 
there is actually something out there to be studied. In addition to those asso- 
ciated with physical states, other density operators and state vectors may be 
introduced for mathematical convenience. For physical states, the density op- 
erator is determined from either the preparation process or inferred from the 
measurement process, and in general it is a statistical mixture of density opera- 
tors for possible preparation processes. Measurement itself constitutes a possible 
preparation process. Following preparation, further experimental processes may 
change the physical state and dynamical equations give the time evolution of the 
density operator between preparation and measurement, the simplest situation 
being where measurement takes place immediately after preparation. A full dis- 
cussion of the predictive and retrodictive aspects of the density operator is given 
in papers by Pegg et al [15] , [35] . Whilst there are often different mathematical 
forms for the density operator that lead to the same predictive results for subse- 
quent measurements, applying the results of the measurements to retrodictively 
determines the preferred form of the density operator that is consistent with 
the available preparation and measurement operators. An example is given in 

2.2 Entangled and Non-Entangled States 
2.2.1 General Considerations 

Here the commonly applied mathematical approach to defining entangled states 
will be described [8]. The formal definition of what is meant by an entangled 
state starts with the pure states, described via a vector in a Hilbert space. The 
formalism of quantum theory treats pure states for systems made up of two 
or more distinct sub-systems via tensor products of sub-system states \<&a) ® 
\&b) €5 \^c) ■■■■ Such products are called non-entangled or separable states. 
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However, since these product states exist in a Hilbert space, it follows that 

linear combinations of such products of the form |$) = ^2 Cq^7.. ^b 

a/37.. 

I $2?) --could also represent possible pure quantum states for the system. Such 
quantum superpositions which cannot be expressed as a single product of sub- 
system states are known as entangled (or non- separable) states. 

The concept of entanglement can be extended to mixed states, which are 
described via density operators in the Hilbert space. If A, B, ...are the sub- 

^B 

systems with p R , p R , being density operators the sub-systems A, B, .then a 
general non-entangled or separable state is one where the overall density opera- 
tor /5 can be written as the weighted sum of tensor products of these sub-system 
density operators in the form 

R 

with J^r Pr — 1 an d Pr > giving the probability that the specific product 
state /3fj = ®7>r\ ®7>r ® --occurs. Entangled states (or non-separable states) 
are those that cannot be written in this form, so in this approach knowing what 
the term entangled state refers to is based on first knowing what the general 
form is for a non-entangled state. 

The alternative operational approach to defining entangled states focuses on 
whether or not they exhibit certain non-classical features such as Bell Inequality 
violation, or whether they can be used in various quantum information protocols, 
or even whether the states cannot be created via local operations and classical 
communications {LOCC). Such an approach is discussed by Bartlett et al (see 
P5] . Section IIB). As will be seen in SubSection l2.9[ the particular definition of 
entangled states based on their non-creatability.via LOCC essentially coincides 
with the approach used in the present paper. On the other hand, Wiseman et 
al |37], [3H] and Reid et al [TU] discuss the concept of a heirarchy of entangled 
states, with states exhibiting Bell nonlocality being a subset of states for which 
there is EPR steering, which in turn is a subset of the entangled states, the 
latter being defined as states whose density operators cannot be written as 
in Eq. (UJ though without further consideration if additional properties are 
required for the sub-system density operators. The operational approach could 
lead into a quagmire of differing interpretations of entanglement dependng on 
which non-classical feature is highlighted, so the present mathematical approach 
is generally favoured [5]. It is also compatible with later classifying entangled 
states in a heirarchy. 



2.2.2 Local Systems and Operations 

As pointed out by Vedral [7], one reason for calling these states separable is 
associated with the idea of performing operations on the separate sub-systems 
that do not affect the other sub-systems. Such operations on such local systems 
are referred to as local operations and include unitary operations Ua, Ub, that 
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change the states via — > UaPrU^ , p R — > UbPrU r , etc as in a time evolu- 

^~A 

tion, and could include processes by which the states p R , p R , are separately pre- 
pared from suitable initial states. We note that performing local operations on 
a separable state only produces another separable state, not an entangled state. 
Such local operations are obviously faciltated in experiments if the sub-systems 
are essentially non-interacting - such as when they are spatially well- separated, 
though this does not have to be the case. The local systems and operations 
could involve sub-systems whose quantum states and operators are just in dif- 
ferent parts of Hilbert space, such as for cold atoms in different hyperfine states 
even when located in the same spatial region. Note the distinction between 
local and localised. If one observer (Alice) is associated with preparing separate 
sub-system A in a physical state ~p R via local operations with a probability P R , 
a second observer (Bob) could be then advised via a classical communication 
channel to prepare sub-system B in state ~p R via local operations. The overall 
quantum state prepared by both observers via this local operation and classical 
communication protocol (LOCC) would then be the bipartite non-entangled 
state /3 = J2r PrPr ® Pit- Multipartite non-entangled states of the form (TTJ) 
can also be prepared via LOCC protocols involving further observers. As will be 
seen, the separable or non-entangled states are just those that can be prepared 
by LOCC protocols. 



2.2.3 Mean Values and Correlations 

For non-entangled states as in Eq. fllj the mean value for measuring a physical 
quantity Ha ® Hb ® He ® — j where Ha, ^b, He, ■■ are Hermitian operators 
representing physical quantities for the separate sub-systems, is given by 



n A ® n B ® n c ® = p R (&a) A (& b ) B (&c) C r ... (2) 



/- \ K - ~K 

where {HkJ = Tr(£ljc Pr), [K = A,B, ..), is the mean value for measuring 

Hk hi the K sub-system when its density operator is p^. Since the overall mean 
value is not equal to the product of the separate mean values, the measurements 
on the sub-systems are said to be correlated. However, for the general non- 
entangled state as the mean value is just the products of mean values weighted 
by the probability of preparing the particular product state - which involves a 
LOCC protocal, as we have seen - the correlation is classical rather than quan- 
tum [5j . In the case of a single product state where p = p A ® /3 s ® p c ® . .we have 

Ha ® Hb ® He ® = {Ha ^ (^ c ) ■•• wmcn i s j us t the product of 

mean values for the separate sub-systems, and in this case the measurements 
on the sub-systems are said to be uncorrelated. For entangled states however 
the last result for (Ha ®Hb® He ® does not apply, and the correlation is 
strictly quantum. 
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2.2.4 Constraints on Sub-System Density Operators 

A key issue however is whether density operators p and p^, p^, in Eq. ([1} 
always represent possible physical states, even if the operators p and p^, p^, 
etc satisfy all the standard mathematical requirements for density operators - 
Hermitiancy, positiveness, trace equal to unity, trace of density operator squared 
being not greater than unity. In this paper it will be argued that there are further 
requirements not only on the overall density operator, but also on those for the 
individual sub-systems that are imposed by symmetrisation and super-selection 
rules. 

2.3 Joint and Separate Measurements, Reduced Density- 
Operator 

2.3.1 Joint Measurements 

Measurements can be carried out on all the sub-systems and the results ex- 
pressed in terms of the joint probability for various outcomes. If Qa is a physical 
quantity associated with sub-system A, with eigenvalues \f and with Hf the 
projector onto the subspace with eigenvalue \f , is a physical quantity asso- 
ciated with sub-system B, with eigenvalues A^ and with 11^ the projector onto 
the subspace with eigenvalue etc., then the joint probability Pab.X^ 

that measurement of ft a leads to result A^ , measurement of f2s leads to result 
A^ ,etc is given by 

P AB ..(i,j,..) = Tr(n.t^f -P) (3) 
This joint probability depends on the full density operator p representing the 
physical state as well as on the quantities being measured. 

2.3.2 Non-Entangled States - Joint Measurements 

In the case of the general non-entangled state we find that the joint probability 
is 

PAB.XhJ,-)=J2 P K P A^ P Bij)~ (4) 
R 

where 

P$ (i) = Tr(fif p£) Pi(j) = Tr(Tlfp B R ) .. (5) 

are the probabilities for measurement results for Q,a, ^b, ..on the separate sub- 

^b 

systems with density operators p R , p R , etc and the overall joint probability is 
given by the products of the probalbilities P^(i), P§(j), -.for the measurement 
results Aj , A • , ..for physical quantities Qa, £Ib, --if the sub-systems are in the 
states p Rl p Rl etc. These products are then weighted by the probability Pr 

-*nj4 -*^B ^~-C 

that the system is prepared in the particular product state p R ® p R ® p R & ..to 
determine the overall joint probability P A B..(h3i ••)• The overall probability is 
of a classical form. Obviously this joint probability depends on the sub-system 
density operators p fl , pf , etc. 
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2.3.3 Single Measurements and Reduced Density Operator 

The reduced density operator ~p A for sub-system A given by 

p A = Tr B,c,...(p) (6) 

and enables the results for measurements on sub-system A to be determined for 
the situation where the results for all joint measurements involving the other 
sub-systems are discarded. The probability -Pa(*) that measurement of Qa leads 
to result \f irrespective of the results for meaurements on the other sub-systems 
is given by 

P A (i) = £ P AB ..( i ,j,..)=Tr(Ut.p) 

j,k,~ 

= Tr A (Ufp A ) (7) 
using H? = 1, etc. Hence the reduced density operator ~p A plays the role of 

3 

specifying the physical state for mode A considered as a separate sub-system, 
even if the original state /3 is entangled. 

For the general non-entangled state, the reduced density operator for sub- 
system A is given by 

R 

A key feature of a non-entangled state is that the results of a measurement on 
any one of the sub-systems is independent of the states for the other subsystems. 
From Eqs.Q and (JgJ the probability -FU(i) that measurement of fi^ leads to 
result \f is given by 

P A (<) = £ P kPaW (9) 

R 

where the reduced density operator is given by Eq. (|8|) for the non-entangled 
state in Eq. (p}. This result only depends on the reduced density operator ~p A , 
which represents a state for sub-system A and which is a statistical mixture of 
the sub-system states /5^, with a probability Pr that is the same for all sub- 
systems. The result for the measurement probability Pa(*) is just the statistical 
average of the results that would apply if sub-system A were in possible states 
/5^. For all quantum states the final expression for the measurement probability 
Pa(£) only involves a trace of quantities Hf , f) A that apply to sub-system A, 
but for a non-entangled state the reduced density operator ~p A is given by an 
expression ([8]) that does not involve density operators for the other sub-systems. 
Thus for a non-entangled state, the probability PA{i) is independent of the 
states p R , p R , associated with the other sub-systems. Analogous results apply 
for measurements on the other sub-systems. 
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2.4 Bell Inequalities 
2.4.1 Hidden Variable Theory 

A key feature of entangled states is that they are associated with violations of 
Bell inequalities [3 and hence can exhibit this particular non-classical feature. 
The Bell inequalities arise in attempts to restore a classical interpretation of 
quantum thory via hidden variable treatments, where actual values are assigned 
to all measureable quantities - including those which in quantum theory are 
associated with non-commuting Hermitian operators. The actual values are 
determined from the hidden variables, which are associated with a classical 
probability distribution (see for example Section 7.1 of [7]). The Bell inequalities 
involve the mean value E{Ai x Bj) of the product of observables Ai and Bj for 
subsystems A, B respectively, for which there are two possible measured values, 
+1 and —1. In hidden variable theory.the mean values E(Ai x Bj) are given by 



where Ai(X) and Bj(X) are the values are assigned to Ai and Bj when the hidden 
variable is A, and P(X) is the hidden variable probability distribution function. 
The possible values for Ai(X) and Bj(X) are +1 and —1. The form given by 
Clauser et al |4] for Bell's inequality is |5| < 2 where 



The minus sign can actually be attached to any one of the four terms. A proof 
of the Bell inequalities may be found in [7] , [8] . 

2.4.2 Non-Entangled State Result 

It can be shown that the Bell inequalities always occur for non-entangled states 
(see Section 7.3 of the book by Vedral [7]). For Bell's inequalities we consider 
Hermitian operators Aj and Bj for subsystems A, B respectively, for which 
there are two eigenvalues +1 and —1, where the operators are given by the 
components A, = ctj • a a an d Bj — bj ■ as of Pauli spin operators a a and a b 
along directions with unit vectors ai and bj . The corresponding quantum theory 
quantity for the Bell inequality is 




(10) 



S = E{Ax x B x ) + E(Ai x B 2 ) + E{A 2 x B x ) - E{A 2 x B 2 ) 



(11) 



S = E{A\ <g> B x ) + E{A X ® B 2 ) + E{A 2 <g> B{) - E(A 2 <g> B 2 ) 



(12) 



where in quantum theory the mean value is given by E(Ai®Bj) = {Ai ® Bjj = 

Tr(j5Ai <E> Bj). For the general bipartite non-entangled state given by [T] it is 
easy to show that 




R 



(13) 
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where (^-i^ = Tr(Aiffy and (Bj^ = Tr(Bjp^) are the expectation values 
of Ai and Bj for the sub-systems A, B in states and respectively. Now 

A I ^ \ B ' i ^ ~ \ B 



Aij and (-Bjy must lie in the range —1 to +1, so that (Bi ± B2J must 
each lie in the range —2 to +2. Hence 



A 1 ^ ^ \ B 1 _ > A / ^ ^ \ S 



|5| < f I (A x ) ||(B 1+ B 2 ) 1 + 1(^)^11(^-5: 



< 



2 , 

H \ / R 



R x 

< 2 (14) 

since to obtain | (l?i + i? 2 ) \ = 2 requires (-Bi) = = ±1 an d then 

^! - I = I (jj^ - (B-if I - 0, or to obtain | (b x - B 2 )^ | = 2 

requires (b^ = - (s 2 ) = ±1 and then |(Si + B 2 ) B | = | (-§1^ + 
B 2 ) 1=0. 

I R 

2.4.3 Bell Inequality Violation and Entanglement 

It follows that for a general two mode non-entangled state \S\ cannot violate 
the Bell inequality limit of 2. Thus, the violation of Bell inequalities proves 
that the quantum state must be entangled for the sub-systems involved, so Bell 
inequality violations are a test of entanglement. For entangled states such as 
the Bell state \^~) (see [8], Section 2.5) the Bell inequality can be violated for 
the choice where a\, a 2 and b±, 6 2 are orthogonal and a±, a 2 are parallel to 
bi — 6 2 , 61 + 6 2 respectively (see [5], Section 5.1). Furthermore, an entangled 
quantum state and measurements on it cannot be described via a hidden variable 
theory - which caused Einstein to consider that quantum theory violated realism. 
Violation of Bell inequalities is clearly a non-classical feature, since classical 
physics envisages all physical quantities such as Ai and Bj having well-defined 
values, even if they are not known and have to be described via a probability 
distribution P(X) of hidden quantities A which determine their values. 

2.5 Non-local Correlations 

Another feature of entangled states is that they are associated with strong cor- 
relations for observables associated with localised sub-systems that are well- 
separated, a particular example being EPR correlations between non-commuting 
observables. Entangled states can exhibit this particular non-classical feature, 
which again cannot be accounted for via a hidden variable theory. 
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2.5.1 Hidden Variable Theory 

Consider two operators fi a and associated with sub-systems A and B. These 
would be Hermitian if observables are involved, but for generality this is not 
required. In a hidden variable theory these would be associated with functions 
Qc(^) (C = A,B) of the hidden variable A, with the Hermitean adjoints Qq 
being associated with the complex conjugates f2£t(A). Correlation functions 
given by the following mean values 



E(sr A xn B ) = j d\p(\)n* A (\)n B (\) 
E(n* A n A x si* B n B ) = J d\p(x)n A (x)n A (x) n* B (x)n B (x) (15) 

satisfy the following correlation inequality 

\E(n* A x n B ) 1 2 < E(n* A ci A x n* B n B ) (16) 

This result is based on the inequality 

JdXP(X)C(X) > (/dAF(A)VC(A)) 2 (17) 

for real, positive functions C(A),P(A) and where J dXP(X) = 1, and which is 

proved in AppendixlHl In the present case we have C{ X) — Vt* A (X)Q, A (X) VL* B {X)Vt B (X), 
which is real, positive. A violation of the inequality in Eq. (|16p is an indication 
of strong correlation between sub-systems A and B. 

2.5.2 Non-Entangled State Result 

It can be shown that the correlation inequalities are always satisfied for non- 
entangled states. In quantum theory the correlation functions are given by 

E(n A ® h B ) = (h A ®h B ) = Tr(pW A <8> h B ) and E(h A h A <gi h B n B ) = 



il A n A ® Q B Q B J = Tr(pQ A Cl B ® tt B Q B ). For a non-entangled state of sub- 
systems A and B we have 

E{h A ®h B ) = ^p R (nl^ A R (h B ) B 



R 



I R 

E(tt A n A ®ci B n B ) = J2 p r ("a"a) (n^n B ) (is) 



Now 



A v B 



\E{n A ®n B )\<Y,PR\{® A y \ \(&bY\ (is) 

R ' 



since the modulus of a sum is always less than the sum of the moduli. Using 

2 C — 



~ V^c) J ( &c - (ficr ) ) > > with (C = A, B), we obtain the Schwarz 
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inequality - which is true for all states - (riffle y > y^c/ y^ c / = I y^ c 
fit,) I 2 , and hence 



\E(h A ®h B )\ <^P H y(^0^^ W^fi s ^ (20) 

-R 



Next we use the inequality 

E^^>fE^V^) (2i) 

for real, positive functions Cr,Pr and where YIrPr = 1- This inequality, 
which was used in the paper by Hillery et al |29j , is proved in Appendix [9l In 

the present case we have Cr = (q a Qa^ (&b^b^ so * na * 



\E{h\®h B )\ 2 <J2 Pr (&a&a) (o^n B ) = E(h\h A ®h B n B ) (22) 

Thus for a non-entangled state we obtain the correlation inequality 

>t A7,6„\l 2 — I /o . (P, ot \ |2 /ot O . (P, fit 



n\ ® ^s; I = I <«a ® ) I < ^aHa ® ^L^s/ ( 23 ) 

where the general result (q a ®Qb^ = ( ® nas been used. Thus 

non-entangled states have correlation functions that are consistent with hidden 
variable theory. 

2.5.3 Weak Correlation Violation and Entanglement 

Hence if it is found that the correlation inequality is violated | (&> A <£> | 2 = 

I (Qa ® &g \ | 2 > (^l A flA ® ^s^s) then the state must be entangled, so the 
correlation inequality violation is also a test for entanglement. Again entangled 
states have features that cannot be explained via hidden variable theory, so 
entangled states are clearly non-classical. .. 

2.6 Identical Particles - Symmetrisation Principle 
2.6.1 Symmetrization and Second Quantization 

Whether entangled or not the physical states for systems of identical particles 
must conform to the symmetrisation principle, whereby the overall density oper- 
ator has to be invariant under permutation operators. Problems arise regarding 
how to define non-entangled states for systems of identical particles. The basic 
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issue is how first to distinguish what are meaningful sub-systems for identical 
particle systems. Some authors ([41], [49].. ) consider states of the form 

R 

to be non-entangled states, where is a density operator for particle i. However 
such a state would not in general be physical, since the symmetrisation principle 
would be violated unless the ~p R were related. For example, consider the state 
for two identical bosonic atoms given by 

p = P a( : a 1 ®f + PerjO 1 (25) 

and apply the permutation P = P(l 2). The invariance of /3 in general 

requires a — £ and 9 = 17, giving p — P a a ® a + P$6 <£> . This is a 
statistical mixture of two states, one with both atoms in state a, the other with 
atoms in state 9. Of course if the atoms were all different (atom 1 a Rb 87 atom, 
atom 2 a Na 23 atom, ..) then the expression (|251) would be a valid non-entangled 
state, but there the atomic sub-systems are distinguishable and symmetrisation 
is not required. What is distinguishable for systems of identical bosons is not the 
individual particles themselves - which do not carry labels, boson 1, boson 2, etc. 
- but the single particle states or modes that the bosons may occupy. For bosonic 
atoms with several hyperfine components, each component will have its own set 
of modes. The same would apply to fcrmionic atoms. For photons the modes 
may be specified via wave vectors and polarisations. Although the quantum 
pure states can be specified via symmetrized products of single particle states 
occupied by specific particles using a first quantization approach, it is more 
convenient to use second quantization. Here, a basis set for the quantum states 
of such sub-systems are the Fock states \n a ) (n a — 0, 1, 2, ..) etc, which specify 
the number of identical particles occupying the mode A, etc., so in this approach 
the mode is the sub-system and the Fock states give different physical states for 
this sub-system. Symmetrization is built into the definition of the Fock states. If 
the atoms were fermions rather than bosons the Pauli exclusion principle would 
of course restrict n a = 0,1 only. Thus in this second quantization approach 
situations with differing numbers of identical particles are different states, not 
different systems. The overall system will be associated with physical states 
with density operators and state vectors in Fock space, which includes states 
with total numbers of identical particles ranging from zero in the vacuum state 
right up to infinity. 

2.6.2 Sub-Systems and Modes 

The point of view in which the possible sub-systems A, B, etc are modes (or 
sets of modes) rather than particles has been adopted by several authors 
[12], [13]), [7], [21], [22] and will be the approach used here. What are or 
are not entangled are modes not particles. Overall, the system is a collection 
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of modes, not particles. Particles are associated with mode occupancies, and 
therefore related to specifying the quantum states of the system, rather than 
the system itself. Note that in this approach states where there is only a single 
atom may still be entangled states - for example with two spatial modes A, B 
the states which are a quantum superposition of the atom in each of these 
modes, such as the Bell state (|1 Q ) |0f,) + |0 a ) |lb))/V2 are entangled states. For 
entangled states associated with the EPR paradox or for quantum teleportation, 
the mode functions may be localised in well-separated spatial regions - spooky 
action at a distance - but spatially overlapping mode functions apply in other 
situations. Furthermore, as well as being distinguishable the modes can act 
as separate systems, with other modes being ignored. For interacting bosonic 
atoms this is much harder to accomplish experimentally than for the case of 
photons, where the relatively slow processes in which photons are destroyed in 
one EM field mode and created in another may require the presence of atoms as 
intermediaries. Two bosonic atoms in one mode may collide rapidly disappear 
into other modes. However, atomic boson interactions can be made very small 
via Feshbach resonance methods. Near absolute zero the basic physics of a BEC 
in a single trap potential is describable via a one mode theory. Hence with A, 
B, .. signifying distinct modes, the general non-entangled state is given in Eq. 
(fT) though the present paper mainly involves only two modes. 

Note however that a different concept of entanglement - particle entangle- 
ment - has also been applied to identical particle systems [T5], [U]. This is 
not the same as mode entanglement so tests and measures for particle entan- 
glement will differ from those for mode entanglement. For completeness a brief 
description highlighting the difference between mode and particle entanglement 
is presented in Appendix 1101 A further discussion about the distinction is given 
in [IB]. 

2.6.3 Multi-Mode Sub-Systems 

As well as the simple case where the sub-systems are all individual modes, the 
concept of entanglement may be extended to situations where the sub-systems 
are sets of modes, rather than individual modes, In this case entanglement or 
non-entanglement will be of these distinct sets of modes. Such a case in con- 
sidered in Subsection 15 .31 where pairs of modes associated with distinct lattice 
sites are considered as the sub-systems. Another example is treated in He et al 
[50] . which involves a double well potential with each well associated with two 
bosonic modes, these pairs of modes being the two sub-systems. Entanglement 
criteria for the mode pairs based on local spin operators associated with each 
potential well are considered (sse SubSection l5.6p . A further example is treated 
by Heaney et al [51] , again involving four modes associated with a double well 
potential. As in the previous example, each mode pair is associated with the 
same well in the potential, but here a Bell entanglement test was obtained for 
pairs of modes in the different wells. The concept of entanglement of sets of 
modes is a straightforward extension of the basic concept of entanglement of 
individual modes. 
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2.7 Super-Selection Rule 

2.7.1 Global Particle Number SSR 

The question of what physical states - entangled or not - are possible in the 
non-relativistic quantum physics of a system of identical bosonic particles - such 
as bosonic atoms or photons - has been the subject of much discussion. Whether 
entangled or not it is generally accepted that there is a super-selection rule that 
prohibits quantum superposition states of the form 

OO OO OO OO 

|4>) = £ G N \N) p= £ \Gn\ 2 \N) (N\+ £ £ (1-5 N , M )C N C* M \N) (M\ 

N=0 N=0 N=QM=0 

(26) 

being physical states when they involve Fock states \N) with differing total 
numbers N of particles. The density operator for such a state would involve 
coherences between states with differing N. Although such superpositions - 
such as the Glauber coherent state \a), where Cn — exp(— |a| 2 /2) a N /#i - do 
have a useful mathematical role, they do not represent actual physical states 
according to the super-selection rule. The papers by Sanders et al [3 6) and 
Cable et al [52] are examples of applying the SSR for optical fields, but also 
using the mathematical features of coherent states to treat phenomena such as 
interference between independent lasers. The super-selection rule indicates that 
the most general physical state for a system of identical bosonic particles can 
only be of the form 

p = £ (|*tf><**|) 

N=0 * 

\$n) = £Cf \Ni) (27) 

i 

where |$jv) is a quantum superposition of states \Ni) each of which involves 
exactly N particles, and where different states with the same N are designated 
as \Ni). This state p is a statistical mixture of states, each of which contains 
a specific number of particles. Such a SSR is referred to as a global SSR, as it 
applies to the system as a whole. Mathematically, the global particle number 
SSR can be expressed as 

[N,p]=0 (28) 
where N is the total number operator. 

2.7.2 Examples of Global Particle Number SSR Compliant States 

Examples of a state vector | <J?jv) for an entangled pure state jT^] and a density 
operator /5 for a non-entangled mixed [53] state for a two mode bosonic system, 
both of which are possible physical states are 

N 

\* N ) = J2C(N,k)\k) A ®\N-k) B (29) 

k=0 
N 

P = J2P(k)\k) A (k\ A ®\N-k) B (N-k\ B (30) 

k=0 
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The entangled pure state is a superposition of product states with k bosons 
in mode A and the remaining N — k bosons in mode B. Every term in the 
superposition is associated with the same total boson number N. The non- 
entangled mixed state is a statistical mixture of product states also with k 
bosons in mode A and the remaining N — k bosons in mode B. Every term in 
the statistical mixture is associated with the same total boson number N. For 
the case of a two mode fermionic system the Pauli exclusion principle restricts 
the number of possible fcrmions to two, with at most one fermion in each mode. 
Expressions for a state with exactly N — 2 fermions are 

1*2) = |1> A ®|1> B (31) 
P = |1)a(1U ®|1>b<1| b (32) 

Neither state is entangled and both are the same pure state since p = |<&2) (3>2|- 
Although the super-selection rules and symmetrisation principle also applies 
to fermions, as indicated in the Introduction this paper is focused on bosonic 
systems, and it will be assumed that the modes are bosonic unless indicated 
otherwise. 

The Bell states for N — 2 bosons provide important examples of four mode 
pure quantum states that are compliant with the global particle number SSR. 
The modes are designated A+, A—,B+,B— and the Fock states are in general 
\nA+,nA-,nB+,n>B-,)- The Bell states may be written 

^singlet) = ^ (| 1 , 0, 0, 1) - |0, 1, 1, 0}) = -L (| A+) ® | B-) - \ A-) ® | B) ) 
|*tripiet,+l) - |1,0,1,0) = |A+)®|B+) 
|*tripJet,0) - -^(|l,0,0,l) + |0,l,l,0)) = -^(|A+)®|B-) + |A-)®|B)) 

|1Wet,-i) - |0,1,0,1) = |A-)®|B-) (33) 

where the second forms may be more familiar. Of these states \^ singlet) and 
\^tripiet,o} are entangled, whilst \V triplet, +1} and triplet -\) are separable. 



2.7.3 Super-Selection Rules and Conservation Laws 

It is important to realise that such super-selection rules are additional con- 
straints to those imposed by conservation laws. For example, the conservation 
law on total particle number only leads to the requirement on the superposition 
state |$) that the |Cat| 2 are time independent, it does not require only one CV 
being non-zero. Super-selection rules are broad in their scope, forbidding quan- 
tum superpositions of states of systems with differing charge, differing baryon 
number and differing statistics. Thus a combined system of a hydrogen atom 
and a helium ion docs not exist in quantum states that are linear combina- 
tions of hydrogen atom states and helium ion states - the super-selection rules 
on both charge and baryon number preclude such states. The basis physical 
states for such a combined system would involve symmetrised tensor products 
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of hydrogen atom and helium ion states, not linear combinations - symmetri- 
sation being required because the system contains two identical electrons. On 
the other hand, super-selection rules do not prohibit quantum superpositions of 
states of systems with differing energy, angular or linear momenta - other phys- 
ical quantities that may also be conserved. Thus in a hydrogen atom quantum 
superpositions of states with differing energy and angular momentum quantum 
numbers are allowed physical states. 

2.7.4 SSR Justification and No Suitable Phase Reference 

There are two types of justification for applying the super-selection rules for sys- 
tems of identical particles. The first approach is based on simple considerations 
and will be outlined below in this subsection. The second approach [33], [35], 
[36] , [37] , [38] , [39] , [26] , [31] , [32] , [16] is more sophisticated and involves linking 
the absence or presence of SSR to whether or not there is a suitable reference 
frame in terms of which the quantum state is described, and is outlined in the 
next subsection and Appendix [11] The key idea is that SSR are a consequence of 
considering the description of a quantum state by an external observer (Charlie) 
whose phase reference frame has an unknown phase difference from that of an 
observer ((Alice) more closely linked to the system being studied. Thus, whilst 
Alice's description of the quantum state may violate the SSR, the description 
of the same quantum state by Charlie will not. In the main part of this paper 
the density operator p used to describe the various quantum states will be that 
of the external observer (Charlie). 

2.7.5 SSR Justication and Physics Considerations 

A number of straightforward reasons have been given in the Introduction for 
why it is appropriate to apply the superselection rule to exclude quantum su- 
perposition states of the form (|2"6"]l as physical states for systems of identical 
particles, and these will now be considered in more detail. 

Firstly, no way is known for creating such states. The Hamiltonian for such a 
system commutes with the total boson number operator, resulting in the |Cjv| 2 
remaining constant, so the quantum superposition state would need to have 
existed initially. In the simplest case of non-interacting bosonic atoms, the Fock 
states are also energy eigenstates, such Fock states involve total energies that 
differ by energies of order the rest mass energy mc 2 , so a coherent superposition 
of states with such widely differing energies would at least seem unlikely in 
a non-relativistic theory, though for massless photons this would not be an 
issue as the energy differences are of order the photon energy Huj. The more 
important question is: Is there a non-relativistic quantum process could lead to 
the creation of such a state? Processes such as the dissociation of M diatomic 
molecules into up to 2M bosonic atoms under Hamiltonian evolution involve 
entangled atom-molecule states of the form 

M 

|$) = £ C m \M - m) mol ® \2m) atom (34) 

rn—O 
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but the reduced density operator for the bosonic atoms is 

Ktoms = E \C m \ 2 (|2m) (2m\) aUm (35) 

m=0 

which is a statistical mixture of states with differing atom numbers with no 
coherence terms between such states. Such statistical mixtures are valid physical 
states, corresponding to a lack of a priori knowledge of how many atoms have 
been produced. To obtain a quantum superposition state for the atoms alone, 
the atom-molecule state vector would need to evolve at some time into the form 

M M 

|$) = £ B m \M - m) mol ® £ A 2n \2n) atom (36) 

m— n—0 

where the separate atomic system is in the required quantum superposition 
state. However if such a state existed there would be terms with at least one 
non-zero coefficient B m A 2n involving product states \M — m) mol ® |2n) atom 
with n 7^ to if the state |$) is not just in the entangled form (fMl) . However, 
the presence of such a term would mean that the conservation law involving the 
number of molecules plus two times the number of atoms was violated. This is 
impossible, so such an evolution is not allowed. 

Secondly, no way is known for measuring all the properties of such states, 
even if they existed. If a state such as (|2"oT) did exist then the amplitudes Cm 
would oscillate with frequencies that differ by relativistic frequencies of order 
mc 2 /h, even if boson-boson interactions were included To distinguish the phases 
of the Cn in order to verify the existence of the state, measurement operators 
would need to include terms that also oscillate at relativistic frequencies, and 
no such measurement operators are known. 

Thirdly, there is no need to invoke the existence of such states in order to un- 
derstand coherence and interference effects.. It is sometimes thought that states 
involving quantum superpositions of number states are needed for discussing 
coherence and interference properties of BECs, and some papers describe the 
state via the Glauber coherent states. However, as Leggett [51] has pointed out 
(see also Bach et al [55], Dalton and Ghanbari [56 J, a highly occupied number 
state for a single mode with N bosons has coherence properties of high order n, 
as long as n <C N. The introduction of a Glauber coherent state is not required 
to account for coherence effects. Even the well-known presence of spatial inter- 
ference patterns produced when two independent BECs are overlapped can be 
accounted for via treating the BECs as Fock states. The interference pattern is 
built up as a result of successive boson position measurements [57] , [36] , [52] . 

Fourthly, the stability of such states against decoherence processes may not 
be great, so even if they could be created, they could rapidly change to other 
states. However, decoherence time scales that are not too short would be ac- 
ceptable. Although BECs are created in high vacuum experiments and are well 
isolated from the external environment in magnetic or dipole traps, they are not 
entirely free from decoherence effects because the bosons do interact with each 
other. Even in a single mode case boson-boson collisions can cause dephasing 
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effects. These could be shown via the decay of the coherence (a). However, 
it may turn out that the lifetime of a coherent state in a single mode BEC is 
quite long - in the case of photons the lifetime could be as long as the inverse 
Townes-Schawlow line width, perhaps of order 10 3 s (see below). If a coherent 
superposition state could be created with a non-zero coherence, this may last 
long enough to carry out further experiments, so this fourth reason for dis- 
carding coherent superposition states is relatively unimportant though further 
studies of their lifetimes would be of some theoretical interest. . 

2.7.6 SSRJustiflcation and Galilean Frames ? 

Finally, in addition to the previous reasons there is an arguement based on the 
requirement that the dynamical equations for such non-relativistic quantum sys- 
tems should be invariant under a Galilean transformation which has been pro- 
posed [55] as a proof of the super-selection rule for atom number. This approach 
is linked to the reference frame based justification of SSR (see Appendix [TTj) . 
However, whilst the paper shows that under a Galilean transformation - corre- 
sponding to describing the system from the point of view of an observer moving 
with a constant velocity v with respect to the original observer, and where the 
two observers have identical clocks - the terms in a superposition state with 
different numbers N of massive bosons would oscillate like expi ^Nmv 2 i) /H, 
such oscillations have non-relativistic frequency differences, and are only to be 
expected if the same quantum state is described by a moving observer. This 
feature alone does not seem to require the super-selection rule, since here the 
moving observer's reference frame has a well-defined velocity with respect to 
that attached to the system so that no twirling operation resulting in the elim- 
ination of number state coherences is involved (see Appendix lllj). and so will 
not be considered further in this paper. 

On the other hand, an approach of this kind involving rotation symmetry 
would seem to rule out such states as quantum superpositions of a boson (spin 
0) and a fermion.(spin 1/2). Let such a state be prepared in the form (\F) + 
\B))/y/2. Consider an observer whose cartesian reference frame is X, Y, Z. This 
is a classical system that can be rotated in space. If the observer rotates with 
his frame through 27r about any axis they are then back in the same position, 
but the observer now sees the state as (— \F) + \B))/ y2. This state is apparently 
orthogonal to the one observed before the rotation, and this is paradoxical since 
the observer would be in the same position. Thus there is a super-selection rule 
excluding states such as (\F) + \ B))/\/2. A similar argument based on the time 
reversal anti-unitary operator was given by Wick et al [24 . 

2.7.7 SSR and Photons 

Though this paper is focused on massive bosonic atoms similar considerations 
also apply to the optical quantum EM field, which involve massless bosons - 
photons. Even in the case of photons, Molmer [55] has argued that the physical 
state for a single mode optical laser field operating well above threshold is not a 
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Glauber coherent state, and the density operator would be a statistical mixture 
of the form (JUJ), with \$ N ) = \N) and Rj> N = cxp{-N)N N /N\. Here the 
density operator is a statistical mixture of photon number states with Poisson 
distribution, or equivalently a statistical mixture of coherent states \a) with 
a = ViV exp(£<^>) and all phases (f> having equal probability. The same general 
reasons for applying the super-selection rule to systems of identical massive 
bosons also apply here, though the details differ. For the free quantum EM 
field there is a conservation law for the photon number in each mode, so in 
this case again |Cjv| 2 would be time independent. Here the Cm would oscillate 
with frequencies that only differ by non-relativistic frequencies of order tuv, 
so the arguement against coherent states based on this feature do not apply. 
However, in the case of the single mode optical laser, the field is generated via 
interactions with incoherently pumped atoms, there is no well defined optical 
phase that can be imposed on the process and the quantum theory for such laser 
processes predicts a quantum state that is a statistical mixture of photon number 
states. In the case of the optical laser field coherent states are not physical 
because there are no optical reference fields with a well-defined phase that could 
be used to determine the phases associated with the expansion coefficients. 
Optical interference and coherence effects can also be explained without invoking 
Glauber coherent states, as [5S] and others such as [55] have shown. However, if 
coherent states could be created they might be relatively stable. In the optical 
laser field case, phase loss via diffusion is related to the laser linewidth, and 
this can be reduced to the Townes-Schawlow limit that varies inversely as the 
mean photon number - which is large. The Townes-Schawlow linewidth can 
be as small as 10~ 3 hz, corresponding to a phase diffusion time of 10 3 s. An 
alternative approach is presented by Wiseman et al [60], [61] . in which the 
optical laser is treated via a master equation, but where monitoring of the laser 
environment (difficult!) is required to determine whether certain pure state 
ensembles - such as those involving coherent states - are physically realisable. 
The conclusion reached is that for finite self energy the coherent state ensemble 
is not physically realisable, the closest ensemble being that involving squeezed 
states, though for zero self energy coherent state ensembles are obtained. 

There is however, an approach (see next sub-section and SubSection II 1 .41 in 
Appendix I lip involving the consideration of phase reference frames, in terms of 
which the quantum state of a single mode laser may be described as a Glauber 
coherent state by an observer (Alice) with one reference frame, but would be 
described as a statistical mixture of photon number states by another observer 
(Charlie) with a different reference frame, whose phase reference is completely 
unrelated to the previous one. 

2.8 Reference Frames and Violations of Superselection Rules 

Challenges to the requirement for physical states to be consistent with super- 
selection rules have occured since the 1960's when Aharonov and Susskind [34] 
suggested that coherent superpositions of different charge eigenstates could be 
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created. It is argued that super-selection rules are not a fundamental require- 
ment of quantum theory, but the restrictions involved could be lifted if there is 
a suitable system that acts as a reference for the coherences involved - [31] , [3S] , 
[55] . [37], [3S], [3H], EH], [31], [33, [H] provide discussions regarding reference 
systems and SSR. 

2.8.1 Linking SSR and Reference Frames 

The discussion of the super-selection rule issue in terms of reference systems is 
quite complex and too lengthy to be covered in the body of this paper. However, 
in view of the wide use of the reference frame approach a full outline is presented 
in Appendix [TT] The key idea is that there are two observers - Alice and Charlie 
- who are describing the same quantum state in terms of their own reference 
systems. The reference systems are macroscopic systems in states where the 
behaviour is essentially classical, such as large magnets that can be used to 
define cartesian axes or BEC in Glauber coherent states that are introduced to 
define a phase reference. The relationship between the two reference systems 
is represented by a group of unitary transformation operators listed as T(g), 
where the particular transformation (translation or rotation of cartesian axes, 
phase change of phase references, ..) that changes Alice's reference system 
into Charlie's is denoted by g. Alice is the internal observer, closely linked 
to the system under study and describes the quantum state via her density 
operator, whereas Charlie is the external observer whose specification of the 
same quantum state via his density operator is of most interest. There are 
two cases of importance, Situation A - where the relationship between Alice's 
and Charlie's reference frame is is known and specified by a single parameter g, 
and Situation B - where on the other hand the relationship between frames is 
completely unknown, all possible transformations g must be given equal weight. 
Situation A is not associated with SSR, whereas Situation B leads to SSR. The 
relationship between Alice's and Charlie's density operators is given in terms 
of the transformation operators (see Eq. (|199|) for Situation A and Eq. (|200p 
for Situation B). In Situation B there is often a qualitative change between 
Alice's and Charlie's description of the same quantum state, with pure states 
as described by Alice becoming mixed states when described by Charlie. It is 
Situation B with the U(l) transformation group - for which number operators 
are the generators - that is of interest for the single or multi-mode systems 
involving identical bosons on which the present paper focuses. An example of the 
qualitative change of behaviour for the single mode case is that if it is assumed 
that Alice could prepare the system in a Glauber coherent pure state - which 
involves SSR breaking coherences between differing number states - then Charlie 
would describe the same state as a Poisson statistical mixture of number states - 
which is consistent with the operation of the SSR. Thus the SSR applies in terms 
of external observer Charlie's description of the state. This is how the dispute 
on whether the state for single mode laser is a coherent state or a statistical 
mixture is resolved - the two descriptions apply to different observers - Alice 
and Charlie. On the other hand there are quantum states such as Fock states 
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and Bell states which are described the same way by both Alice and Charlie, 
even in Situation B. The general justification of the SSR for Charlie's density 
operator description of the quantum state in Situation B is derived in terms 
of the irreducible representations of the transformation group, there being no 
coherences between states associated with differing irreducible representations 
(see Eq. (I224jl ). For the particular case of the U(l) transformation group the 
irreducible representations are associated with the total boson number for the 
system or sub-system, hence the SSR that prohibits coherences between states 
where this number differs. Finally, it is seen that if Alice describes a general non- 
entangled state of sub-systems - which being separable have their own reference 
frames - then Charlie will also describe the state as a non-entangled state and 
with the same probability for each product state (see Eqs. (j229[) and (|230p ). 
For systems involving identical bosons Charlie's description of the sub-system 
density operators will only involve density operators that conform to the SSR. 
This is in accord with the key idea of the present paper. 

2.8.2 Coherent Superposition of Atom and Molecule ? 

Based around the reference frame approach Dowling et al [62] and Terra Cunha 
et al [13] propose processes using a BEC as a reference system that would create 
a coherent superposition of an atom and a molecule, or a boson and a fermion 
|62j . Dunningham et al |63j consider a scheme for observing a superposition of a 
one boson state and the vacuum state. Obviously if super-selection rules can be 
overcome in these instances, it might be possible to produce coherent superpo- 
sitions of Fock states with differing particle numbers such as Glauber coherent 
states, though states with N ~ 10 8 would presumably be difficult to produce. 
However, detailed considerations of such papers indicate that the states actually 
produced in terms of Charlie's description are statistical mixtures consistent 
with the super-selection rules rather than coherent superpositions, which are 
only present in Alice's description of the state. Also, although coherence and 
interference effects are demonstrated, these can also be accounted for without 
invoking the presence of coherent superpositions that violate the super-selection 
rule. As the paper by Dowling et al [55] entitled "Observing a coherent superpo- 
sition of an atom and a molecule." is a good example of where the super-selection 
rules are challenged, the key points are described in Appendix [12J Essentially 
the process involves one atom A interacting with a BEC of different atoms B 
leading to the creation of one molecule AB, with the BEC being depleted by 
one B atom. There are three stages in the process, the first being with the 
interaction that turns separate atoms A and B into the molecule AB turned 
on at Feshbach resonance for a time t related to the interaction strength and 
the mean number of bosons in the BEC reference system, the second being free 
evolution at large Feshbach detuning A for a time r leading to a phase factor 
4> = A r, the third being again with the interaction turned on at Feshbach res- 
onance for a further time t. However, it is pointed out in Appendix fTS] that 
Charlie's description of the state produced for the atom plus molecule system 
is merely a statistical mixture of a state with one atom and no molecules and 
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a state with no atom and one molecule, the mixture coefficients depending on 
the phase <f> imparted during the process. However a coherent superposition is 
seen in Alice's description of the final state, though this is not surprising since a 
SSR violating initial state was assumed. The feature that in Charlie's descrip- 
tion of the final state no coherent superposition of an atom and a molecule is 
produced in the process is not really surprising, because of the averaging over 
phase differences in going from Alice's reference frame to Charlie's. It is the 
dependence on the phase <fi imparted during the process that demonstates co- 
herence (Ramsey intcrfcromctry) effects, but it is shown in Appendix Q2] that 
exactly the same results can be obtained via a treatment in which states which 
are coherent superpositions of an atom and a molecules are never present, the 
initial BEC state being chosen as a Fock state. In terms of the description by 
an external observer (Charlie) the claim of violating the super-selection rule has 
not been demonstrated via this particular process. 

2.8.3 Detection of SSR Violating States 

Whether such super-selection rule violating states can be detected has also not 
been justified. For example, consider the state given by a superposition of a 
one boson state and the vacuum state (as discussed in [53] )■ We consider an 
interferometric process in which one mode A for a two mode BEC interferometer 
is initially in the state a |0)+/3 |1), and the other mode B is initially in the state 
|0) - thus — (a \0}+(3 \1))a®\0) b in the usual occupancy number notation, 

where \a\ 2 + \/3\ 2 = 1. The modes are first coupled by a beam splitter, then a 
free evolution stage occurs for time r associated with a phase difference <fi = At 
(where A = lob — is the mode frquency difference), the modes are then 
coupled again by the beam splitter and the probability of an atom being found 
in modes A, B finally being measured. The probabilities of finding one atom in 
modes A, B respectively are found to only depend on |/3| 2 and (j>. Details are 
given in Appendix [12] There is no dependence on the relative phase between 
a and /3, as would be required if the superposition state a |0) + \1) is to be 
specified. Exactly the same detection probabilities are obtained if the initial 
state is the mixed state p(i) = \a\ 2 (\0) A (0\ A ® |0) B (0| B ) + |/3| 2 (|1) A (1\ A O 
|0) B (0| B ), in which the vacuum state for mode A occurs with a probability 
| a | 2 and the one boson state for mode A occurs with a probability |/3| 2 . In this 
example the proposed coherent superposition associated with the super-selection 
rule violating state would not be detected in this interferometric process, nor in 
the more elaborate scheme discussed in [63] . 

2.9 Super-Selection Rule - Separate Sub-Systems 
2.9.1 Local Particle Number SSR 

We now consider the role of the super-selection rule for the case of non-entangled 
states. The global super-selection rule on total particle number has restricted 
the physical quantum state for a system of identical bosons to be of the form 
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PTp. Such states may or may not be entangled states of the modes involved. 
The question is - do similar restrictions involving the sub-system particle number 
apply to the modes, considered as separate sub-systems in the definition of non- 
entangled states ? The viewpoint in this paper is that this is so. Note that 
applying the SSR on the separate sub-system density operators p R , p R , .. is 
only in the context of non-entangled states. Such a SSR is referred to as a local 
SSR, as it applies to each of the separate sub-systems. Mathematically, the 
local particle number SSR can be expressed as 

Wx,Pr]=0 (37) 

where Nx is the number operator for sub-system X = A, B, ...The SSR re- 
striction is based on the proposition that the density operators p R , p R , .. for 
the separate sub-systems A, B, ..should themselves represent possible physical 
states for each of the sub-systems, considered as a separate system and thus be 
required to satisfy the super-selection rule that forbids quantum superpositions 
of Fock states with differing boson numbers. It is contended that expressions 
for the non-entangled quantum state p in which ~p R , p^, p^.. were not physical 
states for the sub-systems would only be of mathematical interest. 

Applying the local particle number SSR to the sub-system density operators 
for non-entangled states is discussed in papers by Bartlett et al [25], [26] as one 
of several operational approaches for defining entangled states. However, other 
authors such as [27] , [25] state on the contrary that states when the sub-system 
density operators do not conform to the local particle number super-selection 
rule are still separable, others such as [25], [3D] do so by implication, so in this 
paper we are advocating a revision to the widely held notion of entanglement in 
identical particle systems, .the consequence being that the set of entangled states 
is now much larger. This is a key idea in this paper - not only should super- 
selection rules on particle numbers be applied to the the overall physical state, 
entangled or not, but it also should be applied to the density operators that 
describe states of the modal sub-systems involved in the general definition of 
non-entangled states. The reasons for adopting this viewpoint are set out below. 
Apart from the papers by Bartlett et al [25] , [25] we are not aware that this def- 
inition of non-entangled states has been invoked previously, indeed the opposite 
approach has been proposed [27], [25J. However, the idea of considering sub- 
system states which satisfy the local particle number SSR has been presented in 
several papers - [27] , [28] , [25] , [26] , [31] , [32] , [33] , mainly in the context of pure 
states for bosonic systems, though in these papers the focus is on issues other 
than the definition of entanglement - such as quantum communication protocols 
[27], multicopy distillation [25], mechanical work and accessible entanglement 
[51] . [52] and Bell inequality violation [55J. However, there are a number of pa- 
pers that do not apply the SSR to the sub-system density operators, and those 
that do have not studied the consequences for various entanglement tests - as 
is done in the present paper. 
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2.9.2 Local SSR Justification and Independent Local Phase Refer- 
ences 

The more elaborate justification in terms of reference frames for this SSR re- 
quirement on non-entangled states is presented in SubSection lll.9l of Appendix 
[TT1 Essentially the idea is that in the context of separable states, each sub- 
system has its own independent phase reference frames, and those of Charlie 
having an unknown phase in relation to those of Alice. This leads to the local 
particle number SSR. 

2.9.3 Local SSR Justication and Physics Considerations 

The more simple reasons for this assertion are analogous to those for the overall 
multi-mode system and may be summarised as: absence of both a preparation 
process and a measurement process for such states, the lack of need of such states 
to describe single mode interference and coherence effects. Such superposition 
states may also be unstable, though again this is not a fatal problem. 

Firstly, sub-system states incompatible with the SSR cannot be prepared. 
Consider for example a typical preparation process. For the situation of two 
modes A, B physically allowed pure states |$jv) could be prepared which in 
general are entangled states of the form 

JV 

\*N)=Y.A%\k) A \N-k) B (38) 

fe=0 

so that the general mixed physical state for the two mode system is 

oo JV N 

P= E E-FW EE A k (A^y \k) A (l\ A ® |JV - k) B (N - l\ B (39) 

JV=0 * k=0 1=0 

Hence the reduced density operator - which specifies the state for mode A if 
measurements on this mode were carried out and measurements on other modes 
discarded - will be given by 

oo N 

pa = E E p*n E AH (A»y \k) A (k\ A (40) 

N=0 <S> fc=0 

which is a statistical mixture of Fock states \k) A . Thus the quantum state 
for mode A considered separately contains no superposition of states \k) A with 
differing numbers of bosons occupying mode A. As in the example considered in 
the previous section, the evolution of \<&n) into a tensor product of superposition 
states for modes A and B of the form 

JV JV 

I^Jv) = E C? \k) A ® £ D% \N - k) B (41) 

k=0 k=0 

is not possible. The preparation of the state for mode A must have involved first 
preparing a physical state for the full multi-mode system - for which the two 
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mode state in Eq. (I3T)1) is a specific example - from which the state associated 
with a particular mode is then determined as given by the reduced density 
operator. As illustrated by the example just given, the super-selection rule on 
the total number of identical bosons for the overall system produces a reduced 
density operator for the sub-system in which the super-selection rule for boson 
number also applies - that is the state for the sub-system does not involve 
quantum superpositions of mode Fock states with differing boson numbers, it 
only can involve statistical mixtures of such states. 

Secondly, measurement processes may be applied to each separate mode 
and again the lack of measurement systems with well defined rclativistic phases 
would preclude measurements that determine the rapidly varying phase dif- 
ferences between the expansion coefficients. in single mode state vectors of the 

oo 

form \§>a) — C n \ n ) a- Invoking the existence of states whose key properties 

n=0 

cannot be measured is somewhat dubious. 

Thirdly, experimental setups involving single mode BECs and optical sys- 
tems can be created and yet there is no need to invoke coherent superpositions of 
number states to explain coherence and interferometric effects. Thus essentially 
the same reasons that justify applying the super-selection rule to the overall 
many boson system also apply to the separate mode sub-systems. 

2.9.4 Local SSR Justification and Joint Measurements 

A consideration of joint measurements on all the sub-systems leads to other 
fundamental reasons why the individual density operators p^, p^, .in the spe- 
cific situation of the general mixed non-entangled state given in Eq. (JTJ) must 
represent physical states for the sub-systems. This state is a statistical mixture 
of product states p R ® p^ ® p^ ® each product state being an overall state of 
the system that could have been prepared. If sub-sysystem A is prepared by one 
experimenter in state with probability Pr, classical communications to other 
local experimenters to prepare the other sub-systems in states Pr, Pr, etc with 
the same probability will result in the preparation of the overall mixed state. If 
such an overall product state is a physical state, then so must be the states of 
the uncorrelated sub-systems involved. Furthermore, measurements on all the 
sub-systems can be carried out, not just those on one particular sub-system A 
- where the results for the sub-system probabilities Pa(*) are determined from 
the reduced density operator p A - see Eq. ©. We have seen in Eq Q that 
the joint probability Pa_b..(«, j, ■•) for measurements on all the sub-systems is 
determined from the product of the individual sub-system probabilities P A (i), 
Pg(j), ..associated with sub-system density operators /3r, p^, the overall 
product being weighted by the probability Pr that a particular product state 
is prepared. The reduced density operators for all the sub-systems do not de- 
termine this joint probability - what is required are the full set of sub-system 
density operators p R , p R , ...along with the overall probability Pr that a par- 
ticular product state is prepared. As these individual sub-system probabilities 
P4 (i), P§(j), ..must determine actual possible measurements then the density 
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operators p^, p R , ..must correspond to possible physical states for the sub- 
systems, the sub-systems being modes or single particle states in the present 
case. But as we have seen, the possible physical states that can be prepared 
for these sub-systems are those as in Eq. (fltfl) which are a statistical mixture 
of number states with no coherences between Fock states with differing boson 
numbers, so the p^, pf, themselves satisfy the super-selection rule. 

2.9.5 Example of State that Violate Local and Global Particle Num- 
ber SSR 

Finally, an objection to applying the super-selection rule to separate modes 
based on emphasising only measurements on only one mode and its the reduced 
density operator may be raised, and suggest that . 'p R , p^ etc may be allowable 
provided that the overall reduced density operators comply with the super- 
selection rule. However, as will be seen this is not in general possible. As shown 
above, measurements on the subsystems with measurements on the other sub- 
systems discarded - are determined only from the reduced density operators 
~p A = ^2rPr7>r alone. Hence it may seem that providing the reduced density 

^A ^~B *~-G 

operators represent physical states then it does not matter if the p R , p R , p R .. 
do not. Indeed, for special cases we can find density operators 'p R that are 
unphysical even though the reduced density operator p A is physical One such 
example is where 

pf = (^(|o> A + |i> A ))(^«oU + <i| A )) Pi = 5 

P2 = (^=(|0> a -|1>a)) (^«0U-<1U)) P2 = \ (42) 
which yields 

^ = | (10)^ <0U) + | <1| A ) (43) 

This is a valid statistical mixture of two physical states for mode A, namely a 
state with no bosons and a state with one boson, even though the contributions 
~Pi and pjf are non physical states consisting of pure states that are each quan- 
tum superpositions of a zero boson state and a one boson state - in violation of 
the super-selection rule. However even a minute change in the Pr will lead to 
the reduced density operators ~p Al ~p B , .. that are non physical. In the example 
given, changes to Pi = 0.51 and P2 = 0.49 will lead to non physical contribu- 
tions |0) A (1^ and 11)^ (0^ to the reduced density operator p A . Also, as all 
the reduced density operators must represent physical states, then the sums in 
7> A = ^2 R Pr7>ri 7>b = ^2r Pr'Pri ■■■ m ust all lead to physical states. Since the 
probabilities P R depend on the preparation process that generates the mixed 
non-entangled state, and may for example depend on external parameters such 

■^A ^B 

as temperature, it would be extremely unlikely for given p R , p R , ..that all such 
sums will lead to physical states, though for special choices of the mode density 
operators and the Pr this can occur. In addition, the density operators for the 
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other modes must be chosen so that the overall density operator is consistent 
with the super-selection rule. For example in the case where there are only two 
modes, the density operators pf = |0) B (0| B and = |1) B (1| B would lead to a 
physically valid reduced density operator p B = | (|0) B (0| B ) + \ (|1) B for 
mode B, but there would be terms such as \ |0) A ® |0) B (0| B in the over- 
all density operator, and such a term involves a coherence between an N = 
state and an N — 1 state which is disallowed. Indeed, for the 'pf, and Pi, 
P2 as in Eq. (|42|) . there may be no choice for p x and /3 2 that gives rise to 
an overall physical state. In Appendix [13] the situation where 'p 1 and p 2 are 
associated with two general pure orthogonal states of the form a |0) B + (3 |1) B 
and -0* |0) B + a* \1) B with (|a| 2 + |/3| 2 ) = 1, is considered, and we find that no 
choice of a and /3 leads to an overall physical state - although again the reduced 
density operator p B = \ (|0) B (0| B ) + \ (|1) B (1| B ) is physical. 

2.9.6 Example of Global but not Local Particle Number SSR Com- 
pliant State 

However, in some cases sub-system density operators can be chosen in the con- 
text of two mode systems which comply with the global particle number SSR 
but not the local particle number SSR. Such a case is presented by Verstraete 
et al [37], 28 . The overall density operator is a statistical mixture 

P = ^(l^i> <^iI)a <8> |^!> <Vil)s + i(l^> <^IU ® 

+ {i>-l\)A ® (^-x\)B + \(\l>-i) (l>-i\)A ® \i>-i) (lU^ 

where \i/jJ = (|0) +w|l))/v / 2, with oj = The \<ipj are superposi- 

tions of zero and one boson states and consequently the local particle number 
SSR is violated by each of the sub-system density operators (^ u \)a and 
\ij)J) {iP u \)b- O n the other hand, the global particle number SSR is obeyed 
since the density operator can also be wriiten as 

P = i(|0)(0|U®|0)(0|) B + i(|l)(l|U®|l)(l|) B 

+ \(\*+) {*+\)ab (45) 

where \^/ + ) AB — (|0) A |1) B + |0) B )/\/2. This is a statistical mixture of 
N = 0,1,2 boson states. Although the expression in Eq. (l44l) is of the form 
in Eq.Q]), the subsystem density operators \ip u ) (iP u \)a an d \ip u ) (V'wI)b do n °t 
comply with the local particle number SSR, so this paper the state would be 
regarded as entangled. However, Verstraete et al [27], [28] regard it as separable. 
They would call it separable but nonlocal. 

2.9.7 General Form of Non-Entangled States 

To summarise: basically the sub-systems are single modes that the identical 
bosons can occupy, the super-selection rule for identical bosons, massive or oth- 
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erwise, prohibits states which are coherent superpositions of states with different 
numbers of bosons, and the only physically allowable p R , p R , ..for the separate 
mode sub-systems that are themselves compatiible with the local particle num- 
ber SSR are allowed. For single mode sub-systems these can be written as 
statistical mixtures of states with definite numbers of bosons in the form 

Pr = J2 P t M (n A \ P B R =J2 P ns \™b) (n B \ .. (46) 

TLA n B 

However, in cases where the sub-systems are pairs of modes the density 
operators p^, p^, ..for the separate sub-systems are still required to conform 
to the symmetrisation principle and the super-selection rule. The forms for 
Pri Pr> ■■ are now °f course more complex, as entanglement within the pairs 
of modes Ai, A2 associated with sub-system A, the pairs of modes B\, B2 
associated with sub-system B, etc is now possible within the definition for the 
general non-entangled state Eq. JT]) for these pairs of modes. Within each pair 
of modes Ai, Ai statistical mixtures of states with differing total numbers Na 
bosons in the two modes are possible and the sub-system density operators are 
based on states of the form given in Eq. I|29p. We have 

N A 

\®Na) A = E C A *(N A , k) \k) M ® \N A - k) M 

00 

PR = EE P ^I^>a(^Ia (47) 

N A =0 $ 

with analogous expressions for the density operators etc for the other pairs 
of modes. Note that |$at a )^ only involves quantum superpositions of states 
with the same total number of bosons Na. The expression (|157|) in SubSection 
15.31 is of this form. 



2.10 Two Mode Coherent State Mixture 

To further illustrate some of the points made about super-selection rules - local 
and global - it is useful to consider a specific case also presented by Verstraete 
et al [27) . |28) . This mixture of two mode coherent states is represented by the 
two mode density operator 

P = J*7jT K a ) ( a >«l 

= S^{\a)(a\) A 9(\a)(a\) B (48) 

where \a) c is a one mode coherent state for mode C = A 7 B with a = \a\ exp(— i8), 
and modes A, B are associated with bosonic annihilation operators a, b. The 
magnitude \a\ is fixed. 
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This density operator appears to be that for a non-entangled state of modes 
A, B in the form 

P = ZPrPb®Pr (49) 

R 

with J2PR -> / I! and p£ -> (|a) (a \) A and p H — > (|a)(a|) B . However al- 

R 

though this choice of p B satisfy the Hermitiancy, unit trace, positivity fea- 
tures they do not conform to the requirement of satisfying the (local) sub-system 
boson number super- selection rule. From Eq. (|48|) we have 



(n\(\a)(a\)\m) A = exp(-|a| 2 ) 



a" (a)* 



n! V mi 



(p|(|a)<a|)|«) fl = exp(-H 2 )^=M^ (50) 

so clearly for each of the separate modes there are coherences between Fock 
states with differing boson occupation numbers. In the approach in the present 
paper the density operator in Eq. (|48|) does not represent a non-entangled state. 
However, in the papers of Verstraete et al [37], [IB], Hillery et al [25], [30] and 
others it would represent an allowable non-entangled (separable) state. Indeed, 
Verstraete et al [27] specifically state " .., this state is obviously separable, though 
the states \a) are incompatible with the (local) super-selection rule." . Verstraete 
et al [27] introduce the state defined in Eq. (j48|) as an example of a state that 
is separable (in their terms) but which cannot be prepared locally, because it is 
incompatible with the local particle number super-selection rule. 

The mixture of two mode coherent states does of course satisfy the total or 
global boson number super-selection rule. The matrix elements between two 
mode Fock states are 

{(n\ A ®(p\ B )p(\m) A ®\q) B ) = exp(-2H 2 ) ^ + _ / ^ exp(-t(n - m + p - q)6) 



\a\ p+q 

exp(-2|a[ 2 ) '^L 1 S n+p>m+q (51) 

i! vpWq- 



These overall matrix elements are zero unless n +p = m + q, showing that there 
are no coherences between two mode Fock states where the total boson number 
differs. The mixture of two mode coherent states has the interesting feature of 
providing an example of a two mode state which satisfies the global but not the 
local super-selection rule. 

The reduced density operators for modes A, B are 



Pa = I g - (I") H)a Pb = J 7^ (I") W)b 



and a straightforward calculation gives 



\a\ 2n 



PA = exp(-H 2 ) E ~ir d n > HU Pb = exp(-|a| 2 ) ^ (\p) (pi), 
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which are statistical mixtures of Fock states with the expected Poisson distribu- 
tion associated with coherent states. This shows that the reduced density oper- 
ators are consistent with the separate mode local super-selection rule, whereas 
the density operators ~p R = (|a) (a\) A , = (\a) (a\) B are not . Later we will 
revisit this example in the context of entanglement tests. 

2.11 Two Mode System - Coherence Terms 

The general non-entangled state for modes a and b is given by 

p = J2P*pi®P% (52) 

R 

and as a consequence of the requirement that p^ and p^ are physical states for 
modes a and b satisying the super-selection rule, it follows that 

((o)") = Tr(p£(ar)=0 ((&)») a = Tr(p£@) n ) = 
((b) m ) b = Tr(p B R (p) m ) = Q ((6t)™) 6 = Tr(^(?n = 

(53) 

Thus coherence terms are zero. As we will see these results will limit spin 
squeezing to entangled states of modes a and b. Note that similar results also 
apply when non-entangled states for the original modes c and d are considered 

- ((c)") c = 0, etc.. 

2.12 Two Sub-Systems of Pairs of Modes - Coherence 
Terms 

In this case the general non-entangled state where A and B are pairs of modes 

- A\, A2 associated with sub-system A, and modes B\ , B2 associated with sub- 
system B, the overall density operator is of the form ([71|) . with C —> A, D — >• B, 
whilst the sub-system density operators are of the forms given in (|47l) . In this 
case we now have in general 

((2,)")^ - TripifrD^O {(at) n ) A = Tr(p£(atr)?Q 
((bj) m ) B = TrCp B R {b 3 r)^Q ((b]) m ) B =TrCp B R {V j ) m )^{) 

i,j = 1,2 (54) 

so unlike the case where the two sub-systems are single modes, there are non- 
zero coherences when they are pairs of modes. 
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3 Spin Squeezing 



The basic concept of spin squeezing was first introduced by Kitagawa and Ueda 
[40] for general spin systems. These include cases based on two mode systems, 
such as may occur both for optical fields and for Bose-Einstein condensates. 
Though focused on systems of massive identical bosons, the treatment in this 
paper also applies to photons though details will differ. 

3.1 Spin Operators, Bloch Vector and Covariance Matrix 

3.1.1 Spin Operators 

For two mode systems with mode annihilation operators a, b associated with 
the two single particle states \4> a ), and where the non-zero bosonic com- 
mutation rules are [e 1 e'] — 1 (e = a or b), Schwinger spin angular momentum 
operators (£ = x, y, z) are defined as 

S x = (fe 1 a + a^b) /2 S y = (6 f a - a f 6) /2i S z = fib - a f a) /2 (55) 

and which satisfy the commutation rules [S% , 5 M ] = ie^xSx for angular mo- 
mentum operators. For bosons the square of the angular momentum operators 
is given by Si + S 2 + S 2 = (jV/2)(jV/2 + 1), where N = {tt+a^a) is the boson 
total number operator, those for the separate modes being n e = e^e (e — a 
or b) . The Schwinger spin operators are the second quantization form of sym- 
metrized one body operators S x = (0 o (i)l + l0aW) (0t(*)l)/ 2 5 &v = 
EMM (4> a (t)\-\<t> a (i)) (tf> b (f)\)/2i ; S, = £,(|&(i)) (0aWI)/2 
, where the sum i is over the identical bosonic particles. In the case of the two 
mode EM field the spin angular momentum operators are related to the Stokes 
parameters. 

3.1.2 Bloch Vector and Covariance Matrix 

If the density operator for the overall system is p then expectation values of the 
three spin operators (s^) = Tr(pS^) (£ = x, y, z) define the Bloch vector. Spin 

squeezing is related to the fluctuation operators AS% = — i n terms of 

which a real, symmetric covariance matrix C^S^yS^) (£, /x = x,y,z) is defined 
HSU, [56] via 

= (s £ S M +S^)/2-(s e )(3 M ) (56) 

and whose diagonal elements C(S^, Sf ) = ^AS^ 2 ^ gives the variance for the fluc- 
tuation operators. The variances for the spin operators satisfy the three Heisen- 
berg uncertainty principle reations (AS X A {ASy 2 ) > || /§ z \ | 2 ; (aS v 2 \ (aS z 2 
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|| (^) I 2 ; (AS* 2 } (^ASx 2 ^ > || (Sy) | 2 j an d spin squeezing is usually defined 

via conditions such as /AS X 2 ^ < || ^S*^ | with (^AS y 2 ^ > || |, for 

being squeezed compared to S y and so on. However this definition is unsatisfac- 
tory since it ignores the presence of the off-diagonal elements of the covariance 
matrix, so a better definition is required. 

3.2 New Spin Operators and Principal Spin Fluctuations 

The covariance matrix has real, non-negative eigenvalues and can be diago- 
nalised via an orthogonal rotation matrix M(—a, —{3, —7) that defines new spin 
angular momentum operators (£ = x,y,z) via 

and where 

C{JiX) = 5^M CA (-a,- ) 8,-7)C(^A,^)M^(-a J -)8,-7) 

X9 

~* 2 



= ^(AV) (58) 

is the covariance matrix for the new spin angular momentum operators Jj (£ = 
x, y, z), and which is diagonal with the diagonal elements ^AJ K 2 ^ , (^AJy 2 ^ and 

^A J z 2 ^ giving the so-called principal spin fluctuations. The matrix M(a, /3, 7) 

is parameterised in terms of three Euler angles a, /3, 7 and is given in |65j (see 
Eq. (4.43)). 

The Bloch vector and spin fluctuations are illustrated in Figure 1. In Fig 1 
the Bloch vector and spin fluctuation ellipsoid is shown in terms of the original 
spin operators S% (£ — x,y,z) 



Figure 1 near here. 



3.3 Spin Squeezing for New Spin Operators 

3.3.1 Heisenberg Uncertainty Principle and Spin Squeezing 

Since the new spin operators also satisfy Heisenberg uncertainty principle rela- 
tionships 

(ajL 2 )(aj, 2 ) > i|(j,)| 2 
(aj, 2 )(aj z 2 ) > \\{j x )\ 2 

(AJ 2 2 )(AJ K 2 ) > i|(j,)| 2 (59) 



39 



spin squeezing will now be defined via condtions such as 

(aj x 2 ) < ~\(j z )\and (AJ, 2 )>i|(j 2 )| 
(AJ, 2 ) < ±\(7 x )\and (AJ Z 2 ) > 

(AJ Z 2 ) < \\(j y )\and (AJL 2 )>i|(j y )| (60) 

for J^being squeezed compared to J y , and so on. By convention we may choose 
^AJ X 2 ^) < (^AJ y 2 ^ < ^AJ 2 2 ^), so the primary spin operator of interest will be 

J x since this has the smallest fluctuation. Note that here we have chosen prin- 
cipal spin fluctuations, but of course the last Heisenberg uncertainty relations 
apply for any new choice of rotated spin operators - as occurs in the next part 
of this section. 



3.3.2 Alternative Spin Squeezing Criteria 

Other criteria for spin squeezing are also used, for example in the article by 
Wineland et al [66 . To focus on spin squeezing for J x compared to any orthog- 
onal spin operators we can combine the first and third Heisenberg uncertainty 
principle relationships to give 

(AJ X 2 ) ((AJ, 2 ) + (AJ Z 2 )) > \ (l (j y ) | 2 + I (j z ) | 2 ) (61) 

Then we may define two new spin operators via 

J± i = cos 8 J y + sin 8 J z Jxi = — sin 8 J y + cos 8 J z (62) 

where 8 corresponds to a rotation angle in the yz plane, and which satisfy the 
standard angular momentun commutation rules [ J± i , Jj_ 2] — iJx, [Jx 2 , Jx } — 
iJxi, [Jx ,Jxi] = iJx2- It is straightforward to show that ^A J y 2 ^ + {^A J z 2 ^ — 

(AJ ±1 2 ) + (AJ ±2 2 ) and | | 2 + | (j ±2 ) | 2 = | (J y ) | 2 + | (7 Z ) | 2 so that 

(AJ X 2 ) ((AJ ±1 2 ) + (AJ i2 2 )) > \ (| (Jxx) | 2 + I (Jx») I 2 ) (63) 

so that spin squeezing for ^compared to any two orthogonal spin operators 
such as Jj_ 1 or Ji 2 would be defined as 

(aj, 2 ) < y(\(j ±1 )F + \(j ±2 )\') 

and 

(AJ ±1 2 ) + (AJ ±2 2 ) > ^(|(Jli)| 2 + |(Jl 2 )| 2 ) (64) 
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This criterion would apply however the choice of rotation matrix M(—a, —(3, —7) 
is made, so A J x does not have to correspond to the principal spin fluctuation 
with the smallest variance though obviously such a choice is preferable over 
some arbitrary set of new spin operators. For spin squeezing in ^AJ 2 2 ^ we 
require 

(ajA 1 1 

f= 7 ,~ s L v x < j ~T7 (65) 



N 



The last step is an approximation based on the assumption that the Bloch 
vector lies in the yz plane and close to the Bloch sphere, this situation being 
the most conducive to detecting the fluctuation (AJ a 2 ). In this situation 



JJ\ (-/li) I 2 + I (^2) I 2 ) is approximately N/2. The condition £ 2 < 1/N is 
sometimes taken as the condition for spin squeezing [6 7) , but it should be noted 



that this is approximate and Eq. (1641) gives the correct expression. 
3.3.3 Planar Spin Squeezing 

A special case of recent interest is that referred to as planar squeezing [55] in 
which the Bloch vector for a suitable choice of spin operators lies in a plane 
and along one of the axes. If this plane is chosen to be the xy plane and the 
x axis is chosen then /j 2 j = and (jy*j — 0, resulting in only one Heisen- 
berg uncertainty principle relationship where the right side is non-zero, namely 
AJ y 2 ^ (a.J 2 ^ > || (Z^ I 2 . Combining this with ^AJ X 2 ^ (AJ y 2 ^ > gives 

AJ y 2 ^ + (^AJ x 2 Jj (^AJ 2 ^j > || (j^ | 2 . So the total spin fluctuation in 

the xy plane defined as ^AJ|| 2 ^} = i^AJy 2 ^ + (^AJ X 2 ^ will be squeezed com- 
pared to the total spin fluctuation perpendicular to the xy plane given by 



(( 



AJ ± 2 = AJ Z 2 if 



(Ajj| 2 ) < l -\{j x )\and (AJ ± 2 ) > (<><>') 



By minimising ( AJy ) whilst satisfying the constraints \ JzJ = \^v/ = a spin 
squeezed state is found that satisfies with ( AJn 2 ) ~ J 2 / 3 , ( AJj_ 2 ) ~ J 4 / 3 , 



Jxj I ~ J for large J = N "/2 [68]. The Bloch vector is on the Bloch sphere 
and condition (l65l) is also satisfied. 
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3.4 Rotation Operators and New Modes 



3.4.1 Rotation Operators 

The new spin operators are also related to the original spin operators via a 
unitary rotation operator i?(a,/3, 7) parameterised in terms of Euler angles so 
that 

J £ = R(a,f3, 1 )S i R(a,(3^)- 1 (67) 

where 

R{a,P,i) = R z (a)Ry(j3)R z (ry) (68) 

with R^(4>) — exp(i(f>S^) describing a rotation about the £ axis anticlockwise 
through an angle <f>. Details for the rotation operators and matrices are set out 
in [56] . Note that Eq. ([67]) specifies a rotation of the vector spin operator rather 
than a rotation of the axes, so (£ = x, y, z) are the components of the rotated 
vector spin operator with respect to the original axes. 



3.4.2 New Mode Operators 

We can also see that the new spin operators are related to new mode operators 
c and d via 

J x = {d)c + ^d)/2 J y = (cftc-$d)/2i J z = (d)d~c^)/2 (69) 
where 

c = R{a,P,i)aR{a,P,i)- 1 d = R(a, /3,j)b R(a, (S.t) -1 (70) 

For the bosonic case a straight-forward calculation gives the new mode op- 
erators as 

A ■ n ( A A ■ • A , 1 • si. 

c = exp(-i7) I cos(-) exp(-ia)a + sm(-) exp(--ta)b 
d = exp(-ii7) ( -sin(-) exp(^ia) a + cos(-) exp(~ia)b 



(71) 

and it is easy to then check that c and d satisfy the expected non-zero bosonic 
commutation rules are [e, e l \ = 1 (e" = c or d) and that the total boson number 
operator is N = (dtd+tfc). As N is invariant under unitary rotation operators 
it follows that Jl + Jl + Jl = {N/2){N/2 + 1). 
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3.4.3 New Modes 

The new mode operators correspond to new single particle states \<f> c ), \4> d ) 
where 

\<j) c ) = exp(-ii 7 ) ^cos(^) exp(~ia) \<f> a ) +sin(^) exp{^ia) \cj> b ] 
\<f>d) = ex P(^7) f-sin(^) exp(-iia) |</>J +cos(^) exp(iia) |0 b 

'(72) 

These are two orthonormal quantum superpositions of the original single particle 
states \4> a }, \4>b)i an d as such represent an alternative choice of modes that could 
be realised experimentally. 

Eqs. (|71[1 can be inverted to give the old mode operators via 

a = cxp(-iia) (cos(^) exp(-ii7)c-sin(^) exp(+i«7)d 



exp(+— ia) ( sin( — ) exp(— ij) c + cos(— ) exp( — ij) d 



(73) 



3.5 Two Mode System - Coherence Terms 

For our two-mode case we have also seen that the original choice of modes with 
annihilation operators a and b may be replaced by new modes with annihilation 
operators c and d. Since the new modes are associated with new spin operators 
Ji (£ = x i Hi z ) f° r which the covariance matrix is diagonal and where the 
diagonal elements give the variances that are relevant for the definition of spin 
squeezing, it is therefore more relevant to consider entanglement for he case 
where the sub-systems are modes c and d, rather than a and b. Consequently 
the general non-entangled state for modes c and d is given by 

R 

and as a consequence of the requirement that p R and p R are physical states for 
modes c and d satisying the super-selection rule, it follows that 

((cD c = Tr{p C R {c) n )=0 (@) n ) c = Tr(j%@) n )=0 
((d) m ) d = Tr(p%(dy n ) = ( K $r) d = Tr(p%($D=0 

(75) 

Thus coherence terms are zero. As we will see these results will limit spin 
squeezing to entangled states of modes c and d. 
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3.6 Quantum Correlation Functions and Measurements 

Finally, we note that the principal spin fluctuations can be related to quantum 
correlation functions. For example, it is easy to show that 

(AX 2 ) = i(((?) a (8) a ) + ((gt)W)+2(?ctoa) + (?^ + <^) 
-\ (((^) 2 + ((ctd) 2 + 2 ((?c) ((ctrf)) (76) 

showing that (^AJ X 2 ^ is related to various first and second order quantum cor- 
relation functions. These can be measured experimentally and are given the- 
oretically in terms of phase space integrals involving distribution functions to 
represent the density operator and phase space variables to represent the mode 
annihilation, creation operators. 
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4 Spin Squeezing as a Test for Entanglement 



With the general non-entangled state now required to be such that the density 
operators for the individual sub-systems must represent physical states and 
conform to the super-selection rule, the consequential link between entanglement 
in two mode bosonic systems and spin squeezing can now be established. We 
first consider spin squeezing for the principal spin operators J x , J yi J z and 
entangled states of the related new modes c, d and then spin squeezing for the 
original spin operators S x , S y , S z and entangled states of the original modes a, 
b. Examples of entangled states that are not spin squeezed and states that are 
not entangled nor spin squeezed for one choice of mode sub-systems, but are 
entangled and spin squeezed for another choice are then presented. 

4.1 Spin Squeezing Requires Entanglement 

Firstly, the variance for a Hermitian operator in a mixed state 

R 

is always greater than or equal to the the average of the variances for the separate 
components 

^An 2 )>Y,PR(^ 2 ) (78) 

R R 

h 2 ^ = Tr{pAn 2 ) with Aft = and^Aft 2 ^ = Tr(p R AQ R 2 ) 

with AQ,r = ft — (QJ . The proof is straight-forward and given in Ref. [55] . 

4.1.1 Cases of J x and J y 

Next we calculate ^AJ 2 ^} , ^AJ 2 ^) and > (^) ' ("^} ^ or ^ e case 

where 'p R = 'p R <E> 'p R - From Eqs. (|69|) we find that 



i((?) 2 (?) 2 + dJdcc 1 + c*cd$ + {d) 2 (Pf) 



4 

J 2 = -\{0) 2 {c) 2 -d^dcd^ -?^c~d$ + {df{^) 2 ) (79) 

so that on taking the trace with ~p R and using Eqs. (|75[) we get after applying 
the commutation rules [e, &] = 1 (e = c or d) 



J y> R 



= l({^ R + (^R) + l((^R($d) R ) 

= \i{^ H + (^>n) + \((^>H(^ R ) (80) 
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As we also have 

'j a 



(jy) R = 5^)^H-a(^)=0 (81) 
using Eqs. (|75|) and we see finally that the variances are 

\({$2) R + (^c) R ) + i((ctc> fl ($d) R ) (82) 



(AJ 



and therefore from Eq. (|78[) 

( A ^) > E^lK^+^ + ^^H^J (83) 



Now 



^) = E^?((^ r -( 2+2 )h)) ( 84 ) 



so that 

\\ (J z ) I < £^((#3) -<#3k))I < E P «l((^ +<^>fl)) (85) 



and thus for any non-entangled state of modes c and d 

(A^)-i|(l.)l 





Similar final steps show that (AJ 2 ,) — || w x) | > for any non-entangled 



v i 

state of modes c and <i. 

This shows that for the general non-entangled state with modes c and <i as the 
sub-systems, the variances for two of the principal spin fluctuations /AJ^\ and 

(AJy^ arc &o£/i greater than || ^ J z \ |, and hence there is no spin squeezing for 



^ or Jy. Note that as | ( J y ) \ — 0, the quantity J I | ( J± i ) | 2 + | ( Jj_ 



-2/I 2 



4G 



is the same as | \ J zj |, so the alternative criterion in Eq. is the same as 

that in Eq. (IBUl) which is used here. 

We can extend the above to obtain further inequalities for the non-entangled 
state. Using Eq. flSTJl 



j x ) = x> v*)„ = W = E p « (^) B = ( 8? ) 

R 

it is easy to see that 



R \ I \ I R 

R R 



Ml) ~ \\ (Jy) I > (AJ») - \\ (j x ) | > (88) 

for any non-entangled state of modes c and d. This completes the set of inequal- 
ities for the variances of J x and J y . 

4.1.2 Case of J z 

For the other principal spin fluctuation we find that 



(89) 



so that using ([75]) 

y^Efti^^-WJ 1 ) +((^-<^>h) 2 ) c (90) 



a 4 \ v \ <* 



From Eq. ([87]) it follows that 



f 2 

> (91) 
Similarly (aJ^ - ±| | > 0. 

4.1.3 No Spin Squeezing for Non-Entangled States 

So overall, we have for the general non-entangled state of modes c and d 



A./,-) ■ i|^|ond (AJ v 2 )>~|<./ : 
1, /-?\ , . /. <?„\ 1 



i|(j v )|ond (AJL 2 )>i| 



AJ Z 2 ) > £|(j v )|and (AJ K 2 )>i|(J B )| (92) 
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Note that the last two pairs of inequalities are trivially true for the general 
non-entangled state, since (^Jx\ = (^v) = ^* This overall result tells us that 

for any non-entangled state of modes c and d we do not have J x being squeezed 
compared to J y (or vice- versa), J y being squeezed compared to J z (or vice- 
versa), J z being squeezed compared to J x (or vice- versa). That is, there is no 
spin squeezing for the non-entangled state! 

4.1.4 Spin Squeezing Tests for Entanglement 

The key value of these results is the spin squeezing test for entanglement - if 
for a given state we find that 

// (AJ*)<±|(j,)| or (Aj2)<I|(j z )| (93) 

or 

// (AJl)<±\(j x )\ or (AJl)<±\(j x )\ (94) 



If ( K AJ 2 2 )<±\(j v )\ or (AJl)<i|(j„)| (!»•-„ 



then the state must be entangled. Thus we only need to have spin squeezing in 
any of the J x , J y or J z to demonstrate entanglement. No particular component 

need be singled out. Note that one cannot have both ^AJ^ < || /j^ 

and^AJ^\ < || (^J z\ \ etc. due to the Heisenberg uncertainty principle. 
It is then straightforward to show that 



or 



or 



If (*Jl)<^\{J y ±l)\ 2 + \{Jlz)\ 2 07) 



// (Aj;}<i0(jl 1 }|2 + |(ji 2 )p (98) 

that is, if J x , J y or Jj are squeezed compared to any of their two orthogonal 
spin components - then the state must be entangled. Again we only need to 
have spin squeezing in any of the J x , J y or J z compared to any of their two 
orthogonal spin components to demonstrate entanglement. 



4.1.5 Inequality for | ( J 

Of the results for a non- entangled physical state for modes c and d we will later 
find it particularly important to consider the first of I 



AJi)>\\(j z )\ and (AJ^>\\{J : )\ (!.<)) 
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This is because we can show that for any quantum state 

I (Jz) I = I (^(n d ~ n c )) | < \{\ (n d ) \ + \ (n c ) |) = \ (jv) (100) 

there is an inequality involving | 1^3 z ^ | and the mean number of bosons ^-^^ 
in the two mode system. Note that there may be entangled states for which 
AJ^ and I^AJ^j are both greater than || l^J z \ |, since all that has been 

proven is that for non-entangled states we must have both (^AJ^J > i| (jz\ I 

and (aJ*) > i|(j z )|. 

Hence we may conclude that spin squeezing in either of the principal spin 
fluctuations J x , J y or J z requires the quantum state to be entangled for the 
modes c and d as the sub-systems, these modes being associated with the prin- 
cipal spin fluctuations via Eq. (|69|) . Although finding spin squeezing tells us 
that the state is entangled, there are however no simple relationships between 
the measures of entanglement and those of spin squeezing, so the linkage is 
essentially a qualitative one. For general quantum states, measures of entangle- 
ment for the specific situation of two sub-systems (bi-partite entanglement) are 
reviewed in [15] . 



4.2 Spin Squeezing requires Entanglement for Original 
Modes 

It is also of some interest to consider spin squeezing for the original spin opera- 
tors S x , S y , S z with the original modes a and b as the sub-systems, even though 
these spin operators are in general associated with a non-diagonal covariance 
matrix and the concept of spin squeezing is rather problematic in view of princi- 
pal spin fluctuations not being involved. In this case the general non-entangled 
state for the original modes is given by 

P = Y, P *pi®PR (101) 

R 

with the p R and p R representing physical states for modes a and b, and where 
results analogous to Eqs. (|75p apply. The same treatment applies as for spin 
operators J x , J v , J z with the modes c and d as the sub-systems and leads to 
the result for a non-entangled state of modes a and b 

(ASI)>±\(S Z )\ and (AS 2 y )>\\(s z ) (102) 

showing that neither S x or S y is spin squeezed for the general non-entangled 
state for modes a and b given in Eq. (|101 1) . We also have 

(s x )=J2Pr(s*) r = (Sy) = Y. p n(Sy) R = ° ( 103 ) 
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so all the results analogous to Eqs. (li?2"j) also follow. Hence we may also conclude 
that spin squeezing in any of the original spin fluctuations requires the quantum 
state to be entangled when the original modes a and b are the sub-systems. Thus 
the entanglement test is 



// (ASi)<l\(.S 



<A^)<iK.s 



(104) 



// (ASi)<-\(S 



or 



// (ASi)<-\{Sy 



or (ASi)<-\(S x 



ASt) < ~\(S 



(105) 
(106) 



then we have an entangled state for the original modes a and b. 

Hence we have seen that spin squeezing - either of the new or original spin 
operators requires entanglement of the new or original modes - the question then 
is: Does entanglement automatically lead to spin squeezing? The answer is no, 
since cases where the quantum state is entangled but not spin squeezed can be 
found. Thus in general, spin squeezing and entanglement are not equivalent.- 
they do not occur together for all states. Spin squeezing is a sufficient condition 
for entanglement, it is not a necessary condition. 

4.3 Entangled States that are Non Spin-Squeezed 

One such example is the generalised N boson NOON state defined as 



P 
|d>> 



= mm 



cos ( 



(ct) 



/V 



10) 



N N 
7~~2 



. s in0^L|O) 



— , - — + sm 



AT! 
N 

7' 



N 

"7 



(107) 



which is an entangled state for modes c and d in all cases except where cos 9 or 
sin 9 is zero. In the last form the state is expressed in terms of the eigenstates 
for (_J) 2 and J 2 , as detailed in [56] . 
A straight-forward calculation gives 



J X 

AJl 





N 

T 



Jy)=0 



A J 



N 

■ I ) = -— cos 20 



N 

T 



AJ: 



(1 - cos 2 29) (108) 



for N > 1, so that using the criteria for spin squeezing given in Eq. (|60p we see 
that as ^AJ 2 ,^} — || (^J zj I > 5 etc ! an d hence spin squeezing does not occur 
for this entangled state. 
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4.4 Non-Entangled States that are Non Spin Squeezed 

Of course from the previous section any non entangled state is definitely not 
spin squeezed. A specific example illustrating this is the N boson binomial state 
given by 

p = |$)($| 
(-cf) N 

m = -^=-\°) ( 109 ) 

where c and d are given by Eqs. (|7T|) with Euler angles a = —n + x, = —2$ 
and 7 = — it, we find that 

c = — cos9 exp(—ix)a — sin8 exp(——ix) b = —ai 

d = s'm8 exp(—ix)a — cos9 exp(--ix) b — —a,2 (HO) 

where the mode operators a± and a,2 are as defined in [56] (see Eqs. (53) therein). 
The new spin angular momentum operators (£ = x, y, z) are the same as those 
defined in [56] (see Eqs. (64) therein) and in [56] it has been shown (see Eq. 
(60) therein) for the same binomial state as in (|109|) that 

J.) = o (•?„) = « (^) = -y 
(ATI) = f (aJJ)-o (in) 

(see Eqs. (162) and (176) therein). Hence the binomial state is not spin squeezed 

since ^AJ^\ = ^AJ 2 ^ = ^| z ^ |. It is of course a minimum uncertainty 

state with spin fluctuations at the standard quantum limit. Clearly, it is a non- 
entangled state for modes c and d , being the product of a number state for 
mode c with the vacuum state for mode d. 

Note that from the point of view of the original modes a and b, this is an 
entangled state, so the question is: Is it a spin squeezed state with respect to 
the original spin operators S$ (£ = x, y, z) ? The Bloch vector and variances for 
this binomial state are given in |56) (see Eq. (163) in the main paper and Eq. 
(410) in the Appendix). The results include: 

~ \ N 

S Z J = -yCOS26» 

ASl) = — (cos 2 28 cos 2 X + sin 2 X ) (a^ 2 ) = ^(cos 2 26* sin 2 X + cos 2 X ) 

(112) 

This gives (AS 2 ^ (AS 2 ) - i| (S z ^ \ 2 = ^7V 2 (cos 2 28 - l) 2 cos 2 X sin 2 X > 
as required for the Heisenberg uncertainty principle. With x = we have 
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(AS*) = f cos 2 2(9 and (aS 2 ,^ = f, whilst || (s z \ | = f|cos20|. As 
(^AS 2 ^ < i| Wz) I there is spin squeezing in S x for this entangled state of 

modes a and 6, though not of course for the new spin operator J x since this 
state is non-entangled for modes c and d. This example illustrates the need to 
carefully define spin squeezing and entanglement in terms of related sets of spin 
operators and modes. The same state is entangled with respect to one choice 
of modes - and spin squeezing occurs, whilst it is non-entangled with respect to 
another set of modes - and no spin squeezing occurs. 

To summarise - with a physically based definition of non-entangled states for 
bosonic systems with two modes (related to the principal spin operators that 
have a diagonal covariance matrix) being the sub-systems and with a criterion 
for spin squeezing that focuses on these principal spin fluctuations, it seen that 
whilst non-entangled states are never spin squeezed and therefore although en- 
tanglement is a necessary condition for spin squeezing, it is not a sufficient one. 
There are entangled states that are not spin squeezed. Furthermore, as there 
is no simple quantitative links between measures of spin squeezing and those 
for entanglement, the mere presence of spin squeezing only demonstrates the 
qualitative result that the quantum state is entangled. Nevertheless, for high 
precision measurements based on spin operators where the primary emphasis 
is on how much spin squeezing can be achieved, knowing that entangled states 
are needed provides an impetus for studying such states and how they might be 
produced. 



4.5 Entangled States that are Spin Squeezed 
4.5.1 Relative Phase Eigenstate 

As an example of an entangled state that is spin squeezed we consider the 
relative phase eigenstate \^;,0 P ) for a two mode system in which there are N 
bosons. For modes with annihilation operators a, b the relative phase eigenstate 
is defined as 

±1^) = —^= Y exp(ifc^)4= — 10) (113) 

where the relative phase 6 p = p(2ir/(N+l)) withp = -N/2, -N/2+1, +N/2, 
is an eigenvalue of the relative phase Hermitian operator of the type introduced 
by Barnett and Pegg [70] (see [56] and references therein) . Note that the eigen- 
values form a quasi-continuum over the range — it to +n, with a small separation 
between neighboring phases of 0(1/N). The relative phase state is consistent 
with the super-selection rule and is an entangled state for modes a, b. The 
Bloch vector for spin operators S x , S y , S z is given by (see [56] ) 

S x ) = N^cos6 p (S v ) = -N^mx6 p (s,) = (114) 

but the covariance matrix (see Eq. (178) in [56) ) is non-diagonal. 
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4.5.2 New Spin Operators 



It is more instructive to consider spin squeezing in terms of new spin operators 
J x , J y , J z for which the covariance matrix is diagonal. The new spin operators 
are related to the original spin operators via 

Jx = S z 

J y = sin p S x + cos Op S y 

J z — — cos Op S x + sin Op S y (115) 

corresponding to the transformation in Eq. (|57p with Euler angles a = —tt + O p , 
j3 = — 7r/2 and 7 = — ir. 



4.5.3 Bloch Vector and Covariance Matrix 



The Bloch vector and covariance matrix for spin operators J x , J y , J z are given 
by (see Eqs. (180), (181) in 56^ - note that the C(J y ,J y ) element is incorrect 
in Eq. (181)) 



J x) = 



and 



C{Jx, Jx) 

C(jy, J X ) 
C{J Z , Jy) 



J-N 2 

12 iv 






Jy =0 



C{J X ,Jy) 

C{,Jyi Jy) 

C(J Z> Jy) 






J z ) = -Af- 



file) 



C{Jx,Jz) 

C(jy, J Z ) 

C(J Z ,J Z ) 



llniV 



\6 64 J 



N 2 



iV> 1 



(117) 



With 



AJl) = ±N 2 



\\xxN and (AJ?) = 



the only non-zero Bloch vector component being ( J 



that (AJl)(AJ 2 ) 2 



± - j N 2 and 
-N ^ it is easy to see 
J z j 1 2 as required by the Heisenberg Uncertainty 

Principle. The principal spin fluctuations in both J x and J z are comparable 
to the length of the Bloch vector and no spin squeezing occurs in either of 
these components. However, spin squeezing occurs in that J y is squeezed with 

respect to J x - ^AJ 2 ^} only increases as |lnA^ whilst i| I increases as 

YqN for large N. Hence the relative phase state satisfies the test in Eq. (|95p 

to demonstate entanglement for modes c, d. 



4.5.4 New Modes Operators 

To confirm that the relative phase state is in fact an entangled state for modes 
c, d as well as for the original modes a, b, we note that the new mode operators 
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c, d are given in in Eq. flTTj) with Euler angles a = —it + 9 p , /3 = —tt/2 and 
7 = — 7r. The old mode operators are given in Eq. (|T3[) and with these Euler 
angles we have 



exp(-i0p)-^= (c- dj 
exp(~i6 p )-j= (c + d) 



(118) 

This enables us to write the phase state in terms of new mode operators c, d as 

. x a, W/2 AT/4-AS/2 7V/4+fe/2 

y^) = ,/ /v-Li ( 77^ ) X! X! 

/ yiy-rx \V^/ k=-N/2r=-N/4+k/2s=-N/4-k/2 

1 1 



r(-l) 



N/i-k/2+r 



y/(N/2-k)\ y/(N/2 + k)\ 

x (N/2 - k)\ (N/2 + k)l 

(N/4 - k/2 - r)\(N/4 ~ k/2 + r)\ (N/4 + k/2 ~ s)\(N/4 + k/2 + s)l 

x{d r ) N/2-(r+s) ^ )N /2 + (r+s) | Q) 

(119) 

We see that the expression does not depend explicitly on the relative phase 6 P 

when written in terms of the new mode unnormalised Fock states ($) N / 2 ~( r+s> > ( ( ft) JV / 2 +( r + s ) |0). 

This pure state is a linear combination of product states of the form | N/2 — m) c <E) 

\N/2 + m) d for various m - each of which is an N boson state and an eigenstate 

for J z with eigenvalue m, and therefore is an entangled state for modes c, d 

which is compatible with the global super-selection rule. Note that there cannot 

just be a single term m involved, otherwise the variance for J z would be zero 

instead of — N 2 . We will return to the relative phase state again in 

SubSection 15.11 

4.6 Bloch Vector Entanglement Test 

We have seen for the general non-entangled states of modes c and d or of modes 
a and b that 



J x ) = {Jy)=0 (120) 

[S x ) = (S y )=0 (121) 
From Eqs. and (|55|) these results are equivalent to 

[d$^ = (c<fr}=0 (122) 

'btf) = (abA = (123) 
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Hence we find further tests for entangled states of modes c and d or of modes 
a and b 

\(dc^\ 2 > |(cd^| 2 >0 (124) 
\(ba^\ 2 > |(a6 f )| 2 >0 (125) 

As we will see in Section [5j these tests are particular cases with m = n = 1 
of the simpler entanglement test in Eq. (|149[) that applies for the situation in 
the present paper where non-entangled states are required to satisfy the super- 
selection rule. 



55 



5 Other Proposed Tests for Entanglement 



There are a number of inequalities involving not only the variances of the spin 
operators but also other quantities, that have been derived for testing whether 
a state for a system of identical bosons is entangled. These are not always 
associated with criteria for spin squeezing. Some of these are based on the 
implicit assumption that the density operators p R , p R in the expression for a 
non-entangled state are not required to conform to the super-selection rule that 
prohibits quantum superpositions of single mode states with differing numbers 
of bosons. These results are based in effect on a different criterion as to what 
constitutes an entangled state, so of course the resulting inequalities will differ 
from those that would apply if the definition of an entangled state is based on 
the considerations presented here in this paper - which emphasise the require- 

s*.A ^B 

ment that the density operators p R , p R should represent physical states for the 
separate modes and hence satisfy the super-selection rule. Other results are 
based on forms of the density operator for non-entangled states that do not 
satisfy the symmetrisation principle. In this Section we examine a number of 
such entanglement tests and find that some are not valid, though some may be 
revised as tests for entangled states defined in accord with symmetrisation and 
super-selection rules. 



5.1 Hillery et al 2006 

5.1.1 Hillery Spin Variance Entanglement Test 

One such entanglement test is presented in the paper by Hillery and Zubairy 
[2"5] entitled "Entanglement conditions for two- mode states". The paper actu- 
ally dealt with EM field modes and the super-selection rule was not applied, 

^B 

so density operators p R , p R for photon modes allowed for coherences between 
states with differing photon numbers, and hence the conditions in Eq. (|75p did 
not apply. However, even for the situation of EM field modes where massless 
photons are involved, it is argued here that the super-selection rule also should 

be applied. Conditions involving the variances (AS x ) , (AS y ) can be ob- 
tained by applying similar arguements to those in Section [5l It is found that 
for the original spin operators S x , S y , S z and modes a and b 



B 
y I B 



\((Vb) R + W R ) + 1 -((^a) R (M) r ) + \((m 2 ) R ((a) 2 ) R + ((b) 2 ) R <(« f ) 2 )) 

+ m R ) + ^s^* m R ) - \i{@) a ) H (m R + (<fif) R <(« f ) 2 » 



(126) 



where terms such as ((b') 2 ) and ((a) 2 )„ previously shown to be zero have 



In 

been retained. Note that in [23] the operators S x , S y , S z constructed from the 
EM field mode operators as in Eq. ([5"5j) would be related to Stokes parameters 
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Hence 

'si) -h(.s ; 



R \ V lR 



It 



= l{@b)J(aty R + l) + ±i(tfa) R (@b) R + l)) (127) 

where the terms ((&^) 2 ) , ((a^) 2 ) cancel out. This is the same as before. 
However, 

= l^) R ®» + ^«®n> 

is now non-zero, since the previously zero terms have again been retained. Hence 
(s x ) 2 R + (S v ) 2 R = ^) R (b) R ^) R (a) R (129) 

so that we now have 

(ASl) R + (ASl) R 

= l((^b) R (@a) R + 1) + l(@a) R ((Vb) R ) + 1) - @) r (t) R (a) R @) R (a) 

= \ (( n )n + + (^)n ~ 1 (S> « n |2 ) (130) 

But from the Schwarz inequality - which is based on ((a^ — (a^))(a — (a))) > 
for any state 

\(a) R \ 2 <(aty R \(b) R \ 2 <(Pb) R (131) 

so that 

(ASl) R + (ASl) R > 1 -((n) R + (tfa) R ) (132) 
and thus from Eq. (|78p it follows that for a general non entangled state 

[ASl) + (AS 2 ) > Pr\ (fo>* + (n a ) R ) (133) 

R 

However, half the expectation value of the number operator is 

\(n) = \ ((n„ + n b )) = Pr\ ((n b ) R + (n a ) R ) (134) 

R 

so for a non-entangled state 

(aS 2 x ) + (asI)>±(n) (135) 

This inequality for non-entangled states is given in [29] (see Eq. (3)). The 
above proof was based on a different definition of entangled states to that in 
this paper. 
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5.1.2 Validity of Hillery Test for Local SSR Compatible Non-Entangled 
States 

However, it turns out that the same inequality is also valid when the definition 
of entangled states is the same as in the present paper. We would then find that 
S ^ = (^S y ^ = and hence 

ASl) + (AS^ ' ' ' 



R \ V l R 



= 2 ( v^)* + { ^«) + (( M )r ((StS) «) (136) 



instead of Eq. (|130l) . Since the term ($b\ ((a^a) R is always positive we find 
after applying Eq. (|75|) that 

AS 2 x ) + (asI)>±(n) (137) 

which is the same as in Eq. (11350 . Hence, finding that (^AS^j + ^A5^\ < 

| (^N^ would show that the state was entangled, irrespective of whether or not 
entanglement is defined in terms of non-physical unentangled states. The Hillery 
et al |29j entanglement test 

ASI) + (aSI)<~(n) (138) 

is still used in recent papers, for example [50], [71] which deal with the entan- 
glement of sub-systems each consisting of single modes a, b for a double well 
situation (in these papers S x -> J*b j & y ~~ ^ ~ ^ab > & z ~* ~ Jab ) ■ 



5.1.3 Non- Applicable Entanglement Test Involving | I^S 

Previously we had found for a general non-entangled state that is based on 
physically valid density operators p R , p R 

ASl)-\\(s z )\ > 

AS 2 y )-±\(s z )\ > (139) 
so that the sum of the variances satisfies the inequality 



(ASl) + (AS 2 v )>\(S z )\ (140) 

This is another correct inequality required for a non-entangled state as defined 
in the present paper. It follows that if only physical states /3jj, are allowed, 

the related entanglement test involving (^AS^ + (^ASy^ would be 

ASl) + (ASl)<\(s z )\ (141) 
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For any quantum state we have 

\(S*) \ =~l {{%) - (na)) | < \ dm) + (n a )) = \ (n) (142) 



which means that it is now required that (^AS 2 ^ + (^AS 2 ^ be less than a 

quantity that is smaller than in the criterion in (|135p . 

However, it should be noted that all states, entangled or otherwise, satisfy 

the inequality (a§1) + (aS^) > \ (§ Z J | so the inequality in (fT4"0|) - though 
true, is of no use in establishing whether a state is entangled in the terms of the 
meaning of entanglement in the present paper. There are no quantum states, 
entangled or otherwise that satisfy the proposed entanglement test given in Eq. 
(|141[) . This general result was stated by Hillery et al [25]. To show this we write 

the Heisenberg uncertainty principle for {^A5^ , (^AS 2 ^ as ^AS* 2 ,^) (^AS 2 ^ 
eii (S z ) | 2 , where £ > 1, then 



(ASl) + (ASl) lf ^ (AS\ 



F(y) where y 



2 V vJ i\/s 



2 



(143) 



It is straightforward to show that F(y) > 1 for all £, y. The minimum value is 
1, which occurs for £ = 1 and y — 1. Even spin squeezed states with (AS 2 ) < 



|| (S z ) | still have {AS 2 X J + \AS 2 y } > \ z ) | , so it is never found that 

(^AS 2 ^ + (^AS 2 ^ <\{^S z \\ and hence this latter inequality cannot used as a 
test for entanglement. 

Fortunately - as we have seen, showing that spin squeezing occurs via either 

(aS 2 x ^ < i| (s z \ | or (AS 2 y ^ < %\(S,)\ is sufficient to establish that the 

state is an entangled state for modes a, 6, with analogous results if principle 
spin operators are considered. Applying the Hillery et al entanglement test in 
Eq. (|138p involving \ i s & l so a valid entanglement test, but is usually less 

stringent than the spin squeezing test involving either (^AS 2 ^ < \ \ {s ■ 

AS 2 ^ < \ \ z \ |. For the Hillery et al entanglement test to be satisfied at 

least one of (^AS 2 ^ or (^AS 2 ^ is required to be less than | (^^' whereas for 

the spin squeezing test to apply at least one of (^AS 2 ^ or (^AS 2 ^ must be less 

than || (s |. The quantity || z \ | is likely to be smaller than | " 
for example the Bloch vector may lie close to the xy plane, so a greater degree 
of reduction in spin fluctuation is needed to satisfy the spin squeezing test for 
entanglement. 



59 



However, this is not always the case as the example of the relative phase state 
discussed in SubSection 14.51 shows. The results in the current SubSection can 
easily be modified to apply to new spin operators J x , J y , J z , with entanglement 
being considered for new modes c and d. The Hillery et al |29] entanglement 
test then becomes 

AJl) + (AJl)<±(N) (144) 



In the case of the relative phase eigenstate we have from Eq. (|117|) that ( A J 2 



i 

A J I) = ±N 2 + \ + \ In AT w ±N 2 for large N. This clearly exceeds \ (n^ = 
t^N, so the Hillery et al [29] test for entanglement fails. On the other hand, as we 



2 

have seen in SubSection 14.51 ( AJ^\ < || (j*) I ~ fg-^' so ^ ne s P m squeezing 
test is satisfied for this entangled state of modes c and d. 

5.2 Hillery et al 2009 

5.2.1 Hillery Strong Correlation Entanglement Test 

In a later paper entitled "Detecting entanglement with non-Hermitian opera- 
tors" Hillery et al 30 apply other inequalities for determining entanglement 
derived in the earlier paper [33] but now also to systems of massive identical 
bosons, while still retaining .density operators p R , p R that contain coherences 
between states with differing boson numbers. In particular, for a non-entangled 
state the following family of inequalities - originally derived in [29], is invoked. 

| ({a) m fo) n } | 2 < (^) m {a) m (#)"(&)") (145) 



Thus if | (Ja) m (tf) n ) I 2 > (0) m (a) m {tf) n (b) n J then the state is entangled. 

A particular case for n = m = 1 is the test \(ab^\ 2 > (n a nb) for an 
entangled state. To put this result in context, for a general quantum state and 
any operator f2 we have = and /(fit _ (ftt\\ (h - (fi))) > 0, 

hence leading to the Schwarz inequality \ (fl) \ 2 — \ (Q^ ) \ 2 < (fl^fl). Taking 



fl = atf leads to the inequality \(ab^\ 2 < (n a (rib + V)) , whilst choosing 

Q = ba) leads to the inequality | (^b^ \ 2 < ((n a + l)rib) for all quantum 
states. In both cases the right side of the inequality is greater than (n a rib), so 
if it was found that | (lab^\ \ 2 > (n a rib) (though of course still < (n a (rib + 1)) 
and < ({n a + l)rib)) then it could be concluded that the state was entangled. 
However,, as we will see the left side | (ub^ \ 2 actually works out to be zero if 

physical states for p R , p R are involved in defining non-entangled states, so that 
for a non-entangled state defined as in the present paper the true inequality 
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replacing | 1 2 ^ (^a^b) 1S j us t < (n a n~b), which is trivially true for any 

quantum state. The test for entanglement requires modification. 

The derivation of the general inequality in [29] . as in Eq. (|145[) follows 
directly from the inequality in Eq. (|23|) obtained in SubSection 12 . 51 for a general 
non-entangled state of sub-systems A and B. If we choose VIa = (a) m and f2s = 
(&)" then from | (o. A ® | 2 < (n^f2 A ® fi^s)) the result of Hillery et al 
[29] stated in Eq. (| 145|) immediately follows. The Hillery et al [29] entanglement 
test is that if 

| ((a) m (&t)»\ | 2 > /(#)"» (8)"* (#)"(&)") (146) 

then it may be concluded that the state is an entangled state for sub-systems 
A and B. Note that the proof of this result did not depend on the sub-system 
density operators p R , p R being required to satisfy SSR. 

5.2.2 Correlation Test for Entanglement for Local SSR Compatible 
Non-Entangled States 

However, for a non-entangled state based on physical p^ R ,p~ R for modes a and b 
where the SSR is satisfied we actually have 

((a)™ (&t)»\ = J2P R ((a)- @)A = EPfl <(a)"% ((#)") = (147) 

\ / R \ I R R / R 

since from Eqs. analogous to (175"]) ((a) m ) R = ((&^) n \ = 0. Hence for a physical 
non-entangled state as defined in the present paper the inequality becomes 

< ((8 f ) m (8) m (6 f r(6) n ) (148) 

which is trivially true and applies for any state, entangled or not. 

Since ((a) m (6^) n ^> is zero for non-cntangled states it follows that it is merely 

necessary to show that this quantity is non-zero to establish that the state is 
entangled. Hence an entanglement test in the case of sub-systems consisting of 
single modes a and b becomes 

\({a) m ft) n )\ 2 >0 (149) 

for a non-entangled state based on physical p R ,p R - This is a useful criterion 
for entanglement in terms the definition of entanglement in the present pa- 
per, and is different to that of Hillery et al [29]. The Hillery et al [29] entan- 
glement test | ((a)" 1 (&t)™^> | 2 > ((8t) m (8) m ft) n (b) n } is also a valid test for 
entanglement and is actually a more stringent test than merely showing that 
| ((a) m (b^) n ^ | 2 > 0, since the quantity | ((a)" 1 (^) n ^> | 2 is now required to be 
larger. In a paper by He et al [5U] (see Section IIIA) the Hillery et al [2"5] 
entanglement test | ((a) m (&t)»\ | 2 > /(8t) m (8) ro (& f )"(&)") is applied for the 
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case where A and B each consist of one mode localised in each well of a double 
well potential. This test whilst applicable could be replaced by the more easily 
satisfied test | ((a) m (6 t )™^> | 2 > 0. However, as will be seen below in SubSection 

15.61 Hillery et al [2DJ entanglement criterion is needed if the sub-systems each 
consist of pairs of modes, as treated in |72] , |50] . 



5.2.3 Examples Applying Correlation Tests for Entanglement 

As an example of applying these tests consider the mixed two mode coherent 
states described in SubSection 12.101 whose density operator for the two mode 
a, b system is given in Eq. (|4"5)l . We can now examine the Hillery et al [3DJ 
entanglement test in Eq. (1146p and the entanglement test in Eq. (|149l) for the 
case where m = n = 1. It is straight-forward to show that 

|(fi?)| a = H 4 

(j^a)$%\ = \a\ 4 (150) 

so that | (®b^\ | 2 = (^(a'a) (tfb)j. A non-entangled state defined in terms of 

the SSR requirement for the separate modes satisfies | | 2 = 0, whilst for 

a non-entangled state in which the SSR requirement for separate modes is not 
specifically required merely satisfies | (a&M | 2 < \{a^a) (b^b)). Hence the test 



for entanglement of modes A, B in the present paper \ (d,b^ \ 2 > is satisfied, 

whilst the Hillery et al [3DJ test | (ab^ | 2 > ((at a) (&t&)\ is not. 

In terms of the definition of non-entangled states in the present paper, the 
mixture of two mode coherent states given in Eq. (1481) is an entangled state, not 
a separable state. However, in terms of the definition of non-entangled states 
in other papers such as those of Hillery et al [29] , [30] the mixture of two mode 
coherent states would be a non- entangled state. It is thus a useful state for 
providing an example of the different outcomes of definitions where the local 
SSR is applied or not. 



5.3 Sorensen et al 2001 

5.3.1 Sorensen Spin Squeezing Entanglement Test 

In a paper entitled "Many-particle entanglement with Bose-Einstein conden- 
sates" Sorensen et al [H] consider the implications for spin squeezing for non- 
entangled states of the form in Eq. (|2"4")l . As discussed previously, a density 
operator of this general form is not consistent with the symmetrisation princi- 
ple - having separate density operators for specific particles i in an identical 
particle system (such as for a BEC) is not compatible with the indistinguishabil- 
ity of such particles. It is modes that are distinguishable, not identical particles, 
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so the basis for applying their results to systems of identical bosons is flawed. 
However, they derive an inequality for the spin variance ( AS ? 



A^)>^((£,) 2 + (^) 2 ) (151) 

that applies in the case of non-entangled states. Key steps in their derivation 
are stated in the Appendix to [51], but as the justification of these steps is 
not obvious for completeness the full derivation is given in Appendix [T5J of the 
present paper. This inequality (|151l) establishes that if 



1 



* = ^TT5 L,*\ < N < 152 > 
S x ) +(S 



then the state is entangled, so that if the condition for spin squeezing analogous 
to that in Eq. (|6"5[) is satisfied, then entanglement is required if spin squeezing 
for S z to occur. Spin squeezing is then a test for entanglement in terms of their 
definition of an entangled state. Note that the condition (|M|) requires the Bloch 
vector to be in the xy plane and close to the Bloch sphere of radius N/2. 



5.3.2 Revising Sorensen Spin Squeezing Entanglement Test - Lo- 
calised Modes 

The work of Sorensen et al really applies only when the individual spins are 
distinguishable. It is possible however to modify the work of Sorensen et al [UJ 
to apply to a system of identical bosons in accordance with the symmetrisation 
and super-selection rules if the index i is re-interpreted as specifying difffcrcnt 
modes, for example modes localised on optical lattice sites i = 1,2, .., N. Details 
are given in Appendix [15] With two single particle states a, b available on each 
site (these could be two different internal atomic states or two distinct spatial 
modes localised on the site) the modes would then be labelled \4> ai ) with a = a, b. 
The mode orthogonality and completeness relations would then be 

(<P a i\<f>0j) = S a pSij 

ElOfaail = 1 (153) 

ai 

With the particles now labelled K = 1,2,3, ...one can define spin operators in 
first quantization via 

K i 

Sy = J2J2^im{^K)\-\4> ai (K)){cf> bi (K)\)/2i 

K i 

S> = J2J2^i(K)){him-\cl> ai (K))^ ai (K)\)/2 (154) 

K i 
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In second quantization if the annihilation, creation operators for the modes \4> ai ) 
,\4> bi ) are a,i, bi and aj, b\ respectively, then the Schwinger spin operators are 
just 

s x = y.^ + ^)I 2 = H^ 

i i 
i i 

S z = Y,^-a%)/2 = Y^Sl (155) 

i i 

It is easy to confirm that the overall spin operators S a and the spin operators 
for the separate pairs of modes \4> ai ) ,\4>u) ( or f° r short) satisfy the 

same commutation rules as Sorensen et al |41j have for the overall spin operators 
and those for the separate particles. With this modification the non-entangled 
state in Eq. (pM|) could be interpreted as being a non-entangled state where the 
subsystems are actually pairs of modes \4> a i) ^u) an d the density operators 
~p % R would then refer to a subsystem consisting of these pairs of modes. It is to 
be noted that entanglement of pairs of modes is different to entanglement of 
all separate modes. It is an example of a special kind of multimode entangle- 
ment - since the modes \4> ai ) i\4>m) m ay themselves be entangled we may have 
"entanglement of entanglement". In terms of the present paper the density op- 
erators /S^ would be restricted by the super-selection rule to statistical mixtures 
of states with specific total numbers Ni of bosons in the pair of modes \4> a i) 
In terms of Fock states \n a i) , |«.6i) for this pair of modes the allowed 
quantum states for the sub-system will be 

\*n<)= E<* \k) ai \Ni-k) bi (156) 

fc=0 

so at this stage the general mixed physical state for the two mode system could 
be 

oo Ni Ni 

PH= E £JW, E E A k* ( A ?T \ k )ai (Aai ® \ N i ~ k )bi ( N i ~ l \bi 
Ni—Q <S> k=0 1=0 

This state has no coherences between states of the two mode subsystem with 
differing total boson number 2Vj for the pair of modes. However this is still an 
entangled states for the two modes \</> ai ) ,1^), so the overall state in Eq. (|24l) 
is not a non-entangled state if the subsystems were to consist of all the distinct 
modes. 

5.3.3 Revising Sorensen Spin Squeezing Entanglement Test - Sepa- 
rable State of All Modes 

It is possible however to link spin squeezing and entanglement in the case where 
the sub-systems consist of all the distinct modes. To obtain a fully non-entangled 
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state of all the modes \4> ai ) ,\4>bi) the density operator p l R must then be a product 
of density operators for modes \<j) ai ) and 

Pr = Pr®Pr (158) 
giving the full density operator as 

p = E p « ® ?*) ® (^ 2 ® ) ® (% 3 ® ® ■ ( 159 ) 

R 

as is required for a general non-entangled state all 2N modes. Furthermore, 
as previously the density operators for the individual modes must represent 
possible physical states for such modes, so the super-selection rule for atom 
number will apply and we have 

((SiDai = Tr(p a n(a i r) = ($)») = = 

((Sir) = T>o3^n = o ((&ir\ = Tr(pS*®r) = o 

\ i bi \ i b% 

(160) 

The question is whether this reformulation will lead to a useful inequality 
for the spin variances such as (AiS^). This issue is dealt with in Appendix 
1151 and it is found that we can indeed show for the general fully non- entangled 
state (fT59|l that 

(ASI)>±\(S Z )\ and (A£»)>i|(£,)| (161) 

This shows that if there is spin squeezing in either S x or S y then the state must 
be entangled. Note that this result depends on the general non-entangled state 
being non-entangled for all modes and that the density operator for each mode 
Oj or bi being a physical state with no coherences between mode Fock states with 
differing atom numbers. In terms of the revised interpretation of the density 
operator to refer to a multi-mode system with modes \4> ai ) , the statement 
that spin squeezing for systems of identical massive bosons requires all the modes 
to be entangled is correct. However superposition states of the form (|156|) that 
are consistent with the super-selection rule applying to pure states of a two mode 
system are precluded, and such states ought to be allowed if entanglement of 
pairs of modes rather than of separate modes is to be considered. 

5.3.4 Revising Sorensen Spin Squeezing Entanglement Test - Sepa- 
rable State of Pairs of Modes 

It is also possible however to link spin squeezing and entanglement in the case 
where the subsystems consist of pairs of modes, but only if further restrictions 
are applied. The general non-entangled state of the pairs of modes would actu- 
ally be of the form 

R 
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where the ~p l R are now of the form given in Eq. (I157[) and no longer are density 
operators for the ith identical particle. Unlike in (| 160[) we now have expectation 
values ((fli)™) i — Tr(p R (a~i) n ) etc that are non-zero, so considerations of the 
link between spin squeezing and entanglement - now entanglement of pairs of 
modes, willl be different. 

If the density operators ~p R associated with the pair of modes ai, hi are all 
restricted to be associated with one boson states then this density operator is 
of the form 

? R = pL(|i}i <iL®|o) i6 <o| i6 ) + p l afe (|i} ia <o| ia ®|o} i6 <i| i6 ) 
+4a(|o) io UL ® |i> a (o| tt ) + Pi b (\o) ia (o\ ia ® <iy 

(163) 

where the p\^ are density matrix elements. With this restriction the pair of 

modes tij, hi behave like distinguishable two state particles, essentially the case 
that Sorensen et al [JT] implicitly considered. The expectation values for the 
spin operators S l x , S y and S\ associated with the ith pair of modes are then 

Sl) R = \{pl b - P l aa ) (164) 

If in addition Hermitiancy, positivity, unit trace Tr(p R ) — 1 and Tr(p R ) 2 < 
1 are used (see Appendix [13} then we can show that p % hh and p l aa are real and 
positive, pl b = (plJ* and p l aa p\ h - \p l ab \ 2 > 0. The condition Tr(p R ) = 1 leads 
to p l aa + p\ h = 1, from which Tr(p R ) 2 < 1 follows using the previous positivity 
results. These results enable the matrix elements in (|163l) to be parameterised 
in the form 

Paa = sin2 «i Pbb = CC, S 2 "i 

Pat, = \J sin 2 aj cos 2 a, sin 2 exp(+i0 i ) = ^/ sin 2 a, cos 2 a, sin 2 /3j exp(— 

(165) 

where cnj, /3j and 4> i are real. In terms of these quantities we then have 

Si) = — sin 2cti sin 2 B- cos </>- /s?,^ = - sin 2a^ sin 2 5,. sin </>■ 
Ir 2 \ v /r 2 

S*) = i cos 2a, (166) 

'Ft A 

and then a key inequality 

(st) 2 R + (SI) 2 R + (SI) 2 R = \-\ sin 2 2a t (1 - sin 4 & ) < \ (167) 
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follows. This result depends on the density operators p l R being for one boson 
states, as in (jl63|) . The same steps as in Sorensen et al [41] (see Appendix IT4l) 
leads to the result 



AS: 



> 



1 

N 



S x 



S, 



(168) 



for non-entangled pair of modes a,, hi. Thus when the interpretation is changed 
so that are the separate sub-systems are these pairs of modes and the sub- 
systems are in one boson states, it follows that spin squeezing requires entan- 
glement of all the mode pairs. 

A similar proof extending the test of Sorensen et al [41] to appply to systems 
of identical bosons is given by Hyllus et al j44 based on a particle entanglement 
approach. In their approach bosons in differing external modes (analogous to 
differing i here) are treated as distinguishable, and the symmetrization principle 
is ignored for such bosons. 



5.4 Sorensen and Molmer 2001 

In a paper entitled "Entanglement and Extreme Spin Squeezing" Sorensen and 
Molmer |73] first consider the limits imposed by the Heisenberg uncertainty 

principle on the variance (AJ|) considered as a function of | ( J z ) | for states 



where the spin operators are chosen such that (J x ) = \ J y) — 0- Note that 
such spin operators can always be chosen so that the Bloch vector does lie along 
the z axis, even if the spin operators are not principal spin operators. Their 
treatment is based on combining the result from the Schwarz inequality 

2 



J- 



J: 



J, 



< J(J + l) 



where J = N/2, and the Heisenberg uncertainty principle 



J: 



(169) 



(170) 



where £ > 1. In fact two inequalities can be obtained 
AJ 2< 



AJ; 



> 



< - 



J(J + 1) 




J . 



J(J + 1) 



J : 



(171) 



1) — (J 




J(J+l) 



J 



(172) 



which restricts the region in a ^AJ^y versus | (j z j \ plane that applies for 
states that are consistent with the Heisenberg uncertainty principle. The first 
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of these two inequalities is given as Eq. (3) in [73] ■ For states in which J x is 
squeezed relative to J y the points in the/AJ^\ versus | (j z) I plane must also 
satisfy 

(AJ») < \\ (7 Z ) | (173) 

Note that as J z is a spin angular momentum component we always have 
J z) \ < J, which places an overall restriction on I (J z ) I. However, for 



£ > 1 there are values of | ( J z ) | .which are excluded via the Heisenberg un- 



2 ■ ■ a 



certainty principle, since the quantity ^ J( J + 1) — \J z/ j ~ £ z ) then 
becomes negative. This effect is seen in Figure 4. 

The question is: Is it possible to find values for ^AJ 2 ^} and | (jz^ I in 
which all three inequalities are satisfied? The answer is yes. Results showing 
the regions in the ^AJ 2 ^} versus |^J Z ^| plane corresponding to the three 
inequalities are shown in Figures 2 and 3 for the cases where J = 1000 and with 
£ = 1.0 and £ = 10.0 respectively. The quantities for which the regions are shown 
are the scaled variance and mean ^AJ 2 ^} /J and | (jz\ \/J, with ^AJ 2 ^} given 

as a function of | z^J | via (|171[) . (| 1 72[) and (jlT3|) . The spin squeezing region 
is always consistent with the second Heisenberg inequality (|172[) and for large 
J = 1000 there is a large region of overlap with the first inequality (|17ip . For 
small J and large £ the region of overlap becomes much smaller, as the result in 
Figure 4 for J = 1 and with £ = 10.0 shows. As the derivation of the Heisenberg 
principle inequalities is not obvious, this is set out in Appendix llGI 

Sorensen and Molmer [73] also determine the minimum for (AJ 2 ) = (J 2 



as a function of | \ J zj \ for various choices of J, subject to the constraints 

J x^J = {^J y ^ = 0. The results show again that there is a region in the 

AJ 2 ^ versus | z ^ | plane which is compatible with spin squeezing. 
So although these considerations show that the Heisenberg uncertainty prin- 
ciple does not rule out spin squeezing, nothing is determined about whether the 
spin squeezed states are entangled states for modes c, d , where the J a are 
given as in Eq. ([M]) . The discussion in [73] regarding entanglement is based on 
the physically incorrect density operator for non-entangled states of identical 
particles in Eq. ([24]) , discussed in the previous section. 



5.5 Duan et al 2000 

A further inequality aimed at providing a signature for entanglement is set out 
in the papers by Duan et al [73], Toth et al [75]. For simplicity we only set out 
the case for which a = 1 in the former paper. This inequality involves position 
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and momentum like Hermitian operators denned by 



x A = -^=(a + a f ) p A = -^-.{a-a*) 



XB 



-i=(6 + 6+) p B = -L.(b-ft) (174) 



These are essentially quadrature operators and satisfy commutation rules [xa, pa] 
\xb,Pb] — i similar to those for position and momentum. An inequality is ob- 
tained for a general two mode non-entangled state involving the variances for 
the commuting observables xa + xb and pa — Pb 

{A(x a +x b ) 2 ) + {A(pa-pb) 2 )>2 (175) 

which could be used to establish a quadrature variance test for entangled states 
of the mode A and mode B sub-systems, so that if 

{A{x a +x b ) 2 ) + {A(pa-pb) 2 ) < 2 (176) 

then the modes are entangled. Such states are possible - consider for example 
any simultaneous eigenstate of the commuting observables xa + x B and pa — 
Pb- For such a state (A(xa + xb) 2 ) and (A(j>a — Pb) 2 ) are both zero, so the 
simultaneous eigenstates are entangled states of modes A, B. 

To confirm whether the inequality (|175l) applies for non-entangled states 
(jlOip in which the sub-system states p^, Pr are physical, the general variance 
result in Eq. ((75)) plus the factorisations (xaxb) r = (xa) r (xb)r and (j>aPb) r 
= (pa) r (pb) r are first used to show that 

(A(x A + x B ) 2 ) + (A(p A - p B f) 

- Pr ((^ 2 ) R - @*)r + (*B 2 ) r - (xb)r + (PA 2 ) r - (pa)r + (PB 2 ) R 



For sub-system states p^, p^ that are physical we have in addition (xa 



R 

2\ 



(xb) r = (pa) r = (pb) r = 0. Also using (a 2 ) R = (@) 2 ) R = (b : 

(^(tf) 2 ^ =0 for physical states we find for the remaining terms in Eq. (|177p 
that 

= i({a J > li + ((s t ) 2 >„ + i + 2(ats> ii )>i 
(pa 2 ) r = -i((a 2 > JJ + ((at ) =) ;i -i-2{at !i)ji )>i 
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Substituting these results into Eq. (|177p establishes the validity of (|175[) for non- 
entangled states in which the p R , p R are physical sub-system states. As there 
are entangled states that violate the inequality (| 1 T5[) . this inequality is valid for 
determining whether a state is entangled. 

5.6 He et al 2012 

In two papers dealing with EPR entanglement He et al [75], [S0J a four mode 
system associated with a double well potential is considered. In the left well A 
there are two localised modes with annihilation operators Si , a,2 and in the right 
well B there are two localised modes with annihilation operators b\ 1 &2- The 
modes in each well could be associated with different internal states or they could 
be associated with different spatial modes of the same internal state. This four 
mode system provides for the possibility of entanglement of two sub-systems 
each consisting of pairs of modes. With four modes there are three different 
choices of such sub-systems but perhaps the most interesting from the point of 
view of entanglement of spatially separated modes - and hence implications for 
EPR entanglement - would be to have the two left well modes 01, S2 as sub- 
system A and the two right well modes 61, 62 as sub-system B. Consistent with 
the requirement that the sub-system density operators p R , p R conform to the 
symmetrisation principle and the super-selection rule, these density operators 
may now be of the form given in Eq. (|47j) . Hence as discussed in Sub-System 
12.121 when considering non- entangled states for the sub-systems A and B we 
now have as in Eq. ([54)) 

A = Tr(p£(a t r)^0 {(al) n ) A = Tr(pi(alr)^0 
((M m ) B = Tr(p B R %r)^Q ((b]) m ) B =Tr(p^(b]r)^0 

(179) 

in general. Hence in this case where the sub-systems are pairs of modes the 
entanglement test in Eq. (| 149|) for sub-systems consisting of single modes cannot 
be applied. 

5.6.1 Correlation Tests for Entanglement 

However, the inequalities derived by Hillery et al [30] (see SubSection l5.2j) 

1 ((«,)'" (b]r) 1 2 < ((flirt*)™ tyn%r) (m 

that apply for two non-entangled sub-systems A and B can now be usefully 

applied, since in this case the quantities (^(ai) m are in general no longer 

zero. Thus there is an entanglement test for two sub-systems consisting of pairs 
of modes. If 

\(&r(b]r)\ 2 > ((alr&r^nbjr) 

foranyofij = 1,2 (181) 
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then the quantum state for two sub-systems A and B - each consisting of two 
modes localised in each well - is entangled. 

5.6.2 Spin Squeezing Tests for Entanglement 

There are numerous choices for defining spin operators but the most useful 
would be the local spin operators for each well |50) defined by 

S% = (a|oi + a[a 2 )/2 S$ = (a|oi - a\a 2 )/2i = {a\a 2 - a\a x )/2 

= (b%+b%)/2 S B = (b 2 b 1 -V 1 b 2 )/2i S B = (? 2 6 2 -^ip2) 

These satisfy the usual angular momentum commutation rules and those or the 
different wells commute. The squares of the local vector spin operators are 
related to the total number operators Na = a\a 2 + a\ai and N B = b\b 2 + b\bi 
as E(S*) 2 = ( N A /2)(N A /2 + 1) and ) 2 = ( N B /2)(N B /2 + 1). 

a. a. 

For the local spin operators we have in general 

(S£) A = TrC P A R Si) £ (S*) g = TtCp B r Sf ) * (183) 

since the pair of modes Oi, a 2 and/or b\, b 2 may now be of the form given in 
Eq. dUD- 

In SubSection 12.51 it was shown that | (ti^Sls^ | 2 < ^Q^Q^Q^ for a 

non-entangled state, so with Qa = £_ = S£-iSy and Q B = S B = S B ~iS B = 
(SfY to give 

\(s^S B )\ 2 < (s^S A S B S B ) (184) 

for a non-entangled state of sub-systems A and B. For the non-entangled state 
of these two sub-systems we have 

(^S B )=^Pn(si) R m R B (185) 

R 

which in general is non-zero from Eq. (|183p . 

Hence a valid entanglement test involving spin operators for sub-systems A 
and B - each consisting of two modes localised in each well exists, and is if 

\(s A S B )\ 2 > (s A S A S B S B ) (186) 

then the two sub-systems are entangled. A similar conclusion is stated in |50j . 
This test for entanglement involves the local spin operators, though it is not 
then the same as spin squeezing criteria. It is referred to as spin entanglement. 
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6 Experiments on Spin Squeezing 



There are several papers [76] , [77] , [78] which contain the results of measuring 
the spin squeezing parameter analogous to the expression in Eq. (1651) and show- 
ing that spin squeezing occurs. The presence of entanglement is then inferred 
by reference to theoretical papers such as [H] that show that spin squeezing 
only occurs for entangled states - it is an entanglement witness. As no indepen- 
dent measures of entanglement (however defined) are presented, nor are other 
independent tests for entanglement carried out, it cannot be said that these 
paper shows experimentally that spin squeezing requires entanglement. In [78j 
the emphasis is on showing how the spin squeezing can be generated via the 
non-linear terms in the Josephson Hamiltonian. 
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7 Discussion and Summary of Key Results 



This paper is mainly concerned with two mode entanglement for systems of 
identical bosons and the relationship to spin squeezing. These bosons may be 
atoms or molecules as in cold quantum gases or may be photons as in quan- 
tum optics. The analysis applies to both though detailed arguements are often 
framed for the bosonic atom case. 

A careful analysis is first given regarding the proper definition of a non- 
entangled state for systems of identical particles, and hence by implication the 
proper definition of an entangled state. Noting that entanglement is meaningless 
until the subsystems being entangled are specified, it is pointed out that whereas 
it is not possible to distinguish identical particles and hence the individual 
particles are not legitimate sub-systems, the same is not the case for the single 
particle states or modes, so the modes are then the the rightful sub-systems to be 
considered as being entangled or not. In this approach where the sub-systems 
are modes, situations where there are differing numbers of identical particles 
are treated as different physical states, not as differing physical systems, and 
the symmetrisation principle required of physical states for identical particle 
systems will be satisfied by using Fock states to describe the states. 

Furthermore, it is argued that the overall physical states should conform to 
the superselection rule that excludes quantum superposition states of the form 
(|26|) as physical states for systems of identical particles - massive or otherwise. 
Although the justication of the SSR in terms of observers and their reference 
frames formulated by other authors has also been presented for completeness, 
a number of fairly straightforward reasons were given for why it is appropriate 
to apply this superselection rule, which may be summarised as: 1. No way 
is known for creating such states. 2 No way is known for measuring all the 
properties of such states, even if they existed. 3. There is no need to invoke 
the existence of such states in order to understand coherence and interference 
effects. 4. The stability of such states against decoherence processes may not be 
great, so even if they could be created they could rapidly change to other states. 
The last reason is of lesser importance. Invoking the physical existence of states 
that as far as we know cannot be made or measured, and for which there are no 
known physical effects that require their presence seems a rather unnecessary 
feature to add to the non-relativistic quantum physics of many body systems, 
and considerations based on the general principle of simplicity (Occam's razor) 
would suggest not doing this until a clear physical justification for including 
them is found. As two mode fermionic systems are restricted to states with at 
most two fermions, the focus of the paper is then on bosonic systems, where 
large numbers of bosons can occupy two mode systems. 

However, although there is related work involving local particle number 
super-selection rules, this paper differs from a number of others by extend- 
ing the super-selection rule to also apply to the density operators /3^, p^, .. 
for the mode sub-systems A, B, .that occur in the definition (TTJ of a general 
non-entangled state for systems of identical particles. Hence it follows that the 
definition of entangled states will differ in this paper from that which would 
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apply if density operators p R , p R , .allowed for coherent superpositions of num- 
ber states within each mode. In fact more states are regarded as entangled in 
terms of the definition in the present paper. Indeed, if further restrictions are 
placed on the sub-system density operators - such as requiring them to spec- 
ify a fixed number of bosons - the set of entangled states is further enlarged. 
The simple justification for our viewpoint on applying the local particle number 
super-selection rule has three aspects. Firstly, since experimental arrangements 
in which only one bosonic mode is involved can be created, the same reasons 
(see last paragraph) justify applying the super-selection rule to this mode sys- 
tem as applied for the system as a whole. Secondly, measurements can be 
carried out on the separate modes, and the joint probability for the outcomes 
of these measurements determined. For a non-entangled state the joint proba- 
bility ([4]) for these measurements depends on all the density operators f? R , 'p R , 
.. for the mode sub-systems as well as the probability Pr for the product state 
?fl ® Pit ® -.occuring when the general mixed non-entangled state is prepared, 
which can be accomplished by local preparations and classical communication. 
For the non-entangled state the form of the joint probability fU.B..(*>ii ••) f° r 
measurements on all the sub-systems is given by the products of the individ- 
ual sub-system probabilities P^{i) = Tr(Uf'p R ), etc that measurements on 
the sub-systems A, B, ..yield the outcomes \f etc when the sub-systems are in 
states p~ Rl ~p R , ., the overall products being weighted by the probability Pr that 
a particular product state is prepared. If "p R , p~ R , did not represent physical 
states then the interpretation of the joint probability as this statistical average 
would be unphysical Thirdly, attempts to allow the density operators p R , p R , 
. . for the mode sub-systems to violate the super-selection rule provided that the 
reduced density operators ~p A , ~p B for the separate modes are consistent with it 
are shown not to be possible in general. 

As well as the above justifications for applying the super-selection rule to 
both the overall multi-mode state for systems of identical particles and the 
separate sub-system states in the definition of non-entangled states, a more 
sophisticated justification based on considering SSR to be the consequence of 
describing the quantum state by an observer (Charlie) whose phase reference is 
unknown has also been presented in detail in Appendix Qj] for completeness. For 
the sub-systems local reference frames are involved. The SSR is seen as a special 
case of a general SSR which forbids quantum states from exhibiting coherences 
between states associated with irreducible representations of the transformation 
group that relates reference frames, and which may be the symmetry group for 
the system. 

The present paper then defines spin squeezing for two mode systems and 
discusses the desirability of defining spin squeezing in terms of the principal 
spin operators J x , J y , J z for which the covariance matrix is diagonal, rather than 
via the original spin operators S x , S y , S z defined in terms of the original mode 
annihilation and creation operators a, b and , W and for which the covariance 
matrix is non-diagonal in general. It is seen that the two sets of spin operators 
are related via a rotation operator and the principal spin operators are given 
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in terms of new mode operators c, d and , Oft, with c, c? obtained as linear 
combinations of the original mode operators a, b and hence defining two new 
modes. 

The immediate consequence for the case of two mode systems of identical 
bosons of the present approach to defining entangled states is that spin squeezing 
in any of the principle spin operators J x , J y or J z requires entanglement of the 
new modes c, d. Similarly, spin squeezing in any of the original spin operators 
S x , S y or S z requires entanglement of the original modes a, b. A typical test 

for entanglement is (aS^ < | (S z \ |/2 or (AS^ < | (s z ^ |/2. It is noted 
that though spin squeezing requires entanglement, the opposite is not the case 
and the NOON state provides an example of an entangled physical state that 
is not spin squeezed. Also, the binomial state provides an example of a state 
that is entangled and spin squeezed for one choice of mode sub-systems may 
be non-entangled and not spin squeezed for another choice. The relative phase 
state provides an example that is entangled for new modes c, d and is highly 
spin squeezed in J y and very unsqueezed in J x . The connection between spin 
squeezing and entanglement is regarded as well-known, but up to now only 
proofs based on non-entangled states that either disregard the symmetrization 
principle or the sub-system super-selection rules exist, placing the connection 
between spin squeezing and entanglement on a somewhat shaky basis. On the 
other hand, the proof given here is based on a definition of non-entangled (and 
hence entangled) states that is compatible with both these requirements. 

There are several papers that obtain different tests for whether a state is 
entangled from those involving spin squeezing that are obtained in this pa- 
per, the proofs often being based on a definition of non-entangled states that 
ignores symmetrization or SSR. Hillery et al [55] obtain criteria of this type, 

such as the entanglement test {AS^j + IaS^) < | (^/- This test is also 
valid if the non-entangled state definition is consistent with the SSR, but is 
different to the test /AS^ + (aS^) < \ (s z ^ | suggested by the require- 
ment that (^AS'lj + (^ASy\ > | (Sz) I f° r non-entangled states - since both 
(ASl") > |(^)|/2 and /a5^\ > \(s z }\/2. The latter inequality is of 
no use since (^AS^ + (^ASy^ > \ (^S z ^ \ for all states. However as previ- 
ously noted, showing that cither (aS^ < | (s z \ |/2 or (AS^ < | (s z ^ |/2 
- or the analogous tests for other pairs of spin operators - already provides a 
test for the entanglement of the original modes a, b. This test is a different 
test for entanglement than that of Hillery et al [29]. The case of the rela- 
tive phase eigenstate is an example of an entangled state in which the spin 
squeezing test for entanglement succeeds whereas that of Hillery et al [33] fails. 
Other inequalities found by Hillery et al [30 for non-entangled states which 
also do not depend on whether non-entangled states satisfy the super-selection 
rule include | ({a) m (6*)"^ | 2 < ({d)) m {a) m (6 t )"(6) n ^), giving another valid test 
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I {{a) m (tf) n J | 2 > ^(at) m (a) m (tf) n (b) n ) for an entangled state. However, with 

entanglement defined as in the present paper we have | ^(a) m (^) n y | 2 = for 

a non-entangled state, so an entanglement test in the form | (^(a) m (6^)"^ | 2 > 
immediately follows. This test is less stringent than that of Hillery et al [3D]., as 
| ^(a) m (&t) n \ | 2 is then required to be larger. Sorensen et al 41, show that spin 
squeezing is a test for a state being entangled, but define non-entangled states for 
identical particle systems (such as BECs) in a form that is inconsistent with the 
symmetrisation principle - the sub-systems being regarded as individual iden- 
tical particles. However, the treatment of Sorensen et al [41] can be modified 
to apply to a system of identical bosons if the particle index i is re- interpreted 
as specifying diffferent modes, for example modes localised on optical lattice 
sites i = 1,2, ,.,N. With two single particle states \(j) ai ) ^m) wrtri annihilation 
operators <Zj, 6, available on each site, there would then be 2N modes involved, 
but spin operators can still be defined. If the definitions of non-entangled and 
entangled states in the present paper are applied, it can be shown that spin 
squeezing in either of the spin operators S x or S y requires entanglement of all 

the original modes at,bi. Alternatively, if the sub-systems are pairs of modes 
a>i,bi and the sub-system density operators /3^. were restricted to states with 
exactly one boson, then it can be shown that spin squeezing in S z requires en- 
tanglement of all the pairs of modes. With this restriction the pair of modes a*, 
bi behave like distinguishable two state particles, which was essentially the case 
that Sorensen et al (41 j implicitly considered. This type of entanglement is a 
multi-mode entanglement of a special type - since the modes aj, bi may them- 
selves be entangled there is an "entanglement of entanglement". So with either 
of these key revisions, the work of Sorensen et al [41] could be said to show that 
spin squeezing requires entanglement. Sorensen and Molmer |73j have deduced 

an inequality involving ^AJ 2 ^} and | ^ J | for states where (j x \ = J = 

based on just the Heisenberg uncertainty principle. This is useful in terms of 
confirming that states do exist that are spin squeezed still conform to this prin- 
ciple. Duan et al [71], Toth et al [75] devise a test for entanglement based on the 
sum of the quadrature variances (A(xa + x b) 2 ) + (A(p : a — Pb) 2 ), which involve 
quadrature components xa,Pa,xb,Pb constructed from the original mode an- 
nihilation, creation operators for modes A, B. Their conclusion that if the sum 
is less than 2 then the state is entangled is valid both for the present definition 
of entanglement and for that in which the application of the super-selection 
rule is ignored. He et al [72], [50] consider a four mode system, with two modes 
localised in each well of a double well potential. If the two sub-systems A and B 
each consist of two modes - with a\ , a>2 as sub-system A and b\ , 62 as sub-system 
B, then tests of entanglement of the two sub-systems of the Hillery [30] type 

I ((a*)™ (£])") I 2 > ((al) m (ai) m (^) n (6j') n ) for any i, j = 1, 2 or involving 

local spin operators | (sf 5 s ^ | 2 > (s$ S A Sf S B apply. 

Overall then, all of the entanglement tests (spin squeezing and other) in 
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the other papers discussed here are still valid when reconsidered in accord with 
the definition of entanglement based on the symmetrisation and super-selection 
rules, though in one case Sorensen et al [H] a re-definition of the sub-systems 
is required to satisfy the symmetrization principle. However, further tests for 
entanglement are obtained in the present paper based on non-entangled states 
that are consistent with the symmetrizaton and super-selection rules. In some 
cases they are less stringent - the correlation test in Eq. (1149[) being easier to 
satisfy than that of Hillery et al [3D] in Eq. (|146p . They are certainly different 
to others previously discovered. 

At present, experiments demonstrating spin squeezing do not show experi- 
mentally whether spin squeezing requires entanglement, however defined, since 
no results for entanglement measures are presented, nor are other independent 
tests for entanglement carried out. 
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9 Appendix 1 - Inequalities 
9.1 Integral Inequality 

If C(A) is a real, positive function of A and P(A) is another real, positive function 
such that J dXP(X) = 1 then we can show that 

fdXP(X)C(X) > (/ dXP(X)^C(X)\ (187) 
To show this write y = J dXP(X)C(X). Then 

y = Jd^p(fi)jdXP(x)c(X) 

= JdXP(X) 2 C{X)+ffdXdn(l~5(X-n))P{fi)P(X)C(X) (188) 

Also, write z = (j dX F(A)^C"(A)) 2 . Then 

z = fd»P(fi)^/C(ri JdXP(X)VC(Xj 

= J fdXd[iP(n)P(X)^/C(riVC(X) 

= JdX P(X) 2 C(X) +JJdXdfi(l- 5(X - n)) P(n)P{X) y/C(fi)y/C(l$9) 
so that 



y-z = J JdX dp (1 - 6(X - n)) P(p)P(X) (C(X) - VCi^^C{X) 

= I J J dX d/z (1 - S(X - /i)) P(/x)P(A) (C(X) + C(ji) - 2y^)Vc(A) 

= \ J J dx d M (i - 6(x - n)) p(m)p(a) (Vc(xj - JW)) 2 

> (190) 
which proves the result. 

9.2 Sum Inequality 

If Cr is a real, positive quantity for various R and Pr is another real, positive 

quantity such that E Pr — 1 then we can show that 
R 

E Pr c R > ( E p* (i9i) 

This inequality is used in [29]. To prove this write y = E Pr C-R- Then 

ft 

i/ = EPsEPfl^ 

S ft 
= EE p sPr Cr 

R S 

= EPrCr + EE(^-Srs)PsPrCr (192) 

R R S 
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Also, write z = ( E P R VCr ) . Then 
R 



E Ps VCsj (E Pr VCr 



= EE PsPrVCsVCr 

R S 

= £ c« + £ £ (i - s Rs) Ps Pr Vc^Vcr 

R R S 



so that 



y - z 



= EE PsPrO-- Srs) (c r - Vc^Vc^) 

R S V ' 

= JEE PsPr(1-6 rs ) (Cr + C s -2^/C~s\ 

Z R S V 
= ^EE P S Pr(1~5r S ) (y/C^-y/Cs 



> o 

which proves the result. 
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10 Appendix 2 - Particle and Mode Entangle- 
ment 



It is useful to contrast the two meanings of entanglement - mode and particle 

- in terms of three examples. The first is from the textbook by Peres ([5], see 
ppl26-128). A system with N = 2 identical particles has one particle in a single 
particle state (mode) \u), the other in an orthogonal single particle state \v). In 
first quantization the synmmetrized quantum pure states for identical bosons 
or for identical fermions are (my notation) are 

l*) 6oson = ^(Ki))®K2)) + K2))<»Ki))) 

m fermion = ^(K 1 )) ® K 2 )> " K 2 )> ® K 1 ))) ( 195 ) 

which consequently means that " two particles of the same type are always entan- 
gled" . Peres obviously considers such entanglement is a result of symmetrization. 
In second quantization the state in both the fermion and boson cases is |1) ® 
jl)^ which is a separable state for modes u, v, and not a (mode) entangled state. 

The second example is taken from the paper of Hyllus et al [44] , specifically 
a case illustrated in Fig 1(b) which shows a state with N = 5 identical bosons. 
The bosons may occupy differing spatial states (eg harmonic oscillator states) 

- referred to by Hyllus et al as external degrees of freedom - and each bosonic 
particle has two distinct internal states (eg hyperfine states) - internal degrees 
of freedom. Fig 1(b) shows two spatial states and two internal states (it, d 
say) , with only the lower spatial state (0 O say) being occupied by N = 5 
bosons. From the Hyllus et al viewpoint (see last para on p 012337-4) "For 
indistinguishable particles, only two possibilities are allowed in this case: either 
ALL the particles are in a separable state (that is, product \4>)® N ) state, or all 
particles are entangled due to the symmetrization." Hyllus et al describe the 
states in terms of first quantization but for purposes of comparison we will also 
describe them via second quantization. What they mean by the separable state 
is in full 

\<f>f N = \<h) I02> I03> I04> I05> (196) 

where for the ith particle the single particle space-spin state would of the form 

\(t>i) = (cosfl \ Ui ) + sin 6 expix |d»>) ® \(j> Qi ) (197) 

in which a particular internal state is chosen The separable state in Eq. (|196[) is 
just a tensor product of single particle states for the five bosons. It is symmetric, 
so the symmetrization principle is satisfied. There is of course one other orthog- 
onal separable state = |^) |£ 2 ) IC3) (£4) with an orthogonal single 
particle space-spin state |^) = (— sin# |ttj) + cos# expzx \di)) ® \cj> 0i ) in which 
the internal state is orthogonal to the previous one. If one of the bosons is taken 
from a state \<j>) and placed in the orthogonal state then representing it in 
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the form of a single tensor product such as j^-J \<f> 2 ) \4>3) 104) |£s) would not sat- 
isfy the symmetrization principle. If one such product as \<j>-y) \<f> 2 ) \<t>%) ^4) IO 
is subjected to an operator which is the sum of all permutation operators P , 
then apart from normalising factor the result will represent the situation where 
one of the five bosons is in the state |£) rather than \<j>). Hence such a state is 
given by 

|*4,l) = MY,P (I0l)l0 2 ) |0 3 ) |0 4 ) l£ 5 )) 
P 

kf# ( I0l) I02> I0 3 ) I04> IO + I'M \<t>2) l^3> I0 5 ) l&> + I0l) I0 2 ) |0 5 ) |0 4 ) \U) 
I + I0l) I0 5 ) I0 3 ) I04> l£ 2 ) + I05> I02> !03> I04> Kl) 

ir# f I0l) 102) 103) 104) |&> + I0l) 102) 103) 1^4) 105) + I0l) 102) |&> I04> I0p)n 
" I +I0l) 1^2)103)104)105) + l^l) 102) 103) 104) 105) 

which are where the sum is over the 5! permuation operators and the Af's are 
normalising factors. However, Hyllus et al refer to this as entanglement by 
symmetrization and regard this state as being entangled. From this point of 
view it is symmetrization via P that is responsible for entanglement in that 

p 

it creates contributions to the state vector which becomes no longer just a simple 
product. There is a term 1^) |0 2 ) \<j> 3 ) |0 4 ) |£ 5 ) followed by |<^) \<f> 2 ) |0 3 ) |£ 4 ) \4> 5 ) 
in which particles 4 and 5 are in different single particle states. 

However, from the opposing point of view in which it is modes, not par- 
ticles that are entangled, and the state just described would not be regarded 
as being entangled. The Fig 1(b) case would be seen as a two mode situation 
in which the two modes are \U) — \u) <g> \(j) ) and \D) = \d) ® |(/> ). In second 
quantization the Fock states are \nu ,Ud) = \nu) <8> \no) with n\j ,nn being 
the mode occupancies. It is these two modes that may or may not be entan- 
gled, and there are six separable pure states (not two) with a total of N = 5 
bosons, namely |5,0), |4,1), |3,2), |2,3), |1,4), and |0,5). The states |5,0) 
and |0 , 5) are of course equivalent in first quantization to \<j>^) \<j) 2 ) 103) 104) |0s) 
and |^) |£ 2 ) ^3) IC4) ^5)) whilst the state in the last equation is just the sepa- 
rable state |4, 1). The general mode entangled pure state with N = 5 bosons 
is given by 

I*) = 2? 5 ,o|5,0)+X>4,i|4,l> 
+2? 3 , 2 |3,2)+2? 2 , 3 |2,3) 
+Pi >4 |l,4)+D ,5|0,5) 

where the T> are expansion coefficients, which is of course equivalent to various 
first quantization expressions. But now we would say it is the two modes \U) 
and \D) that are entangled, not the 5 bosons! Entanglement for TV = 5 boson 
pure states is associated with there being six distinct Fock states that occur 
for five bosons being split between two modes. If there were four modes then 
for iV = 5 boson pure states there would be many more distinct Fock states 
available depending on how the bosons are divided amongst the modes. It is 
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more a question of combinatorics rather than symmetrization which is relevant 
in determining the dimension of the space of entangled states. A quite different 
picture of what is meant by an entangled state occurs when entanglement refers 
to modes rather than particles. 

Finally, the concept of particle entanglement is introduced somewhat differ- 
ently by Wiseman et al [T§] . 
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11 Appendix 3 - Reference Frames and Super- 
Selection Rules 

Several papers such as [37] , [3J5] , , [S] j [H] ; [33 , [33j explain the link between 
reference frames and super-selection rules (SSR). In this Appendix we present 
the key ideas involved. 

11.1 Two Observers and Reference Frames 

The first point to appreciate is that there are two observers - Alice and Charlie 

- who are involved in describing the same state of a particular quantum system. 
Charlie is the external observer, Alice the internal observer - perhaps closely 
linked to the system. It is important to realise that it is Charlie 's description of 
the quantum state which is of most interest, in particular how this description 
may differ from what Alice may regard as the system state. The system could 
be a multi-mode system involving identical particles, it could just be a single 
mode system or it could even be a single particle with or without spin. Alice and 
Charlie each describe quantum states in terms of their own reference frames, 
which might be a set of coordinate axes for the case of the spin or position 
states for the single particle system, or it could be a large quantum system with 
a well-defined reference phase in the case of multi-mode or single mode systems 
involving identical particles. Alice and Charlie may each choose from a set of 
possible reference frames - for the single particle case there are an infinite number 
of difference choices of coordinate axes for example, related to each other via 
rotations and/or translations. In Situation A - which is not associated with SSR 

- Alice and Charlie do know the relationship between their two reference frames 
(and can communicate this relationship via classical communications) - such as 
in the case of the single particle system when the relative orientation of their 
two different coordinate axes are known. In Situation B - which is associated 
with SSR - Alice and Charlie do not know the relationship between their two 
reference frames - such as in the multi-mode or single mode system involving 
identical particles when the relative phase between their two large quantum 
phase reference systems is not known. Alice and Charlie describe the same 
state via density operators a and /S, and the key question is the relationship 
between these two operators in situations A and B and for various types of 
reference frames. In terms of the notation in [26 p — > a and ~p — > /?. 

11.2 Symmetry Group 

A particular relationship going from Alice's to Charlie's reference frame is speci- 
fied by the parameter g, which in turn defines a unitary transformation operator 
T(g) that acts in the system space. Particular examples will be listed below. 
If there was a third observer - Donald - and the relationship going from Char- 
lie's to Donald's reference frame is specified by the parameter h, which in turn 
defines a unitary operator T(h), then if we symbolise the relationship going 
from Alice's to Donald's reference frame by the parameter hg, it follows that 
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T(hg) — T(h)T(g). This shows that the unitary operators satisfy one of the 
requirements to constitute a group, referred to generally as the transformation 
group. The other requirements are easily confirmed. The unitary operator 
T(0) = 1 corresponding to the case where no change of reference frame oc- 
curs (specified by the parameter 0) exists, and satisfies the requirement that 
f(0g) = f(0)f(ff) = f(g0) = T (g)f (0). The unitary operator T^T 1 ) = f ( ff )t 
corresponding to the relationship specified as g~ x that converts Charlie's ref- 
erence frame back to that of Alice exists, and satisfies the requirement that 
f (0) = f{g- 1 )T{g) = f(g)f(g- 1 ). Hence all the group properties are satisfied. 
A few examples are as follows: 

1. Translation group - single spinless particle system, with p, x.the momen- 
tum, position vector operators. Here is a vector giving the translation of 
Charlie's cartesian axes reference frame from that of Alice, thus g =.a^. The 

unitary translation operator is T(a^) = exp(ip- a^/h). 

2. Rotation group - single particle system, with J the angular momentum 
vector operators. Here u is a unit vector giving the axis and rotation angle 
4> for rotating Alice's cartesian axes reference frame into that of Charlie, thus 
g - ., (f>. The unitary rotation operator is T(_u<_, <ft) = exp(i</>J • u^/K). 

3. Particle number U(l) group - single mode bosonic system, with a the 
mode annihilation operator and N a = a 1 * a the mode number operator. Here 9 a 
is the phase change Alice's to Charlie's reference frame. The unitary operator 
isf(0„) =exp(iA> a # a ). 

4. Particle number U(l) group - multi-mode bosonic system, with a as a 
typical mode annihilation operator and N = ^ a) a the total number operator. 

a 

Here 9 is the phase change from Alice's to Charlie's reference frame. The unitary 
operator is T(0) = exp(iN9). 

In these examples the system operators p, J, N a , N etc are the generators 
of the respective groups. In many situations the generators commute with the 
Hamiltonian for the system (or more generally with the evolution operator that 
describes time evolution of the quantum state), in which case the group of 
unitary operators T(g) is the symmetry group, and the generators are conserved 
physical quantities. 

11.3 Relationships - Situation A 

In Situation A, where the relationship between the reference frames for Alice 
and Charlie is known and specified by a single parameter g, Alice's description 
of the state a is related to Charlie's description p for the same state via the 
unitary transformation 

p = T(g)dT(gy 1 (199) 

Note that this is a passive transformation - no change of state is involved, just 
the same state being described by two different observers. 
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As an example, consider the spinless particle and the translation group. If 



x. — 



is a position eigenstate then T( ap y = 4 — a^J. A pure quantum po- 
sition eigenstate described by Alice as a = |4>) ($| with state vector |$ 

would be described by Charlie as p = \ty) 1 but now with \ty 
which is also a pure quantum position eigenstate but with eigenvalue x — Qy 
This is as expected since Alices's cartesian axes have been translated by 4 to 
the origin of Charlie's axes without change of orientation. In the case of mo 
mentum eigenstates 



exp(i p ■ a^/K) 



4 
4 



so a pure 
would be 



we have T( 

quantum momentum eigenstate described by Alice with |3>) = 
described by Charlie with \ty) = exp(i p ■ a/h 

mentum eigenstate with the same eigenvalue p . Alice and Charlie describe the 

pure momentum eigernstate with the same density operator p = a, the phase 
factor cancels. 

For more general pure states, consider a quantum state described by Alice 



p j , which is also a pure mo- 



as a = 1$) ($1 with state vector 1$) = 



dx^ i 



14; 



States of this form 



can represent localised states when ^(4) is only significant in confined spatial 
regions, or they can represent delocalised states such as momentum eigenstates 

p ^ when 4>{x^) = (2irh)~ 3 / 2 exp(i p ■ x^/h). We see that Charlie also describes 



-4 W4- 



states |$) = / dp^ 



a pure quantum state p = |^) (*| but now with |^) = T(a^) |$) = J d 

4) | 4 ) = / ^ ^(4) so the wavefunction is now ^ = ^ + 4). 

Note that if Alice's state vector was written in terms of momentum eigen- 
(4) 4^, then Charlie's state vector | = / dp^ip{p^) 4 

has a momentum wave function ip(p) — exp(ip • 4 /ft) 4>(p) related to that 
of Alice by a phase factor. Note that a state which is a quantum superposition 
of momentum eigenstates as described by Alice is also described as a quantum 
superposition of momentum eigenstates by Charlie. A similar feature applies in 
all situation A cases, and is related to SSR not applying in situation A. 

The case of the particle with spin and the rotation group is outlined in Ref. 



11.4 Relationships - Situation B 

In Situation B, where on the other hand the relationship between frames is 
completely unknown, all possible transformations g must be given equal weight, 
and hence the relationship between Alice's and Charlie's description of the same 



S!) 



state becomes 

P = J w(g)dgf(g)df(g)- 1 

= G[a] (200) 

where J w(g)dg is a symbolic integral over the parameter g, which includes a 

weight factor w(g) so that J w{g)dg = 1. This linear process connecting a to p 

is the "Q- twirling" operation. Again, this is a paassive transformation. 
It is straightforward to show that for any fixed parameter h that 

f (typfih)- 1 = p (201) 

showing that Charlie's density operator is Q invariant under the transformation 
group - unlike the case for Situation A. 

As an example, consider the single mode bosonic system and the U(l) group. 
If \n a ) is a Fock state then T(9 a ) \n a ) — exp(in a a ) \n a ). Consider a pure quan- 
tum state described by Alice as the Glauber coherent state a = |$) ($| with state 
vector | $(/?)) =^ C(n a ,(3) \n a ), where C(n a , (3) = exp(-|/3| 2 /2) /y/fnJL 

It is straightforward to show that 

f(0 a ) |$(/3)> = |$(/3exp(ifl ))) (202) 

so that the Glauber coherent state is transformed into another Glauber coherent 
state, but with j3 changed via a phase factor to (3exp(i9 a ). The quantum state 
described by Charlie is given by 

P = |^|$(/3exp^ a )))(<D(/3exp^ a ))| (203) 

/flf) 
~£ EE C(n a ,/3)C(m a ,pyT(6 a ) \n a ) (m a | f(0 o )t 

-^exp(i[n Q -m a ]0 a ) 

= ^ |C(n a ,/J)| 2 |n a ) <n„| 

= ^ CX p(-|/3| 2 )i^: |n a )(n a | (204) 

which is a mixed state consisting of a Poisson distribution of Fock states with 
mean occupation number n a — \/3\ 2 . In view of the first expression for p it 
can also be thought of as a mixed state consisting of Glauber coherent states 
each with the same amplitude \(3\ = \/W^, but with all phases (arg/3 + 9 a ) 
equally probable. Thus, whereas Alice describes the state as a pure state that 
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is a quantum superposition of Fock states with differing occupancy numbers, 
Charlie describes the same state as a mixed state involving a statistical mixture 
of number states. The former violates the SSR whereas the latter does not. A 
similar feature applies in all situation B cases, and is related to SSR applying 
in Situation B. Whether Alice could ever prepare such a state in the first place 
is controversial - see the discussion presented above in SubSections 12.71 and 12.91 
However, assuming she could, the quantum state as described by Charlie is a 
mixed state. 

The situation just studied relates of course to the debate [55] regarding 
whether the quantum state for a single mode laser operating well above thresh- 
old should be described by a Glauber coherent state or as a Poisson statistical 
mixture of photon number states. The first viewpoint (Alice) describes the 
state from the point of view of an internal observer with a reference frame, 
the second (Charlie) describes the same state from the point of view of an ex- 
ternal observer for whose reference frame relationship to that of the internal 
observer is unknown. The debate is regarded by [3j5] as settled on the basis that 
both viewpoints are valid, they are just at cross purposes because they refer to 
descriptions of the same quantum state by two different observers. 

It should not be thought however that the quantum state would always be de- 
scribed in such a fundamentally different manner for all Situation B cases. As an 
example, consider the multi-mode bosonic system and the U(l) group. Consider 
the pure quantum state described by Alice as the multi-mode N boson Fock state 
a = |$) ($| with state vector |$(iV)) = \n 1 n 2 ...n a ...\ N) = Y\ ■■ \n a ) 

a 

where N = n a . We have T{9) |n 1 n 2 ...ri a ...; N) = e~xp(iN9) \nin 2 --.n a ...\ N), 

a 

so that the same state would be described by Charlie as p = |\&) (^| and with 
\^f) = \nin2...n a ...; N). This is also a multi-mode N boson Fock state with 
exactly the same occupancies. The product exp(iN9) exp(~iN9) of phase fac- 
tors averages out to unity and here p = a, so Alice and Charlie both describe 
the multi-mode Fock states in the same way. Another example for two mode 
bosonic systems and the U(l) group is provided by the one boson Bell states 
(the BS notation used here is non-conventional). These are entangled two mode 
states that Alice would describe via the state vectors 1$*) = (| 10) ± |01))/V2. 
We have T(9) l^) = exp(i9) l^), so that the same state would be described 
by Charlie with I**) = (| 10) ± |01))/V2. Again the product of phase factors 
averages to unity and p = a, so Alice and Charlie both describe the quantum 
states as Bell states, and in the same form. 

11.5 Dynamical and Measurement Considerations 

Discussions of the relationship between equations governing the dynamical be- 
haviour of Alice's and Charlie's density operators depend on whether the evolu- 
tion is just governed by a Hamiltonian or whether master equations describing 
evolution affected by interactions with an external environment are involved. 
Such matters will not be treated in detail here, nor will the issue of relating 
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Alice's and Charlie's measurements. The latter issue is dealt with in |26| . 

However, in the case where Alice describes the Hamiltonian evolution of her 
density operator via the Liouville - von-Neumann equation 

d 

ih—d=[H,d] (205) 

where in Alice's frame the Hamiltonian is H, and where in addition the trans- 
formation group is also the symmetry group so that T(g)HT(g)^ 1 = H for all 
g, it is easy to see that for both Situations A and B, Charlie's density operator 
will evolve via the same LVN equation 

ih^p=[H,p] (206) 

Thus both Alice and Charlie will describe the same dynamical evolution, though 
of course the initial (and hence evolved) states may differ in the two cases. 



11.6 Nature of Reference Frames 

Reference frames of differing types are involved for the various transformation 
groups. The common feature is that they are thought of as actual physical sys- 
tems themselves which are either macroscopic classical systems or macroscopic 
quantum systems in states associated with the classical limit. They are intended 
to be essentially unaffected by the presence of the systems for which they are 
acting as reference frames. In some cases relatively uncontroversial examples 
exist, such as for the cartesian axes associated with the translation and rotation 
groups associated with the single particle system. The physical reference system 
may be a large magnet whose magnetic field points in a well defined direction 
and defines a z axis, combined with an electrostatic generator whose electric 
field is in another well defined direction at right angles that defines an x axis. 
In other cases the existence of suitable reference frames is less clear. 

In this SubSection we will describe possible phase reference frames as if they 
are entirely separated (or uncorrelated) with the system of interest. In terms 
of the treatment by Bartlett et al [39], |26) these are non-implicated reference 
frames. In the next SubSection and in the next Appendix phase reference frames 
that are correlated with the system of interest will be described - these are the 
so-called implicated reference frames of Bartlett et al. 

For the large quantum system with a well-defined reference phase associated 
with the U(l) group in the case of multi-mode or single mode systems involving 
identical particles, the usual choice is a single mode bosonic system such as a 
single mode BEC or a laser with a large mean occupancy, and which is thought 
of as being prepared in a Glauber coherent state |3>(a)) in order to provide 
the phase reference frame, the reference phase being arg a. Whether such a 
reference frame really exists is controversial. The discussion presented above in 
SubSections 12.71 and 12.91 raises the question of whether such a phase reference 
state could ever be prepared, so this choice of a physical phase reference is rather 
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unsatisfactory. However, from the point of view of this presentation we assume 
it does, so that - as in the previous example - Alice can describe the reference 
state as another coherent state. Again, whether Alice could ever prepare such 
a state is questionable. 

Another possibility for a physical phase reference is a macroscopic low fre- 
quency harmonic oscillator, whose quantum energy eigenstates \n) - with n — 
0, 1, ..,n max and energies nhu> can be used to construct phase eigenstates \9 P ) 
with p = 0, 1, .., n max and 8 p — px 2tt/ (n max + 1), and which are defined by [70) 

1 n max 

\9 P ) = - exp(m6g \n) (207) 

These states are orthonormal. The separation between the equally spaced phase 
angles A9 = 2n/(n max + 1) can be made very small if n max is large enough. 
Under the effect of the harmonic oscillator Hamiltonian H = fruN, where N is 
the number operator, the phase state \9 p ) evolves into \9 p — a; At) during a time 
interval At, so if the time intervals are chosen so that uiAt = 27r/(n max + 1), the 
phase angle 6 P changes into 6 v —\. Thus the system behaves like a backwards 
running clock [79) . the phase angles 9 p defining the positions of the hands. If 
the clock initially has phase 9 P the probability of finding the clock to have phase 
9 q after a time interval At is given by 



1 sin 2 ((n max + l)A/2) 

(^max ~t~ 1 ) 2 sin (A/2) 



p(e q ,9 P , m) = - \ „r w :^:,T' ' <m 



where A = 9 p — 9 q — toAt. For times At such that toAt <C 27r/(n max + 1) the 
probability of the phase remaining as 9 p is close to unity. Thus if the phase state 
\9 p ) is used as a phase reference, it will remain stable for a time At satisfying 
the last inequality. For At ~ lOO^js and n max ~ 10 4 so that phase is defined to 
~ 10 -3 radians, an oscillator frequency uj ~ 10° s _1 would suffice for this phase 
reference standard. Such macroscopic oscillators do exist, though the process 
to prepare them in the phase reference quantum state \9 P ) would be technically 
difficult. Whether such a system would be useful as a phase reference for optical 
fields or a BEC is another issue 



11.7 Relational Description of Phase References 

In this SubSection phase reference frames that are correlated with the system 
of interest will be described - these are the so-called implicated reference frames 
of Bartlett et al [39]. [26]. 

One such approach to describing phase references in the U(l) group case 
is via the concept of maps. For simplicity consider a one mode system S, the 
basis vectors for which are Fock states \m) s , where it is sufficient to restrict 
m = 0, 1, .., rn max . The reference system R, will also be a one mode system 
with Fock states \n) R , where n is large. Product states \m) s <g> \n) R for the 
combined modes exist in the Hilbert space Hs ® Hr and are eigenstates of the 
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various number operators, including the total number operator Nt = N$ + Nr 
- where the eigenvalue is I = to + n. The product states may be listed via 
m — 0, 1, .., TOmax and n = 0, 1, ..or m = 0, 1, .., TO max and I = to, to + 1, .... 
Here we will describe how a coherent superpostion of number states, such as a 
Glauber coherent state can be represented. 

In the so-called internalisation or quantisation of the reference frame the 
product state \m) s <g> |ra) H is mapped onto the product state \m) s ®\n — m) R 
where n > m max . Thus 

\m) s <g> |ra) fl -> |m) s <g> |n - ra) fl (209) 
Hence for a linear combination of system states given by 

l*>s= E C -»Ws (210) 
we have for the state |$) s (g> |n) fl in Hs ® ifit 

"Irnax mmax 

|$) s ® |n) fl = ^ C ro |m) s ® |n) B -> ^ C ™ l m )s i n ~ m )fl = l*n>AS 

m— m—Q 

(211) 

The mapping results in an entangled state where there are n bosons distributed 
betweeen the two modes. This state \^ n )ns 1S a P ure state which is compatible 
with the SSR and is in one-one correspondence with the original system state 
|3>)s- Note that to create this state the reference state \n) R must have more 
bosons in it than TO max . The density operator for the original pure system 
S state would be as — \&) s (&\ s , and we note that this state violates the 
SSR. The state |<J>) g would be essentially a Glauber coherent state if C m = 
exp(— |a| 2 /2)a m /(v / TO!), with TO max ^> \a\ 2 . However, for the mapped state 
\^ n )ns t ne reduced density operator p s is given by 

Ps = Tr R (\* n ) RS (*\ RS ) 

rn max 

= J2 l^l 2 HsHs (212) 

m=0 

This is a mixed state and is compatible with the SSR. For the Glauber coherent 
state this is the Poisson distribution of number states. Hence the original 
SSR violating superposition of number states for system S is mapped onto 
a state in the combined system for which the reduced density operator is a 
statistical mixture and is consistent with the SSR. as would correspond to 
Alice's description of the state, p s to Charlie's. 

In the alternative so-called externalisation of the reference frame the map- 
ping is between product states, and is the reverse of the previous mapping. The 
product state \m) s (S> \n) R is mapped onto the product state \m) s (S> \m + n) R 
in the Hilbert space Hs <8> Hr where the former is spanned by vectors \m) s and 
the latter by vectors \m + n) R , and where n > m max . Thus 

\m) s (g> \n) R — > \m) s ® |m + n) R (213) 
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The mapping of the H$ Hr state \^ti}rs ^ nen * s 

m max 

l*«)fls = E C m \m) s ®\n-m) R 

m max /^max \ 

-> £ C m |m) s ® \n) R = I E C ™ \m) s ® |n) fl - |£ n )££L4) 

m=0 \m=0 / 

The mapping results in a non-entangled state which is incompatible with the 
SSR. The state in the subspace Hs is a coherent superposition of number states, 
whilst that in Hr is a Fock state. The reduced density operator in Hs is erf 
given by 

a* = Tr R (\~ n ) RS (~ n \ RS ) 

m max m max 

= E E C ™°1 \ m )s(k\s (215) 

m=0 fe=0 

which is the same as as — \&) s (&\s an d involves coherences between different 
number states in contradiction to the SSR. Clearly this second mapping just 
reverses the first one. 

Of these two treatments of phase reference frames, the internalisation ver- 
sion has a closer link to physics in that the pure state \^ n )ns can m prin- 
ciple be created and does lead to a way of creating a state that is in one- 
one correspondence with any SSR violating pure state | <1>) ^ , though it is in 
the form of an entangled state of the S, R sub-systems rather than just S 
alone. This is an important point to note - the original SSR violating state 
does not exist as a state of a separate system, all that exists is an SSR com- 
patible entangled state that is in one-one correspondence with it. However, 
the general process for creating a state such as \^ n }ns ls no * explained. For 
simple cases such as |<j>) s = (|0) s + |l) s )/\/2 the creation of the required state 

l^™)flS = (|0)s ® l n )fi + |l)s ® \ n ~~ l)i?)/v / 2 ; where n > 1 would seem feasible 
via the ejection of one boson from a BEC in a Fock state \n) R into a previously 
unoccupied mode. . 

11.8 Irreducible Matrix Representations and Super-selection 
Rules 

If \i) (i = 1,2,..) are a set of orthonormal basis vectors in the system state 
space, then the group of unitary operators T(g) is represented by a group of 
unitary matrices D(g) 

f(g)\i)=Y / D ji (g)\j) (216) 

3 

with elements Dji(g), and such that D(hg) — D(h)D(g) etc corresponding to 
the group properties of the operators. This is a matrix representation of the 
transformation group. 
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The theory of such group representations and their application to quan- 
tum systems is well established, following the pioneering work of Wigner in the 
1930s. We can just use the results here. A key concept is that of irreducible 
representations. Within the system state space we can in general choose so- 
called irreducible sub-spaces, denoted as T a of dimension d a and spanned by 
new orthonormal basis vectors |r Q A) (A = 1, 2, .., d a ) such that 

d a 

f(g) |r Q A) = £^ A ( 3 ) \T a p) (217) 

For each irreducible sub-space T a there is no smaller sub-space for which the 
operation of all T(g) just leads to linear combinations of vectors within that 
sub-space. The d a x d a matrices D a (g) then form an irreducible matrix repre- 
sentation for the transformation group. For different a the representations are 
said to be inequivalent. 

The irreducible matrices satisfy the so-called great orthogonality theorem [80j 

w{g)dgD« x {g)Dl T (g)* = ^-S aP S^S XT (218) 

Oa 

The proof of this result is based on Schur's lemma. 

The importance of the irreducible representations and the consequent or- 
thogonality theorem lies in its application to Situation B cases, where we have 
seen that Charlie's density operator p is invariant under any of the transforma- 
tions T(h)'pT(h)~ 1 = p. Suppose we represent p in terms of the basis vectors 
|r Q A) associated with the irreducible representations 

?=EE< |r„A)<I>T| (219) 

where R will be a Hermitian, positive definite matrix with unit trace since it 
represents a density operator. Applying the transformation gives 

T^pTih)- 1 = J2J2 R frD^(h)\r af i)(T^\Dl T (h)* 

= p (220) 
Averaging over h and using the great orthogonality theorem gives 

p = E E ( E y RTx ) |raAi> {Ta ^ (221) 

This is in the form of a mixed state involving irreducible state vectors |r a p) 
each occuring with a probability given by 
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which is the same for all /i associated with a given irreducible representation 
T a . This is clearly a positive real quantity and since 



A A 



= Trp = 1 (223) 

the probabilities sum to unity as required. 

The final result for Charlie's density operator 

P = EE pa \ T ^)(r a ^\ (224) 

demonstrates the presence of a super-selection rule. In Charlie's description 
of the quantum state there are no coherences between states |r Q /i) associated 
with differing irreducible representations of the transformation group. This 
represents the general form of the SSR for all transformation groups in Situation 
B cases. 

As an example, consider the U(l) group and the single mode bosonic system. 
Since the Fock states satisfy T(9 a ) \n a ) = exp(m a # a ) \n a ) they form the basis 
for the irreducible representations of the U(l) group, the occupation number 
n a specifying the irreducible representation and the lxl matrices exp(m a # a ) 
being the unitary matrices. Hence Charlie will describe the quantum state as 

p = ]TPK) \n a )(n a \ (225) 

Tla 

which is a statistical mixture of Fock states with no coherences between different 
Fock states. This result is of the same form as in Eq. (l4l))) and is in accord with 
the SSR on boson number. 

As another example, consider the U(l) group and the multi-mode bosonic 
system. Here sums of products of Fock states 

\nin 2 ...n a ... ■N) = nK>K> ■>«>... N = J2^ a (226) 

a a 

such that the total occupancy is N — n a can be used to form irreducible 

a 

representations for the transformation group in terms of linear combinations of 
the products with the same N. Writing these linear combinations as 

\K)= E C Sn,..» a ... } l™...n a ...;iV) (227) 

{n 1 n2-..n a ...} 

we have since T(6)\nin2--.n a ...;N) — exp(iN8) |niri2...n a ...; N) we see that 
f{6) \^%) = exp{iN9) |*^) also, so the |*^) define the irreducible basis states. 
The total occupancy N specifies the irreducible representation, but here there 
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are many irreducible representations with the same N depending on the various 
/i. In this case Charlie will describe the state as 



? = EE P ^ (228) 

N M 

which is a statistical mixture of multi-mode states all with the same total 

occupancy N. Although there are coherence terms between individual modal 
Fock states, there are no coherences between states with different total occu- 
pancy. This result is of the same form as in Eq. (|27|) and again is an example of 
a super-selection rule operating in terms of Charlie's description of the quantum 
state. 

Finally, we note that in situation A where the relationship between the 
frames is known and there is no invariance for Charlie's density operator, we do 
not have SSR applying. For the single particle case and the translation group 
the momentum states p \ define the irreducible representations, each specified 
by p , and as we saw Charlie's description of the quantum state involved linear 
combinations of these irreducible basis vectors, in contradiction to the SSR. 



11.9 Non-Entangled States 

The essential feature of an non- entangled or separable state is that the sub- 
systems are considered to be unrelated to each other. Hence, both for Alice 
and Charlie there will be separate reference frames for each sub-system, with 
transformation groups - TA(g a ) for sub-system A, Tb (<?;,) for sub-system B, etc 
which relate the reference systems of Alice to those of Charlie. The transfor- 
mations g a , Cfb, .. are different. The overall transformation operator would be 
of the form T(g a ,gi,, ..) = Ta((? ) ® ^b(s i 6) ® — Alice would describe a general 
non-entangled state as having a density operator 

a = ^P R a R ®a R ®d R ® ... (229) 

R 

It then follows for Situation B where the reference frames for Alice and Char- 
lie are unrelated, that Charlie would describe the same state via the density 
operator 

R 

where 

Pr = J w(g c )d 9c f c (g c ) d c R fc(gc)- 1 C = A,B, .. (231) 

Note that separate twirl operations are applied to the different sub-systems, as 
explicitly shown in the papers by Vaccaro et al [31] (see Section IIIA, Eqn. 3.3 
therein) and Paterek et al [33] (see Section 6). This leads for general trans- 
formation groups to the local group super- selection rule, where the p R involve 
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no coherences between states associated with differing irreducible representa- 
tions of the transformation group. We see that Charlie also describes a non- 
entangled state and with the same mixture probability Pr as for Alice. Thus 
non-entanglement or separability is a feature that is the same for both Alice 
and Charlie, as ought to be the case. 

In the context of sub-systems consisting of modes (or sets of modes) occupied 
by identical bosons, the case of interest is Situation B, with each transformation 
group being U(l). Here the relationship between Charlie's and Alice's phase 
reference frames arc unknown. Hence irrespective of Alice's description of the 
sub-system states a R , a R ,... we see from the previous section that Charlie will 
describe the separate sub-system states p R , as statistical mixtures of number 
states for the separate modes (or total number states for the sets of modes in 
each sub-system). Thus from Charlie's point of view the separate mode density 
operators will satisfy the SSR. Thus we see that the introduction of reference 
frames and two observers - Charlie being the external one whose description of 
the quantum states is of primary interest - leads to the same SSR outcome as 
the simpler considerations set out in SubSections 12.71 and 12.91 Essentially the 
same considerations have been used in |25j . [31] and the other papers to justify 
the local photon number supers election rule. 
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12 Appendix 4 - Super-Selection Rule Violations 

? 



12.1 Preparation of Coherent Superposition of an Atom 
and a Molecule ? 

A key paper dealing with the coherent superposition of an atom and a molecule 
is that by Dowling et al [52] , entitled "Observing a coherent superposition of an 
atom and a molecule" . Essentially the process involves one atom A interacting 
with a BEC of different atoms B leading to the creation of one molecule AB, 
with the BEC being depleted by one B atom. 

12.1.1 Hamiltonian 

The Hamiltonian is given by 

H = hw A b A b A + hw M b\fi M + huj 2 blb 2 + ^ (b\ I b A b 2 + b M b A b\) (232) 

where b A ,bu and b 2 are standard bosonic annihilation operators for the atom, 
molecule and BEC modes respectively, uj a ,ujm and co 2 are the corresponding 
mode frequencies and k defines the interaction strength for the process where a 
molecule is created or destroyed from/to an atom A and a BEC atom B. A is 
the frequency difference between the molecular state AB and the two separate 
states for atoms A and B - this is zero on Feshbach resonance - and is given by 

A = u>m — u A — lu 2 (233) 

The Hamiltonian commutes with the total number operator _/V tot , where 

N tot = 2 V M b M + b\b A + b% (234) 

where the molecule number operator is multipled by two. 

12.1.2 Initial State 

Initially the state of the system is given by the density operator Eqs (10) and 
(11) in the paper 

W QL = J ^exp(-iiV tot 0)|*) OL (*| OL exp(+iJv tot 0) (235) 
I*)ol = WW (236) 

where \A) is a state with one atom A and |/3) is a Glauber coherent state for 
the BEC of atoms B.The super-operator acting on the pure state \^) 0L (^Iql ^ s 
called the twirling operator, the group of unitary operators exp(— iN to tO) depend 
on a phase variable 9 and are a unitary representation of ?7(1), the generator 
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being N tot . These operators act as a symmetry group for the system and leave 
the Hamiltonian invariant. The initial state is also given by 



W L = 

Pa-m(°) = 



(237) 
(238) 




/ 27 ex PH" 2<? ) \@) C 8 ! exp(+m 2 0) 
^p„(<n>)|n) (n| 



(239) 



(240) 



n 




(241) 



where ni — b 2 b2 is the number uperator for the BEC mode and p n {< n >) — 
{exp(— < n >) < n >" /n\} is a Poisson distribution, whose mean is < n > = 
|/3| 2 . Initially then there is one atom A and the BEC is in a statistical mixture of 
number states with a Poisson distribution, which is mathematically equivalent 
to a statistical mixture of Glauber coherent states |/5exp(— i6)) with the same 
amplitude y/< n > but with all phases (arg (5 + 9) being equally weighted. 

12.1.3 Implicated Reference Frame 

In the paper by Dowling et al |62j the BEC is acting as an implicated phase 
reference frame (see [39], |26j). The state of the reference frame as described 
by Charlie is given by 



and from Eq. (|232|) , there is an interaction between the reference BEC and the 
separate atom A and molecule M systems. However, because < n >= |/3| 2 is 
very large, the BEC is essentially unchanged during the process, as reflected in 
the use of approximations in eqs (27), (28) of the paper. Another implicated 
phase reference frame situation, but involving a two mode reference frame is 
discussed in the paper by Paterek et al [33] 

Overall, in terms of the discussion in Appendix HT1 Wp r. would be Charlie's 
description of the initial state, whereas Alice would describe it as \^) 0L (^Iol- 
Presumably in the paper by Dowling et al [62] what is referred to as the " state 
of the laboratory" be Charlie's reference frame, and what they refer to as the 
"internal reference frame" would refer to that of Alice. However, whether Alice 
could actually prepare such a state as \^) 0L {^\ol * s controversial - see SubSec- 
tions 12.71 and 12. 9\ though here this is assumed to be possible. 

12.1.4 Process - Alice and Charlie Descriptions 

There are three stages in the process, the first being with the interaction that 
turns separate atoms A and B into the molecule AB turned on at Feshbach 



Pref = Pa(0) = 




— exp(-in 2 9) \P) (/3| exp(+m 2 0) 

Z7T 



(242) 
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resonance for a time t = n/ (2k < n >), the second being free evolution at large 
Feshbach detuning A for a time r leading to a phase factor <f> = At, the third 
being again with the interaction turned on at Feshbach resonance for a further 
time t = 7t/(2k < n >). The typical initial state \^) 0L given by \A) (eq 
(11)) evolves into |*) 3i given by (see eq. (32) of paper) 



Wsl = ( sin(|) \A) - exp(z arg/3) cos(|) \M) 1 |/3) (243) 

using approximations set out in eqs (27), (28) of the paper that depend on 
< n > being large. Here \M) is a state with one molecule AB. Thus it looks 
like a coherent superposition of an atom state \A) and a molecule state \M) 
has been prepared, the atom plus molecule system being disentangled from the 
BEC. Alice would describe the final state of the system as \^) 3L | 3i , so from 
her point of view a coherent superposition of an atom and a molecule has been 
prepared. 

However, for Charlie the final state of the system is described by a den- 
sity operator W^l which is reconstructed by applying the twirling operator to 

l*>3£ (*lsi • Notin g that 

exp(-iN tot 9)\^) 3L = (exp(-tfl) sm(~) \A) - exp(-2i6) exp(iarg/3) cos(^) \M)^j \(3exp{-i9)) 

(244) 

and using 

Tr 2 (\0exp(-i6)) (/3exp(-i0)|) = (/3 exp(-i9)\/3 exp(-i9)) = 1 (245) 

we see that Charlie's final reduced density operator for the atom-molecule system 
is 

Pa-m(3) = Tr 2 W 3L 

/d6 
— exp(-iN tot 9) |*) 3L (y\ 3L exp(+iN tot 6) 

= [ ( j— (exp(-i6) sin(^) \A) - exp(-2i6) exp(iarg/3) cos(~) \M) 
J 2n V 2 2 



x ^exp(+i6») sin(p (A\ - cxp(+2i6>) exp(-i arg/3) cos(^) (M| 
= sin 2 (^) |A)(A|+cos 2 (^) |M)(M| (246) 

Thus the coherence terms like \A) (M\ and \M) (A\ do not appear in the final 
density operator when the average over 9 (not 0) is carried out. 

For Charlie the density operator for the atom and molecule is of course a 
statistical mixture of a state with one atom and no molecule and a state with 
no atom and one molecule. The authors of [62] actually point this out in the 
paragraph after eq (35) where (presumably for the case 4> — 7r/4) it is stated "the 
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state is found to be . . . an incoherent mixture of an atom and a molecule." . The 
probabilities for detecting an atom A or a molecule AB are as in eq (33) of the 
paper. In terms of Charlie's description, the density operator at the end of the 
preparation process does not signify the existence of a coherent superposition of 
an atom and a molecule, as the title to the paper might be taken to imply. The 
existence of such a coherent superposition would of course be present in Alice's 
description, but it is Charlie's (laboratory) description that is more relevant. 

12.1.5 Coherence Effects Without SSR Violation 

Note that coherence effects are still present since the atom or molecule detection 
probabilities depend on the phase cj) associated with the free evolution stage of 
the process. However, as in many other instances, the presence of coherence ef- 
fects does not require the existence of coherent superposition states that violate 
the super-selection rule. The authors actually point this out in the paragraph 
after eq (35), where it is stated "we have clearly predicted the standard oper- 
ational signature of coherence, namely Ramsey type fringes, but the coherence 
is not present in our mathematical description of the system." What they are 
referring to is Charlie's description of the final state - which indeed shows no 
such coherence, but the belief that coherent superposition states are needed to 
predict coherence effects is mistaken. 

To drive this point home, the process can be treated with the initial state for 
the BEC being given as a Fock state \N). With the interaction being given as in 
Eq. (|232p (eq (14) in the paper) the state vector is a simple linear combination 
of two terms 

|¥(t)) = A(t) \A) \N) + B(t) \M) \N - 1) (247) 

This is of course an entangled state. Coupled equations for the two amplitudes 
A(t) and B(t) can easily be obtained and simple solutions obtained for stages 
where the Feshbach detuning is either zero or large. The state vector is continu- 
ous from one stage to the next , and the reduced density operator at the end of 
the three stage process for the atom plus molecule sub-system can be obtained. 
It is of the form 

?a-m(3) = Tra(|*(3)><*(3)|) 

= sin 2 (^) |A)(A|+cos 2 (|) \M)(M\ (248) 

which is of course a statistical mixture of a state with one atom and no molecule 
and a state with no atom and one molecule - and is exactly the same result as 
obtained in the paper by Dowling et al. [62] . Note that coherence effects in regard 
to the interferometric dependence on <fi for measurements on the final state has 
been found without invoking either the description of the BEC via Glauber 
coherent states or the presence of a coherent superposition of an atomic and a 
molecular state. The result can easily be extended for the case where the BEC 
is initially in a statistical mixture of Fock states with differing N occuring with 



103 



a probability Pn- Each initial state \A) \N) evolves as in Eq. (12471) . We then 
would have 

Pa-m(3) = Tr 2 {J2 p N |*iv(3)) (*jv(3)|) 

N 

= EPn (sin 2 (|) \A) (A\ + cos 2 (^) \M) (M|) 

= sin 2 (^) \A) (A\ + cos 2 (^) \M) (M\ (249) 

which is the same as before. Allowing for a statistical mixture of Fock states 
makes no difference to the interferometric result. 

12.1.6 Conclusion 

Dowling et al [62] state in their abstract that "we demonstrate that it is possible 
to perform a Ramsey-type interference experiment to exhibit a coherent super- 
position of a single atom and a diatomic molecule" . However the interferometric 
or coherence effects (involving the dependence on <p) cannot be said to exhibit 
the existence of such a coherent superposition, since the same interferometric 
results can be obtained without ever introducing such a quantum state. There is 
not a convincing case that quantum states that violate the super-selection rule 
forbidding the creation of coherent superpositions of Fock states with differing 
particle numbers can be created, even in Alice's reference system. The fact that 
an SSR violating state \^) 3L is created in Alice's reference system is not 

surprising, because in the process considered the initial state |/3) for the BEC 
was assumed as a factor in Alice's initial state, and this was itself inconsistent 
with the SSR. Furthermore, such SSR violating states are not needed to de- 
scribe coherence and interference effects, so that justification for their physical 
existence also fails. 

12.2 Detection of Coherent Superposition of a Vacuum 
and a One-Boson State ? 

Whether such super-selection rule violating states can be detected has also not 
been justified. For example, consider the state given by a superposition of 
a one boson state and the vacuum state (as discussed in 63J). Consider an 
interferometric process in which one mode A for a two mode BEC interferometer 
is initially in the state a |0) +/3 |1), and the other mode B is initially in the state 
|0) - thus |tf (»)) = (a \0)+/3 |1))a®|0) B in the usual occupancy number notation, 
where |a| 2 + |/3| 2 = 1. Modes A, B could refer to two different hyperfine states of 
a bosonic atom with non-relativistic energies ftw^.and hu)g, mode annihilation 
operators a, b. The modes are first coupled by a beam splitter, which could be 
a resonant microwave pulse that causes transitions between the two hyperfine 
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states and which can be described via a unitary operator Ubs such that 

Ubs(\1) a ®\0) b ) = (|l) A ®|0) B -i|0)^®|l) i3 )/V2 
Ubs(\0) a ® |1> B ) = H|1) A ®|0) B + |0) A ®|1) B )/^ 
C>B5(|0) A ® |0) B ) = (|0> A ®|0) B ). (250) 

After passing through the beam splitter the system is allowed to evolve freely 
for a time r, the Hamiltonian being Hf ree = (mc 2 + ^^4)0^0+ (mc 2 + fiiJB)b^b - 
where collisional effects have been ignored and the rest mass energy included for 
completeness. Following the free evolution stage, the modes are then coupled 
again via a beam splitter, and the probability of an atom being found in modes 
A, B then being measured. A straightforward treatment of the evolution shows 
that the final state is given by 

|*(/)> - a(|0) A ®|0) B ) 

+j3 exp(— i{mc 2 /h + uia}t) 

(1 — exp(— iAr) ., . . . . 1 + exp(— iAr) ,, , , . ,\ 
^ ® |o) B ) - * ^ -(|o) A ® |1) B )J 

(251) 

where A = lub — lo 'a is the detuning. The probabilities of finding one atom in 
modes A, B respectively are 

P 10 = \P\ 2 sin 2 (Ar/2) P 01 = |/3| 2 cos 2 (Ar/2) (252) 

Thus whilst coherence effects occur depending on the phase difference 4> = At 
associated with the interferometric process, the overall detection probabilities 
only depend on the initial state via |/3| 2 . There is no dependence on the rela- 
tive phase between a and /3, as would be required if the superposition state 
a|0) + j3\l) is to be specified from the measurement results. Exactly the 
same detection probabilities are obtained if the initial state is the mixed state 
p(i) = \a\ 2 (\0) A (0\ A ® |0) B (0| B ) + |/3| 2 (|1) A (1\ A ® |0) fl (0| B ), in which the vac- 
uum state for mode A occurs with a probability \a\ 2 and the one boson state 
for mode A occurs with a probability |/3| 2 . In this example the coherent super- 
position associated with the super-selection rule violating state would not be 
detected in the interferometric process. The paper by Dunningham et al [6 3) 
considers first a detection process that involves using a Glauber coherent state 
as one of the input states. Similar interference effects as in Eq. (|252j) are ob- 
tained. A second detection process in which the single term Glauber coherent 
state is replaced by a statistical mixture with all phases equally weighted in 
considered next, leading to the same interference effects. This again confirms 
that it is not necessary to invoke the existence of coherent superpositions of 
number states in order to demonstate interference effects. 
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13 Appendix 5 - Non-Physical Two Mode States 



We now consider some possible states for the second mode B - to be combined 
with , an d Pi, P2 as in Eq. (|42j) . These states are two general pure 
orthogonal states of the form a |0) B +/3 \1) B and — j3* \0) B + a* \l) B with (|a| 2 + 
|^| 2 ) = 1. We have 

pf = ((a\0) B +(3\l) B ))((a*(0\ B +(3*(l\ B )) 

P2 = ((-/8*|0) B + a*|l) B ))((-)8<0| B + a(l| fl )) (253) 

This gives the reduced density operator 

P B = ^(|0)b<0| b ) + 5(|1) b <1| b ) (254) 

A straightforward calculation gives for the overall density operator for the 
two mode non-entangled state as in Eq. ([TJ 

P = Pl+?2 

where the 'p 1 , /5 2 are contributions that are consistent with or inconsistent with 
the super-selection rule. We have 

Pi = \\0) a (0\ a ®\0) b (0\ b + ~\1) a (1\ a ®\0) b (0\b 
+ \ |0) A (OU ® 11} s (lis + \ |1>a <!U ® |1>b (11b 
+ \a*p \0) A (1U ® |1>b (0|b + ^/3* (OU ® |0>b (11b (255) 

and 

? 2 - J(M 2 - l/?l 2 ) |o) A (iU ® |o>b (o|b + J(M 2 - l/?l 2 ) (OU ® |0>b (01b 
+J(I/?I 2 - M 2 ) n A (iU ® |i)b (11b + J(l/?l 2 - M 2 ) (OU ® |1)b (11b 

+1<**P \1)a (°U ® |1>b (01b + \<*F \®) A (!U ® |0)b (11b (256) 

Now to make p 2 = requires |a| 2 = | /? | 2 so that the first four terms in ~p 2 
are zero. This in combination with \a\ 2 + \(3\ 2 ) = 1 leads to |a| = |/3| = 
However this results in the remaining two terms in ~p 2 ~ which are coherences 
between N = and N = 2 states - always being non-zero. Overall then, no 
choice of a, (3 will lead to a overall density operator which is physical. Adding 
further states |2) B , |3) B , ...does not rectify the problem. 
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14 Appendix 6 - Derivation of Sorensen et al 
Results 

Sorensen et al [41] derive a number of ineqalities from which they deduce a 
further inequality for the spin squeezing parameter in the case of a non-entangled 
state. From this result they conclude that spin squeezing implies entanglement. 
The final inequality they obtain for a non-entangled state is 

AS 2 z )>^(^S x y + (s y ) 2 ^ (257) 

Their approach is based on writing the density operator for a non-entangled 
state of N identical particles as in Eq. (fM|) 



P 



n 



The spin operators are defined as S x =J2^x = X«(l0bW) {4> a Wl+I^a(*)) (^hWI)/ 2 

i 

;S y = j:si = Ei(l&(0) (0.(i)l-l^«(O> (4>M)/u ;S z = J2Sl = Ei(l0 6 (i)> 

z z 

|</> a (i)) ((/> a (i)|)/2 , where the sum i is over the identical atoms and each atom is 
associated with two states \<j) a ) and |</> h ). Clearly, the spin operators satisfy the 
standard commutation rules for agular momentum operators. 

Sorensen et al [H] state that the variance for S z satisfies the result 

ASl) (si) 2 R + £ Pn (S z )l - (S Z ) 2 (259) 

To prove this we have 

i i 

2 



R 



TV 



where we have used 
, 2 1 



(5) = -^M^M-^ai^ium 2 



\{\<t>S)) (hiiMM W)\ - (I0 6 (i)> <&(i)l*.(i)> 

+ ~HI*«(0) (^a(i)l^(i)) + (ka«> <^a(*)l0a(*)) (0aWI) 

j((l^«X^«l + (l^(i)X^«l) 

il 4 (261) 
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a result based on the orthogonality, normalisation and completeness of the states 
|4,(i)>, |&(i)>. Also 

S z ) o = Tr(p R j:Sl) 

I H i 

- 

so eliminating the term r Pr . (^) (^*) ^ gives the required expres- 
sion for ( A5?\ = (s*\-(Ss 



Next, Sorensen et al [H] state that 

2 \ 2 



S x ) <NJ2P R J2(Si) (Sy) <Nj2PRE\(s l y )\ 2 (263) 

' R i ' R \ I j \ I R 

To prove this we have 

[S x ) = ^PflTrG^ESi) 



'$*)\ < e p «?K^)j ( 264 ) 

since the modulus of a sum is less than or equal to the sum of the moduli. Now 



S x ) — I \ S X 



* (Eft? I <*>,!) 



2 

using the general result that I ^ VC-R ) < E Pr Cr> where Pr : = 1 with 

V -R J R R 

here v / Cr = | ( ) I ■ Next consider 

ei(^)j) 2 = (ei(^; 

= Z(\(s l x)\-\(si)\) 2 > (266) 
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so that 

sS <nJ2PrE\(s 1 x ) J 2 (Sy) 2 <nJ2PrJ2\(s;)\ 2 (267) 

R i ^ ' R x ' ^ i x/ " 

which is the required result. The inequality for (Sy) i s proved similarly. 
Another inequality is stated [4T for (S z \ . This is 



s z ) 2 <E P « (s*) 2 R ( 268 ) 



R 



To show this we have 



S z ) = J2 P « Tr (PR^ 

R. 

R 1 

= E p *(^ 

R 

s,)\ < E p «i(^)J ( 269 ) 



so that 



< E p «i(^ 



using the general result that ( Y Pr VCr ) < Yl Pr Cr, where Y Pr — 1 with 
here = I (s a 



R 

Finally, we hnd that 



(271) 



109 



To show this we use the properties of the density operator p l R for the ith particle 
of Hermitiancy, positiveness, unit trace Tr(jj l R ) = 1 and Tr(p R ) 2 < 1. In terms 
of matrix elements of the density operator p R between the two states \(j> a (i)), 

\<f>l,(i)) the quantities (^S^ , (^Sy) and (sij are 

si) R = rr^id^cox^WI + I^COX^WI)) 

S\) r = l(pU-pL) (272) 

where p l cd = (^ c (i)|Pii l0d(O)- The Hermitiancy and positiveness of ~p R show 
that p 4 bb and p l aa are real and positive, p l ah = (pjj* and p l oa p£ b - \p l ab \ 2 > 0. The 
condition Tr{jj l R ) — 1 leads to p aa + p\ h = 1, from which Tr(p R ) 2 < 1 follows 
using the previous positivity results. Taken together these conditions lead to 
the following useful parametrisation of the density matrix elements 

Pla = sin2 a i Pbb = COs2 °i 



p l ab — y sin ai cos 2 a, sin ft exp(+i0 i ) p\ a — y sin cos 2 a* sin ft exp(— i<f>^) 

(273) 

where a^, ft and <^ are real. In terms of these quantities we then have 

~"\ 1 2 

Sxf — -sm2ai s'm /3 i cos(j> i 

= ^ sin 2«j sin 2 sin ^ 

Si) = \ cos 2a, (274) 
It is then easy to show that 

Si) 2 R + (si) R + (St)l = J-Jsin 2 2a* (1- sin 4 ft) 

< \ (275) 

and the final inequality (12711) then follows by taking the sum over particles i 
and then using ^2 r Pr = 1- If only the Schwarz inequality is used instead of 
the more detailed consequences of Hermtiancy, positiveness etc it can be shown 

that /SM +(^) + (^) -J' wllidl though correct is not useful. 

Combining the inequalities in Eqs. (J263]), PSgfl and (|27Tj) into Eq. (f259| 
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shows that 



N / — \ 2 

>- T -$*5(si) R 



> ^((^) 2 + (s,) 2 ) ( 276 ) 

for the case of a non-entangled state. This result is that in Sorensen et al.|41j. 



Ill 



A ^ 2 ) - E p « ( a ^ 2 ) ( 278 ) 



15 Appendix 7 - Revised Sorensen et al 

15.1 Variance ^AS 2 X 

Here we will see if the modified approach to Sorensen et al can lead to a useful 
inequality for (AS%\ or.^AS 2 ) that applies when non-entangled states are 

those when all the separate modes a* and hi are the sub-systems . We will 
attempt to follow the approach used for the simple two mode case in Section 2J 
Firstly, the variance for a Hermitian operator Q in a mixed state 

P = Y,PrPr (277) 

R 

is always greater than or equal to the the average of the variances for the separate 
components 

> A \ " P. / \02^ 

R R 

where (AO 2 ) = Tr(pAtt 2 ) with AO = S7-(ft) and^Aft 2 )^ = Tr(p R AQ R 2 ) 

with Att R = O — (fj) . The proof is straight-forward and given in Ref. [69] . 

Next we calculate (AS 2 ) , (AS 2 ) and (s x j , (s y ) , (s z ^ for the 
case where 

p = E p n {pb 1 ®pr)® {pr ® #0 ® ® ) ® • ( 279 ) 

as is required for a general non-entangled state a/Z 27V modes. Furthermore, 
the density operators for the individual modes must represent possible physical 
states for such modes, so the super-selection rule for atom number applies and 
we have 

((ai) n ) ai = =0 ((a\r) =Tr(T R i (alr) = 

\ I at 

(ft)*"), = Tr(p^(6 i ) m ) = ($)"») =Tr(^n = 

(280) 

The Schwinger spin operators are 





= 

i 


+ atS,-)/2 = 


E^ 

i 


Sy 




-a&)/2i = 


= 22 s l 

i 


A 


= E(^ 


-2toi)/2 = 





(281) 
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where 2j, bi and a], &j respectively are mode annihilation, creation operators. 
From Eqs. (|28T1) we find that 

Sl=E(S l x ) 2 + ES l x Sl (282) 

so that on taking the trace with ~p R and using Eqs. (|279[) we get after applying 
the commutation rules [e, — 1 (e = a or b) 

As we also have 



using Eqs. (|279|) and we see finally that the variance (^AS 2 ,^ is 

(ASl) R = E((Siy) R ~E(si) 2 R (285) 
all the terms with i ^ j cancelling out. and therefore from Eq. (|278[) 

[ASl) >J2 p nE (((St) 2 ) - (286) 



But using pSUj) 



(S z x ) 2 = j (b]a$la t + b%a% + a^tya, + afaafa) 
& )2 )r = ^(((6%) fl + <(« t S),) fl ) + ^(((S t a) J ) fl ((6%) i? )(287) 



and 



Sl) R = (288) 



so that 



(289) 

Now using (12"50|) 
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S z ) = ^PnEiSi 

R * 
R i 

< E p «iSi(((^%) fl -((« t «),)«))i 

R 1 

< E p * i +<(^> R )) ( 291 ) 



and thus 



A.s':')-||(.v 



R 1 

= E p «^?(((« t «) i >«((6% 



> (292) 

A similar proof shows that (^AS 2 ^ — g| I > f° r the non-entangled state 

of all 27V modes. 

This shows that for the general non-entangled state with all modes a* and bi 



as the sub-systems, the variances for two of the spin fluctuations (^AS 2 J and 
AS 2 ) are both greater than ^ | ( S z ) | , and hence there is no spin squeezing for 



S x or S y . Note that as | (S y ) \ = 0, the quantity J l\ (S ±1 ) \ 2 + \ (S ±2 )\ 2 



is the same as | z/ {■> so the alternative criterion in Eq. is the same as 
that in Eq. ff5U|) which is used here. 

For the other spin fluctuation (AS 2 ) since we have 



S x ) = J2 Pr (s*) r = {Sy)=J2 p R{Sy) R = ° ( 293 ) 

R. R 



then the other two uncertainty relationships just give \ AS y j (AS Z J > 0; 

AS Z 2 ^ (^AS X 2 ^ > 0, so spin squeezing in S z is meaningless. 
Hence we have shown that for a non- entangled physical state for all the 2N 



modes at and h 

f AS 2 \ > _, 

2' \ "/ 1 \ V ~ 2 



2\^ 1 |/3\l „„J / A C2\ - 1 



AS X )>-\{S Z )\ and (AS 2 )>-\(S Z )\ (294) 
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so that spin squeezing in either S x or S y requires entanglement. 



15.2 Variance ( AS 



Here we will see if the modified approach to Sorensen et al can lead to a useful 
inequality for ( AS* 2 \ that applies when non-entangled states are those when 



the pairs of modes and bi are the separate sub-systems . We will attempt 
to follow the approach used by Sorensen et al when identical particles i were 
regarded as the sub-systems. 

Now the general non-entangled state will be 



P 



= Y, p kpr®pr®pr®- ( 295 ) 



where the ~p R are now of the form given in Eq. (|157[) and the conditions in Eq. 
(|280[) no longer apply. The Fock states are of the form \Ni a ) (g> \Ni b ) for the pair 
of modes cii and bi , and fior this Fock state the total occupancy of the pair of 
modes is Ni = N ia + No,. From the super-selection rule the density operator p^, 
for the ith pair of modes and bi is diagonal in the total occupancy. For Ni 
= there is one non zero matrix element ((0| ia (g) {0\ ib ) Pr (|0) ia ® |0) ib ). For Ni 
— 1 there are four non zero matrix elements, which may be written 



a a 



«i| te ®<oi tt )? fl (ii> ia ®io> tt ) = P i 

(<ili„®<ou?k(|o) io ®|i) i6 ) = pU 
«oL®<iU)?i(|i> io ®|o) i6 ) - P \ a 

«0L®<iy^(|0) ia ®|l) i6 ) = p\ b (296) 

ForiVj = 2 there are nine non zero matrix element ((2| ia i8>(0| if) )^j (|2) ia ®|0) i6 ), 

— ) ((°L ® ( 2 Ub) Pr (|0)m ® l 2 )ib) an d the number increases with Ni. 

If we restrict ourselves to general entangled states where Ni = 1 for all pairs 
of modes, then the density operator ~p R is of then form 

Pr = PL(|l)i (lL®|0) i6 (OU) + p l a6 (|l) io <0| ia ®|0) i6 (iy 

+Pla(\0) ia (l\ ia ® \l) lh (0| i6 ) + 4 6 (|0) io (0| M ® \l) ib (297) 

In addition Hermitiancy, positivity, unit trace Tr(p R ) = 1 and Tr(p R ) 2 < 1 can 
be used as in Eq (|273l) to parameterise the matrix elements in (j296[) . 

.9 j 2 

p bb = COS Cti 



p\ b = \J sin 2 (Xi cos 2 a; sin 2 /3j exp(+i0 i ) = y sin 2 cti cos 2 a* sin 2 /3, exp(— i^J 

(298) 
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The expectation values for the spin operators S % x , S l y and S\ associated with 
the ith pair of modes are then 

(SI) R = Tr(p l R l -(pfi l+ Z%) 

= 2 (P*b + Pba) 

(si) R = \(pU-fL) (299) 

which are of exactly the same form as in Eq. (|272l) as in the Appendix [TJ] 
derivation of the original Sorensen et al [41] results based on treating identical 
particles as the sub-systems. The proof however is now different and rests on 
restricting the states p l R to each containing exactly one boson. 

The remainder of the proof is exactly the same as in Appendix [T?] and we 
find that 

(&Sl)>±((s x ) 2 + (s y ) 2 ^) (300) 

for non-entangled pairs of modes Sj and 6,;. Thus when the interpretation is 
changed so that are the separate sub-systems are these pairs of modes, it follows 
that spin squeezing requires entanglement of all the mode pairs. 
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16 Appendix 8 - Heisenberg Uncertainty Prin- 
ciple Results 



Here we derive the results in SubSection 15.41 leading to inequalities for the vari- 
ance (Aj|) considered as a function of | \J Z ) \ for states where the spin op- 



erators are chosen such that (^J X J = (Jy/ = 0- 

From the Schwarz inequality ( J z ) < \J 2 Z ) so that 



J >) + \ J v) + \ J >) <(Jl) + (Jl) + (Ji) = J(J+i) (301) 

giving Eq. (I169p . Subtracting (jx) — (j = from each side gives 

A,/*) + (AJ2) + (j z ) 2 < J(J+1) (302) 

Substituting for ^AJ^\ from the Heisenberg uncertainty principle result in 
Eq. ([rfD]) gives 

AJ^ 2 _ ^j(j + i) _ ( A J|) + y (J z ) 2 < (303) 

The left side is a parabolic function of /AJ|^ and for this to be negative 
requires ( A J | \ to lie between the two roots of this function, giving 



aji) > ±\(j(j+i)-(j z ) 2 )-J(j(j+i)-(j z ) 2 ) -e(j ; x " 



(304) 



1 / \2\ // ,^\2\ z 1^ i J 



AJ v - 2 U J(J+1) ~\ Jz / j+v(/< J+1 >-vw J ~ € \ J 

(305) 

which are the required inequalities in Eq. (|171|) and (|172|) . 
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17 Figure Captions 



Figure 1. Bloch vector and spin fluctuations shown for original spin operators. 

Figure 2. Regions in the < AJ^r > versus | < J z > | plane (shown shaded) 
for states that satisfy (a) the spin squeezing inequality Eq. (|173|) (b) the smaller 
Heisenberg uncertainty principle inequality Eq. (llTlf) and (c) the larger HUP 
inequality Eq. (|172l) . The case shown is for J = 1000 and HUP factor £ = 1. 
Both < Aj| > and | < J z > | are in units of J. The spin operators are chosen 
so that < J x >—< J y >= 0. 

Figure 3. As in Figure 2, but with J = 1000 and HUP factor £ = 10.0. 

Figure 3. As in Figure 2, but with J = 1 and HUP factor £ = 10.0. 
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Figure 3: 
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This figure "SpinSqgFigl.jpg" is available in "jpg" format from: 



http://arxiv.org/ps/1305.0788vl 



This figure "SpinSqgFig2.jpg" is available in "jpg" format from: 
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